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ABSTRACT 


The  research  was  undertaken  to  answer  issues  related  to  cavity  expansion  and  in 
situ  testing  using  the  pressuremeter.  Emphasis  of  the  research  was  placed  upon  the 
fundamental  concq)ts  underlying  the  cylindrical  cavity  expansion  theory. 

Tests  were  performed  in  a  cuboidal  shear  device  (CSD)  at  strain  rates  of  0.01  %, 
0.05%,  0.10%,  0.50%,  1.00%  and  5.00%  per  minute  on  artificially  sedimented  kaolin 
clay  and  kaolin-silica  mixture  whose  plasticity  indices  range  similar  to  the  most  of  the 
clays.  It  was  found  from  the  experiments  that  the  normalized  shear  strength  (with 
respect  to  0.01  %/min)  increases  linearly  with  the  logarithm  of  strain  rate.  The  increase 
in  undrained  shear  strength  in  the  pressuremeter  stress  path  is  about  14.3%  per  log  cycle 
for  kaolin  clay  and  15.3%  for  kaolin-silica  mixture.  The  undrained  shear  strength  in  the 
conventional  triaxial  test  was  found  to  increase  about  8  to  10%  for  a  tenfold  increase  in 
strain  rate.  Therefore,  it  can  be  concluded  that  the  undrained  shear  strength  increases 
about  40-50%  more  in  pressuremeter  stress  path  tests  than  in  the  triaxial  stress  path  tests. 

A  new  state  variable,  ip,  indicative  of  the  fabric  of  clays  is  introduced.  Based  on 
these  concepts  and  general  failure  criterion  a  simple  model  to  predict  failure  parameters 
of  anisotropic  clays  for  many  commonly  encountered  stress  paUis  is  developed.  The 
model  capability  to  interpret  in  situ  strength  measured  under  a  given  stress  path  and 
transfer  it  to  another  stress  path  is  illustrated.  Finally,  the  ability  to  obtain  failure 
parameters  for  any  stress  path  using  data  from  a  single  CIUC  test  is  demonstrated. 

The  model  developed  based  on  cavity  expansion  theory  is  able  to  incorporate  the 
influence  of  decreasing  strain  rate  along  the  surrounding  soil  mass.  The  findings  from 


the  strain  rate  test  for  the  pressuremeter  was  used  in  the  model  to  estimate  the  difference 
in  the  interpretation  of  the  imdrained  shear  strength.  Even  though  the  actual  undrained 
shear  strength  increases  with  the  strain  rate,  the  interpreted  undrained  shear  strength  from 
the  PMT  expansion  curve  shows  a  higher  strain  softoiing  behavior  for  the  higher  value 
obtained  for  PMT  stress  path.  From  the  parametric  studies  it  was  determined  that  the 
level  of  upper  yield  does  not  have  any  significant  influence  on  the  predicted  strength. 

From  the  PMT  simulated  in  the  CSD  it  was  found  that  for  kaolin  clay, 
irrespective  of  the  initial  strain  level  or  strain  rate,  the  relaxation  time  was  about  200  to 
250  minutes  for  undrained  and  drained  conditions.  Theoretical  models  have  been 
included  to  show  how  to  obtain  creep/relaxation  parameters  from  PMT  creep/relaxation 
tests  conducted  in  a  borehole.  Several  correlations  are  also  identified  in  order  to  obtain 
relaxation  parameters  from  creep  parameters  and  empirically  from  plasticity  index. 
These  creep/relaxation  parameters  are  required  to  predict  the  time  dependent  behavior 
of  soil. 


Model  pressuremeter  expansion  in  calibration  chamber  (CC)  is  numerically 
simulated  using  Finite  Element  Method  (FEM).  The  elasticity  based  FEM  model  was 
developed  using  commercially  available  program  ANSYS.  The  elasticity  based  analysis 
was  helpful  to  further  investigate  the  boundary  effects  in  chamber  testing.  Numerical 
simulation  of  pressuremeter  testing  using  critical  state  finite  element  program  CRISP  will 
be  performed  and  its  validity  will  also  be  evaluated. 

The  test  data  from  the  CSD  and  CC  are  used  to  train  a  three  layer  feed  forward 
neural  network  using  a  back  propagation  algorithm.  The  potential  of  artificial  neural 
network  will  then  be  evaluated  for  incorporating  strain  rate  effects  in  pressuremeter 
testing  and  for  similar  applications  in  the  area  of  computational  mechanics. 


INTRODUCTION 


The  geotechnical  oigineering  profession  today  is  increasingly  using  in  situ  devices 
to  determine  soil  parameters.  In  situ  testing  is  intended  to  reduce  disturbance  and  stress- 
relief  associated  with  sampling  and  laboratory  testing.  A  major  drawback  of  most  in'  situ 
devices  is  that  the  boundary  conditions  are  poorly  defined,  as  a  result,  the  test  data  are 
usually  interpreted  empirically  to  correlate  the  soil  properties.  Of  all  the  in  situ  tests 
now  in  use,  the  pressuremeter  offers  the  greatest  possibility  to  interpret  the  test  results 
analytically  because  of  its  well  defined  boundary  conditions.  Furthermore,  load 
deformation  relation  is  directly  obtained  from  the  test,  and  this  information  may  with 
proper  interpretation  yield  the  constitutive  relationship  for  the  soil  and,  in  situ  horizontal 
stress,  deformability  modulus  and  in  some  cases,  the  consolidation  characteristics  of  the 
soil. 


However,  in  spite  of  the  tremmdous  potratial  of  the  pressuremeter,  it  has 
encountered  several  problems.  One  of  the  major  problem  was  that  the  pressuremeter  test 
(PMT)  predicts  higher  undrained  shear  strength  and  modulus  values.  There  are  also 
inconsistencies  in  the  results  obtained  by  diffoent  user,  mainly  because  of  differences 
in  rate  of  probe  expansion,  probe  advancing  techniques,  and  interpretatitni  methods. 
Unfortunately,  there  are  no  formal  standards  for  performing  pressurem^  tests. 

One  of  the  major  diffoence  between  the  PMT  and  laboratory  tests,  for  ocample, 
triaxial  test,  is  the  rate  of  straining  the  soil  specimen.  Other  differences  betweoi  the 
PMT  and  laboratory  tests  include  the  stress-relief  and  stress  relaxation,  which  is  not  a 
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usual  concern  in  the  standard  laboratory  testing.  The  main  target  of  the  research 
undertaken  is  to  study  how  the  time  dependent  factors  affect  pressuremeter  tests. 
Emphasis  is  givoi  to  quantifying  the  strain  rate  effect  in  the  PMT,  and  studying  the 
influence  of  stress  relief  and  stress  relaxatitm  in  the  PMT  by  simulating  pressuremeter 
stress  path  in  cuboidal  shear  device,  as  well  as  conducting  model  pressuremeter  test  in 
calibration  chamber. 

Chapter  2  describes  computer  controlled  calibration  chamber  equipment  to 
perform  model  pressuremeter  tests.  A  control  system  having  electro-pneumatic  control 
and  differential  pressure  transducer  has  been  developed  to  maintain  the  same  internal  and 
external  cell  pressures  throughout  the  K,  consolidation  phase.  This  assures  a  rigid 
system  in  the  horizontal  direction.  Other  developments  such  as  automatic  data 
acquisition  and  reduction,  improved  slurry  preparation  technique,  and  an  effective 
piezometer  design  are  also  described  among  other  details  including  some  experimental 
results.  The  calibration  chamber  syyparatus  is  used  to  perform  model  pressuremeter  tests. 

A  detail  description  of  the  cuboidal  shear  device,  slurry  preparation  technique, 
properties  of  the  soils  used,  servo-controlled  system  to  perform  K,  consolidation  and 
strain  controlled  tests  are  provided  in  Ch^ter  3.  Simulation  of  pressuremeter  test  strain 
path  in  the  cuboidal  shear  device  is  discussed  in  Chapter  4.  The  main  thrust  of  this 
chapter  is  to  compare  various  strain  rate  related  tests  with  the  strain  rate  tests  in 
pressuremeter  test.  Major  conclusions  are  made  from  the  extensive  experimental 
program.  A  numerical  scheme  is  developed  to  incorporate  the  strain  rate  effect  in 
pressuremeter  test  and  the  ^rifluence  of  several  parameters  that  affect  the  shear  strength 
are  also  studied  in  Chapter  S. 

Chapter  6  discusses  the  effect  of  stress  relief  and  stress  relaxation  in 
pressuremeter  test.  Tests  were  conducted  in  cuboidal  shear  device  to  study  the  stress 
relaxation  and  normalization  period.  Theoretical  models  have  been  included  to  show  how 


to  obtain  creep/relaxation  parameters  from  PMT  creep/relaxation  tests  conducted  in  a 
borehole. 

The  main  conclusions  of  the  research  are  re-evaluated  in  the  final  chapter. 
Appendix  A  contains  a  listing  of  the  technical  papers  and  discussions  already  published 
or  in  preparation  on  the  results  of  this  research. 


MODEL  PRESSUREMETER  TESTING 
USING 

AN  AUTOMATED  FLEXIBLE  WALL  CALIBRATION  CHAMBER 


2.1  Introduction 

An  automated  control  and  data  logging  system  has  been  developed  to  perform 
model  pressuremeter  tests  in  clays  using  a  flexible  wall  calibration  chamber.  Several  of 
the  important  questions  related  to  cavity  expansion  and  pressuremeter  testing  in  clays  are 
being  investigated  using  the  newly  developed  calibration  chamber  test  set  up.  Methods 
for  slurry  preparation  and  consolidation  under  K„  conditions  are  presented.  Special 
considerations  have  been  given  to  prqxare  high  quality  specimens  and  perform  pore 
pressure  measurements.  Preliminary  test  results  indicate  successful  performance  of  the 
system. 


A  series  of  one  eighth  scale  model  pressuremeter  tests  are  being  conducted  to 
evaluate  the  effects  of  strain  rate  and  stress  disturbance  in  cavity  expansion  of  cohesive 
soils.  The  initial  phase  of  the  research  involved  the  development  of  an  automated  control 
and  data  logging  system  to  perform  model  pressuremeter  tests  in  a  double  wall 
calibration  chamber.  The  calibration  chamber  built  for  a  previous  research  study 
(Huang,  1986)  was  modified  to  work  with  the  new  set  up. 

The  test  set  up  includes  a  double  wall  calibration  chamber,  a  slurry 
consolidometer,  differential  piston  pump,  model  pressuremeter,  piezometer,  electro- 
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pneumatic  control,  data  acquisition  system  and  control  panels.  The  flexible  wall 
calibration  chamber  system  and  the  procedure  involved  in  preparing  uniform  specimens 
obtained  by  consolidating  a  slurry  of  kaolin  or  kaolin-silica  mix,  are  described.  The 
performance  of  newly  designed  piezometers  is  also  discussed.  Using  the  same  soils,  tests 
were  performed  in  a  cuboidal  shear  device  to  simulate  pressuremeter  stress  paths,  with 
strain  rates  varying  from  0.01  %  per  minute  to  5.0%  per  minute.  The  experimental  data 
indicated  an  increase  in  undrained  strength  of  15%  for  each  tenfold  increase  in  strain  rate 
(Skandarajah  et  al.,  1991).  One  of  the  main  objectives  of  this  research  is  to  confirm  the 
strain  rate  effects  measured  in  the  true  triaxial  simulation.  Calibration  chamber  test 
results  will  also  be  used  to  calibrate  the  anisotropic  model  developed  for  the  study. 

2.2  Strain  Controlled  Pressuremeter  Tests  in  a  Calibration  Chamber 

Most  of  the  scaled  in  situ  tests  using  a  calibration  chamber  has  been  performed 
in  granular  soils.  In  the  past,  pressuremeter  tests  in  calibration  chamber  were  performed 
in  granular  soils  by  Jewell  et  al.,  1980  and  in  cohesive  soils  by  Huang  et  al.,  1988.  A 
calibration  chamber  provides  a  controlled  environment  to  perform  this  type  of  testing 
with  the  unique  capability  of  subjecting  soil  samples  to  known  stress  history  and 
boundary  conditions.  It  also  makes  it  possible  to  work  with  uniform  and  reproducible 
samples. 

The  calibration  chamber  approach  for  pressuremeter  testing  in  cohesive  soils  has 
been  used  by  Huang,  1986  prior  to  this  research.  Other  studies  concentrated  field 
pressuremeter  tests  and  comparison  with  conventional  laboratory  tests  on  samples  from 
the  same  site.  It  is  believed  that  calibration  chamber  technique  for  cohesive  soils  is  a 
desirable  alternative. 

The  two  types  of  calibration  chambers  used  in  practice  are  the  rigid-wall  and 
flexible  wall  chambers.  A  rigid-wall  chamber  has  a  rigid  wall  that  ensures  no  lateral 
strain.  The  main  disadvantage  in  a  rigid-wall  chamber  is  that,  to  avoid  boundary  effects, 


the  size  of  the  chamber  must  be  much  larger  than  the  in-situ  device  that  is  being 
calibrated.  This  makes  the  testing  expensive  and  time  consuming.  In  a  flexible  wall 
chamber,  it  is  possible  to  have  independent  control  on  vertical  and  lateral  stresses,  which 
makes  the  simulation  of  field  tests  possible  with  relatively  small  samples.  The  flexible 
wall  chamber  used  in  the  current  research  has  two  walls,  the  inner  wall  being  slightly 
thinner  than  the  outer  one  so  that  if  the  cell  pressure  at  some  stage  exceeds  its  yield 
value,  it  would  burst  inward  ensuring  the  safety  of  the  working  environment.  The  two 
types  of  boundary  conditions  used  are  BCl:  Constant  vertical  stress  and  zero  lateral 
strain,  and  BC2:  Constant  vertical  stress  and  constant  lateral  stress.  It  can  accommodate 
a  specimen  of  size  200  mm  in  diameter  and  367  mm  in  height.  The  design  allows  for 
Ko  consolidation  as  well  as  pressuremeter  testing  on  that  specimen  at  the  end  of 
consolidation.  In  the  flexible  wall  chamber,  the  sample  is  hydraulically  confined  around 
a  soft  membrane  which  makes  it  possible  to  have  independent  control  on  lateral  and 
vertical  stresses.  In  order  to  make  the  consolidation  and  testing  time  manageable,  the 
specimen  dimensions  were  scaled  down  to  the  above  mentioned  values.  The  following 
sections  give  a  detailed  account  of  the  salient  features  of  the  test  set  up,  techniques  of 
sample  preparation,  chamber  consolidation,  and  strain  controlled  model  pressuremeter 
testing. 

2.2.1  Chamber  Top  Platen 

The  chamber  top  platen  is  shown  in  Fig.  2.1.  It  is  a  28  mm  thick  aluminum 
platen  and  200  mm  in  diameter.  The  model  pressuremeter  is  attached  to  it  at  its  center 
through  1/4"  NPT.  Piezometers  are  also  attached  to  this  platen.  They  are  arranged 
diagonally  opposite  at  20  mm,  27  mm,  and  44  mm  respectively.  It  is  thus  possible  to 
obtain  two  porepressure  readings  at  the  same  distance  from  the  center  of  the  probe.  The 
chamber  top  platen  becomes  the  bottom  platen  of  the  slurry  consolidometer.  The  tips  of 
the  piezometers  extend  approximately  to  the  center  of  the  probe.  The  chamber  platen 
provides  connectors  for  back  pressuring  and  also  for  initial  flushing.  The  pore  and  probe 


Fig.  2.1 


Isometric  View  of  Chamber  Top  Platen 
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pressures  are  measured  using  high  sensitivity  Model  AB  type  pressure  transducers  (Data 
Instruments,  Inc.)  and  the  chamber  top  platen  has  all  the  mounting  adapters  attached  to 
it.  Prior  to  using  the  platen  in  the  slurry  consolidometer,  all  the  piezometers  and  model 
pressuremeter  are  filled  with  deaired  distilled  water  and  thoroughly  flushed  to  remove 
any  entrapped  air  bubbles.  The  model  pressuremeter  has  a  custom  made  9.5  mm  ID,  0.8 
mm  thick  latex  membrane.  Since  light  has  detrimental  effect  on  latex,  when  not  in  use, 
the  platen  is  stored  in  a  dark  place. 

2.2.2  Piezometers 

Evaluating  pore  pressures  is  an  important  and  integral  part  of  calibration  chamber 
test  and  is  essential  to  the  success  of  undrained  testing.  It  is  important  to  saturate  the 
piezometers  completely  so  that  negative  pore  pressures  resulting  from  cavitation  are 
avoided.  However,  excess  positive  pore  pressure  should  also  be  avoided.  The 
piezometer  shown  in  Fig.  2.2  was  designed  based  on  these  considertions.  The 
piezometer  is  typically  made  of  19  gauge  stainless  steel  hypodermic  needle.  The  needle 
is  attached  to  the  1/8"  pipe  union  through  brass  solder.  The  total  length  is  160  mm  and 
the  tip  is  plugged  with  fine  polyethylene  porous  plastic  having  a  mean  pore  size  of  10  to 
20  microns.  At  the  end  of  each  test,  porous  plastic  plugs  are  removed  using  0.56  mm 
plain  steel  acoustic  string.  The  small  diameter  of  the  pressure  sensitive  area  makes  it 
ideal  for  instantaneous  response.  Initial  test  results  indicate  their  acceptable  performance 
and  authors  are  considering  to  improve  their  performance  by  using  silicone  oil  instead 
of  water  for  initial  saturation  of  the  piezometers. 

2.2.3  Slurry  Tank 

The  slurry  preparation  procedure  must  be  systematically  followed  to  obtain 
uniform  and  reproducible  specimens.  The  slurry  preparation  involves  three  phases: 
slurry  mixing,  slurry  vacuuming,  and  transfer  to  consolidometer.  A  tank  shown  in  Fig. 

2.3  has  been  designed  to  perform  the  these  tasks.  It  is  343  mm  in  diameter  and  559  mm 
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Cross-Section  of  Piezometer 


Cross-Section  of  Slurry  Tank 


in  height.  Known  amounts  of  water  and  kaolin-silica  mix  are  added  to  the  tank  and  are 
thoroughly  mixed  using  conventional  heavy  duty  mixer  with  1/3  H.P.  motor  and  type  316 
stainless  steel  blades,  at  1725  rpm,  until  a  uniform  mix  is  obtained.  Subsequently,  the 
mix  is  subjected  to  full  vacuum  for  a  period  of  six  hours  by  which  time  the  air  entrapped 
during  mixing  is  removed.  A  pressure  of  5  psi  is  applied  to  the  surface  of  the  mix  to 
transfer  the  slurry  slowly  into  the  consolidometer  through  25.4  mm  tygone  flexible 
tubing.  The  mix  is  placed  gradually  from  bottom  to  top  in  the  consolidometer  and  the 
chances  of  air  entrapment  are  minimized.  Since  very  viscous  mixing  is  involved,  it  is 
important  to  properly  select  the  suitable  motor  and  blades  for  obtaining  a  uniform  slurry. 

2.2.4  Slurry  Consolidometer 

The  soils  used  for  testing  to  date  are  Georgia  Kaolin  and  crushed  Ottawa  silica. 
Deionized  and  deaired  water  added  to  the  soil  and  uniform  slurry  mix  is  obtained  using 
the  technique  mentioned  above.  The  water  content  corresponds  to  twice  the  liquid  limit 
of  the  soil.  The  slurry  is  then  transferred  to  the  slurry  consolidometer  shown  in  Fig. 
2.4.  The  device  is  used  to  consolidate  the  slurry  from  an  initial  height  of  800  mm  to  a 
final  height  of  approximately  350  mm.  The  slurry  consolidometer  has  2  steel  pipe 
compartments  having  an  inside  diameter  of  200  mm.  The  chamber  top  platen  becomes 
the  bottom  platen  for  this  device,  with  all  the  piezometers  and  model  pressuremeter 
attached  to  it.  A  porous  stone  is  attached  to  the  piston  and  the  slurry  is  consolidated  by 
207  kPa  air  pressure  applied  to  the  piston.  Double  drainage  is  allowed  for  the  slurry  to 
consolidate  under  K^  conditions.  It  takes  approximately  12  days  for  the  com.pletion  of 
primary  consolidation.  The  lower  compartment  has  same  height  as  the  soil  sample.  It 
is  split  longitudinally  in  two  halves  and  bolted  together.  The  interior  is  lined  with 
sandpaper  which  is  required  to  prevent  slippage  of  the  membrane  caused  by  adhesion 
between  clay  and  the  membrane  during  the  process  of  consolidation.  The  additional 
space  for  the  slurry  during  the  initial  phase  of  the  consolidation  is  provided  by  the  upper 
compartment  which  is  bolted  to  the  lower  one.  The  0.64  mm  thick  custom  made  latex 
membrane  for  the  specimen  extends  out  of  the  lower  compartment  and  provides  a  seal 
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between  the  two  cylinders  (Fig.  2.4).  This  completes  the  first  phase  of  the  consolidation 
process.  The  sample  confined  in  the  lower  compartment  and  encased  in  the  latex 
membrane  is  then  transferred  to  the  calibration  chamber  with  the  split  mold.  Since  the 
sample  is  transported  without  touching  it,  and  no  mechanical  extrusion  is  involved, 
disturbance  is  essentially  eliminated. 

2.2.5  Differential  Piston  Pump 

A  differential  piston  pump  is  used  to  expand  the  probe  at  a  constant  rate  of  strain 
(Fig.  2.5).  It  has  two  pistons  of  diameters  12.738  mm  and  9.525  mm,  respectively. 
They  are  attached  to  stepper  motor  and  a  DCDT.  Through  the  stepper  motor  control 
box,  four  different  speeds  are  preset  with  a  switch  box.  By  controlling  the  rate  of 
movement  of  the  stepper  motor  and  with  the  input  valve  closed,  the  volume  of  water 
injected  to  the  probe  can  be  regulated  at  a  constant  rate.  Thus,  with  the  preset  switch 
box,  it  is  possible  to  obtain  four  constant  strain  rates  for  the  probe  expansion.  A  device 
of  this  accuracy  is  necessary  because  a  radial  strain  of  12%  for  the  model  pressuremeter 
corresponds  to  a  differential  piston  movement  of  only  37  mm. 

2.2.6  Electro-pneumatic  Control 

During  the  second  phase  of  consolidation  inside  the  calibration  chamber  it  is 
essential  to  create  a  rigid  system  with  no  lateral  deformation  in  the  horizontal  plane.  To 
ensure  this,  the  cell  and  wall  pressures  are  balanced  throughout  the  consolidation  process. 
A  differential  pressure  transducer  monitors  the  pressure  difference  between  cell  and  wall, 
sends  the  resulting  voltage  signal  to  an  electro-pneumatic  transducer  which  converts  the 
input  voltage  to  an  appropriate  output  pressure  which  is  applied  to  the  annular  space 
between  the  internal  and  external  walls.  This  electro-pneumatic  control  assures  a  zero 
lateral  deformation  condition,  which  is  essential  for  successful  Kg  consolidation. 
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2.2.7  Flexible  Wall  Calibration  Chamber 

The  schematic  of  the  flexible  wall  calibration  chamber  is  shown  in  Fig.  6  (after 
Huang,  1986).  Upon  completion  of  the  slurry  consolidation  in  the  consolidometer,  the 
sample  confined  in  the  lower  compartment  (bolted  split  mold)  is  transferred  to  the 
chamber  and  placed  upside  down.  The  bottom  platen  of  the  slurry  consolidometer 
becomes  the  chamber  top  platen.  During  the  transfer  operations  it  is  critical  to  center 
the  sample  on  the  piston  of  the  calibration  chamber  before  placing  it.  Once  it  is 
supported  on  the  piston  platen,  it  is  very  difficult  to  re-align  it  without  disturbing  the 
sample.  The  piston  providing  the  vertical  stress  inside  the  chamber  has  the  same 
diameter  as  that  of  the  sample  and  has  a  maximum  vertical  travel  of  62  mm.  Since  the 
sample  is  encased  in  a  membrane  and'is  attached  to  top  and  bottom  platens  with  two  158 
mm  O'rings,  both  ends  of  the  sample  are  isolated  from  the  cell  water.  Hence  the  stresses 
in  horizontal  direction  through  cell  pressure,  and  vertical  direction  through  piston 
pressure,  can  be  independently  controlled.  By  balancing  the  pressure  between  cell  and 
wall,  and  by  maintaining  constant  piston  pressure  and  preventing  volume  change  in  the 
cell-water  system,  conditions  simulating  K„  consolidation  are  obtained.  The  sample  at 
the  end  of  the  first  phase  of  the  consolidation  process  in  the  slurry  consolidometer  does 
not  have  uniform  water  contents  along  the  vertical  direction  due  to  friction  between  clay 
and  rigid  wall.  By  subjecting  this  sample  to  a  second  phase  of  K„  consolidation  and 
avoiding  rigid  boundary,  a  very  uniform  sample  is  obtained.  Fig.  7  shows  the  three 
dimensional  profile  of  the  variation  of  water  contents  for  a  typical  kaolin-silica  mix 
specimen.  Porous  discs  are  placed  at  both  ends  of  the  sample  to  permit  double  drainage. 
During  consolidation,  loss  of  sample  volume  is  replaced  by  the  piston  movement, 
therefore  allowing  displacements  in  the  vertical  direction  only.  K„  consolidation  under 
an  effective  vertical  stress  of  276  kPa  was  performed  using  a  single  load  increment 
method  proposed  by  Campanella  and  Vaid,  1972. 

In  summary,  the  five  step  procedure  that  accomplishes  chamber  K„  consolidation 
is  made  of  the  following  steps; 


Fig.  2.6 


Cross-Section  of  Calibration  Chamber 
(After  Huang,  1986) 


1)  Close  the  pore  water  drainage  lines. 

2)  Increase  cell  and  axial  pressures  to  back  pressure  +  effective 

consolidation  pressure  (690  kPa  +  276  Id’a)  simultaneously. 

3)  Maintain  a  constant  cell-water  system. 

4)  Open  the  drainage  lines  and  permit  drainage  against  high  back  pressures. 

5)  Monitor  axial  deformation  and  cell  pressure  changes,  and  electro- 

pneumatically  control  the  wall  pressure. 

Because  of  the  limited  vertical  movement  of  the  piston,  very  soft  samples  could  not 
be  tested  using  the  present  calibration  chamber. 

2.3  Boundary  Effects 


The  sample  is  eighteen  times  larger  than  the  diameter  of  the  pressuremeter  probe. 
In  order  to  evaluate  the  radial  distance  at  which  stress  increase  diminishes  to  zero  a 
simple  elasticity  based  approach  can  be  used.  The  equation  of  equilibrium  for  a 
cylindrical  cavity  is  (tension  being  positive): 


da  a,~a 
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Introducing  compatibility  and  plane  strain  conditions,  the  variation  of  radial  and 
circumferential  stresses  are: 
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For  a  maximum  probe  pressure  of  1 173  kPa  (e,  =  12%)  and  a  cell  pressure  of  830  kPa, 
elasticity  based  analyses  indicate  that  the  stress  increase  is  essentially  negligible  at  a 


radial  distance  equivalent  to  9  times  the  probe  radius  (Fig.  8).  Numerical  analysis  by 
Carter  et  al.  (1979)  have  shown  that  in  a  soil  mass  which  extends  laterally  to  infinity,  the 
stress  increase  diminishes  to  essentially  zero  at  approximately  20  times  the  cavity  radius 
as  radial  strain  reaches  about  25%.  In  the  chamber  pressuremeter  tests,  maximum  radial 
strain  was  limited  to  12%.  Elasticity  based  arguments  coupled  with  the  numerical 
analysis  by  Carter  et  al.,  1979,  thus  indicate  that  the  size  of  the  probe  to  soil  specimen 
ratio  is  satisfactory  within  the  strain  rate  of  interest. 

2.4  The  Model  Pressuremeters 

Based  on  the  design  concept  of  a  single  cell  lateral  load  tester  made  by  Oyo 
Corporation,  Tokyo,  Japan  (Suyama  et  al.,  1982),  one  eighth  scaled  water  inflated  model 
pressuremeters  were  built  as  a  part  of  an  earlier  research  program.  The  probe  has  a 
diameter  of  11,1  mm  and  is  112  mm  long.  Custom  made  0.8  mm  thick,  9.5  mm  I.D. 
latex  membranes  are  used  as  the  probe  membrane.  The  fluid  pressure  is  measured  at  the 
center  of  probe,  accounting  for  the  head  loss.  Strain  controlled  pressuremeter  tests  are 
performed  by  injecting  constant  amount  of  fluid,  using  the  differential  piston  pump. 
Details  of  the  model  pressuremeter  are  shown  in  Fig.  9  (after  Huang,  1986). 

2.5  Strain-Controlled  Testing 

The  complete  control  system  involved  in  the  strain-controlled  probe  expansion  is 
shown  in  Fig.  10.  Initial  calibration  chamber  tests  using  the  model  pressuremeter  have 
been  performed  in  a  50/50  blend  of  kaolin  and  very  fine  crushed  silica.  The  slurry  was 
Ko  consolidated  under  a  pressure  of  207  kPa  around  the  model  pressuremeter  and  six 
piezometer  needles,  and  then  transferred  to  the  calibration  chamber.  The  performance 
of  the  piezometers  was  evaluated  during  the  B  parameter  check.  After  saturation,  the 
specimen  was  subjected  to  a  second  phase  Ko  consolidation  in  the  chamber  under  an 
effective  vertical  stress  of  276  kpa  against  a  back  pressure  of  690  kPa.  A  strain 
controlled  pressuremeter  test  was  performed  at  the  end  of  chamber  consolidation.  All 
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the  control  and  acquisition  during  the  testing  was  performed  by  a  Keithly  series  500  data 
acquisition  system.  The  probe  expansion  curve  obtained  at  the  end  of  the  test  represents 
the  combined  effect  of  system  compliance,  membrane  stiffness  and  the  soil  resistance. 
The  system  compliance  was  measured  by  inserting  the  model  pressuremeter  in  a  tight- 
fitting  thick  walled  brass  tube.  The  membrane  stiffness  was  evaluated  by  performing  the 
pressuremeter  test  in  an  empty  chamber  with  cell  pressure  equal  to  horizontal  stress 
obtained  at  the  end  of  chamber  consolidation.  The  net  pressuremeter  curve  is  obtained 
from  probe  expansion  curve  by  subtracting  the  system  compliance  and  the  membrane 
stiffness. 

The  results  indicate  that  small  radial  strains  occurred  before  exceeding  the  true 
soil  lateral  stress  which  is  the  chamber  cell  pressure  at  the  end  of  consolidation. 
Hence  with  reasonable  accuracy,  true  "lift-off  point  can  be  evaluated  using  the  model 
pressuremeter  under  controlled  conditions.  Essentially  elastic  behavior  was  observed  for 
the  unload-reload  cycle. 

A  typical  model  pressuremeter  expansion  curve  for  kaolin-silica  mix  is  shown  in 
Fig.  2.11.  The  compliance  curve  and  the  membrane  stiffness  are  also  shown  in  the  same 
figure.  The  variation  of  probe  pressure  with  respect  to  radial  strain, reported  in  Fig. 
2.11,  was  not  very  smooth  and  has  step  response.  This  was  found  to  be  due  to  the 
limitation  on  resolution  of  A/D  converter  in  the  data  acquisition  system  used.  Hardware 
gain  was  added  to  the  existing  system  for  subsequent  testing  and  much  better  response 
was  observed  for  kaolin  samples  (Fig.  2.12).  The  variation  of  pore  pressure  for  kaolin- 
silica  mix  during  chamber  consolidation  is  shown  in  Fig.  2.13.  The  excess  pore 
pressure  devclooed  during  pressuremeter  expansion  is  as  shown  in  Fig.  2.14  for  kaolin 
specimens. 

From  the  net  pressuremeter  curve  obtained  using  the  raw  data  for  kaolin-silica 
mix  specimen,  shown  in  Fig.  2.11,  ■  horizontal  stress  was  interpreted  as  820  kPa  and 

the  true  horizontal  stress  was  observed  to  be  830  kPa.  From  Fig.  2.12,  the  lateral  stress 
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Fig.  2.12 


Pressuremcter  Expansion  Curve  -  Test  5  (Kaolin) 
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Pore  Pressure  Dissipation  During  Chamber  K„  Consolidation 
(Test  5,  Kaolin) 


27 


Fig.  2.14 


Excess  Pore  Pressure  Developed  During  Model  Pressuremeter 
Expansion  (Test  5,  Kaolin) 
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was  interpreted  for  the  kaolin  specimen  to  be  815  kPa  and  the  true  lateral  stress  (cell 
pressure  after  K,,  consolidation)  was  820  kPa.  Hence  excellent  lateral  stress  predictions 
were  obtained  using  a  model  pressuremeter  test  in  the  calibration*  chamber  which  further 
validate  the  immense  potential  of  pressuremeter  to  predict  true  lateral  stresses  when 
performed  with  no  disturbance.  The  pore  pressure  dissipation  observed  in  two 
piezometers  while  consolidating  the  kaolin  specimen  under  Ko  cosolidation  is  shown  in 
Fig.  2.13  and  the  variation  of  pore  pressure  observed  on  two  different  piezometers  is 
quite  similar  and  validates  their  fiincitoning.  Excess  pore  pressures  were  observed  during 
the  probe  expansion  for  kaolin  sample  (Fig.  2. 14)  and  conclusive  results  on  pore  pressure 
behavior  can  made  after  additional  calibration  chamber  testing.  This  pore  pressure  data 
is  very  important  for  validating  many  of  the  existing  numerical  models  and  good  quality 
test  data  in  the  related  area  is  not  available. 

2,6  Summary 

An  automated  calibration  chamber  system  which  includes  a  double  wall  chamber, 
slurry  consolidometer,  piezometers,  electro-pneumatic  control,  model  pressuremeter,  data 
acquisition  and  control  system  has  been  developed.  From  the  testing  done  so  far,  the 
following  are  concluded; 

a)  By  employing  standard  slurry  preparation  techniques  and  two  phase 
consolidation  process,  a  very  uniform  sample  can  be  obtained. 

b)  True  lateral  stress  predictions  can  be  made  using  pressuremeter.  This  also 
reflects  the  disturbance  free  state  obtained  at  the  end  of  chamber  consolidation. 

Currently  model  pressuremeter  testing  using  the  calibration  chamber  is  being 
performed.  Conclusive  results  on  several  aspects  related  to  pressuremeter  testing  will 
be  made  after  performing  statistically  admissible  number  of  tests.  Since  only  couple  of 
tests  at  the  maximum  can  be  performed  for  each  month,  more  time  is  required  to  perform 
additional  testing.  Detailed  account  of  results  and  conclusions  based  on  calibration 
chamber  testing  will  thus  be  reported  at  a  later  date. 


LABORATORY  TESTING  USING  CUBOIDAL  SHEAR  DEVICE 


3.1  Introduction 

By  simulating  pressuremeter  stress  path  tests  in  the  cuboidal  shear  device*  the 
strain  rate  effects,  stress  relaxation  time,  oversize  probe  and  disturbance  effects  in  the 
pressuremeter  test  were  studied. 

A  slurry  consolidometer  was  used  to  prepare  102  mm  cubical  samples  from 
artificially  sedimented  mix.  The  cubical  specimens  were  reconsolidated  one 
dimensionally  in  the  cuboidal  shear  device  to  simulate  the  field  condition  and  then  to 
perform  pressuremeter  test  with  various  strain  rates.  The  measurements  were  made 
automatically  using  a  data  acquisition  system.  This  chapter  describes  the  slurry 
consolidometer,  cuboidal  shear  device,  measuremoit  systems,  membrane  prqiaration, 
properties  of  the  soils  used,  and  experimoital  procedures  for  various  tests  p^ormed. 
The  experimental  results  of  the  strain  rate  effects,  error  in  pressuremeter  test 
interpretation  due  to  strain  rate,  and  stress  relaxation  effects  are  presented  and  discussed 
in  Chapters  4,  S  and  6,  respectively. 

3.2  Slurry  Consolidometer 

A  slurry  consolidometer  made  of  plexiglass  was  used  to  prq)are  102  mm  cubical 
specimens  by  sedimentation  and  consolidaticm  undo*  K,  condition  from  the  powdered 
soil-water  mix.  An  isometric  view  of  the  slurry  consolidometer  is  depicted  in  Fig.  3.1. 
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Fig.  3.1 


Isometric  View  of  Slurry  Consolidometer  (Modified  from  Sivakug, 
1987) 


The  cross-sections  of  the  upper  and  lower  chambers  are  shown  in  Figs.  3.2  and  3.3, 
respectively. 

The  consolidometer  consists  of  two  chambers,  a  piston  and  a  base.  The  internal 
cross-section  of  the  chambers  is  a  102  mm  x  102  mm  square.  The  heights  of  the  upper 
and  lower  chambers  are  229  mm  and  102  mm,  respectively.  The  walls  of  the  lower 
chamber  are  externally  reinforced  in  order  to  minimize  lateral  deformation  and  to  assure 
that  the  sample  obtained  was  consolidated  one  dimensionally  under  K,  condition.  The 
interior  surface  of  the  lower  chamba  is  lined  with  teflon  in  order  to  minimize  the  wall 
friction  in  the  sedimentation,  consolidation  and  extrusion  stages.  The  consolidation  load 
is  applied  through  a  plexiglass  piston  (Fig.  3.4)  consisting  of  a  32  mm  diameter  292  mm 
long  rod  and  a  32  x  102  x  102  mm  base.  A  6.35  mm  thick  sintered  bronze  porous  stone 
is  attached  to  the  bottom  of  the  piston  to  facilitate  drainage  through  the  top  of  the  sample. 
There  are  16  holes  3.2  mm  in  diameter  connecting  the  porous  stone  to  the  top  of  the 
piston.  An  RTV  silicone  rubber  seal  is  attached  to  the  middle  of  the  piston  to  avoid  the 
slurry  from  being  squeezed  up  between  the  piston  and  the  chamber  walls  during  the  early 
stages  of  the  consolidation.  Another  identical  porous  stone  is  embedded  in  the  base  of 
the  consolidometer  (Fig.  3.5)  to  allow  bottom  drainage.  Provision  of  drainage  at  the 
bottom  and  top  reduces  the  drainage  path  by  half  and  thus  reduce  the  consolidation  time 
by  four  times. 

3.3  Cuboidal  Shear  Device 

Under  field  conditions,  the  soil  is  subjected  to  three-dimensional  states  of  stress. 
Hence,  for  appropriate  simulation  of  the  field  conditions  and  for  predictions  based  upon 
sophisticated  constitutive  models,  it  is  desirable  to  test  soil  specimens  under  truly  triaxial 
states  of  stress.  The  cuboidal  shear  device  is  such  an  apparatus  where  three  stresses  can 
be  applied  independently  under  any  stress  path. 
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Fig.  3.4 


Cross-Section  of  Piston  (After  Sivakugan,  1987) 


The  boundary  stresses  can  be  applied  to  a  cuboidal-shaped  specimen  either  via 
rigid  flat  platens  or  flexible  membranes,  or  a  combination  of  both.  The  first  attempt  to 
test  sand  in  a  cubical  triaxial  apparatus  by  Kjellman  (1936)  and  then  by  Jakobson  (1957) 
using  a  segmental  rigid  platen  had  limited  applicability  due  to  basic  mechanical 
difficulties.  Later  on,  the  mechanical  problems  were  overcome  and  the  non  uniform 
stresses  induced  in  the  rigid  platen  loading  were  minimized  by  using  layers  of  silicone 
grease  (Green,  1971),  and  flexible  rubber  membranes  on  each  of  the  three  pairs  of  faces 
for  tests  on  sand  (Ko  and  Scott,  1967).  Green  (1971)  and  Sture  and  Desai  (1979)  have 
described  in  detail  the  subsequent  developments  of  flexible  boundaries  and  the  various 
combinations  of  flexible  and  rigid,  and  ligid-lubricated  boundaries  on  two  pairs  of  faces. 
A  comparison  of  the  merits  and  drawbacks  of  the  three  types  o/ boundary  conditions  used 
in  true  triaxial  apparatus  is  given  in  Table  3.1  (after  Sture  and  Desai,  1979). 

Loads/pressures  applied  through  rigid  boundaries  assure  strain  controlled 
conditions,  and  uniform  strains  which  can  be  measured  precisely,  at  the  expense  of  non 
uniform  stresses  induced  in  the  sample.  The  major  issue  with  this  type  of  equipment  is 
the  interference  of  loading  platens.  Several  creative  techniques  were  proposed  (Pearce, 
1971)  to  overcome  this  problem,  however  they  consist  of  complicated  and  expensive 
mechanical  systems. 

The  flexible  membranes  can  be  used  for  stress  controlled  tests  where  uniform 
stress  distribution  is  possible,  however  the  uniformity  of  large  strains  is  difficult  to 
maintain.  The  interference  of  flexible  membranes  can  be  avoided  with  proper 
precautions  (Sture  and  Desai,  1979).  Plain  strain  or  any  other  strain  controlled  test  can 
be  performed  with  flexible  boundary  loadings  through  stepwise  corrections  using  a  trial 
and  error  approach. 
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The  mixed  boundary  apparatuses  avoid  the  boundary  interference  but 
heterogeneous  stress  and  strain  fields  occur  near  the  boundaries.  Usually  they  are  large 
devices  with  complicated  test  operations.  However,  there  is  not  any  one  type  of 
apparatus  most  suitable  for  testing  all  types  of  soil  over  a  large  range  of  stress  and  strain 
levels  and  paths. 

The  true  triaxial  apparatus  used  for  this  research  is  of  the  flexible  boundary  type, 
based  on  the  design  by  Sture  and  Desai  (1979).  A  detailed  description  of  the  device  was 
given  by  Sivakugan  (1987)  and  Sivakugan  et  al.  (1988).  An  isometric  view  of  the  CSD 
is  shown  in  Fig.  3.6.  The  dimensions  and  a  cross  sectional  view  of  the  space  frame  are 
shown  in  Figs.  3.7  and  3.8,  respectively. 

The  space  frame  and  the  cylindrical  casings  of  the  cuboidal  shear  device  were 
machined  from  solid  forged  billet  of  aluminum.  A  102  mm  cubical  sample  "floats”  in 
the  1 14  mm  cubical  cavity  of  the  space  frame  and  is  confined  within  six  identical  silicone 
rubber  membranes.  The  specimen  is  loaded  by  compressed  air  or  compressed  nitrogen 
applied  to  the  membranes  through  the  cylindrical  casings.  A  thin  coating  of  silicone  oil 
is  always  applied  on  the  surfaces  of  the  membrane  and  space  frame  in  order  to  minimize 
the  friction  between  the  sample  and  the  membrane,  and  between  the  membrane  and  the 
space  frame.  In  the  absence  of  friction,  the  three  pairs  of  orthogonal  stresses  are 
principal  stresses. 

The  cylindrical  pressure  casings  contain  the  linear  variable  differential 
transformers  (LVDTs)  and  the  pressure  transducers  as  shown  in  Fig.  3.9.  They  also 
serve  as  holders  of  the  flexible  membranes.  The  casings  on  opposite  sides  are  connected 
together  by  stainless  steel  tubings  and  thus  equal  pressures  can  be  applied  on  opposite 
sides  of  the  specimen.  All  three  directional  pressures  can  be  controlled  independently, 
either  manually  through  pressure  regulators  or  automatically  through  servo-controlled 
solenoid  valves. 
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Fig.  3.7 


Dimensions  of  Space  Frame  (After  Sivakugan,  1987) 
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Fig.  3.8  Cross-Sectional  View  of  the  Space  Frame  (After  Sivakugan,  1987) 


3.3. 1  Total  and  Pore  Pressure  Measurement 


Pressures  are  applied  through  the  pressure  casings  to  the  flexible  membranes, 
using  "in  house"  compressed  air  up  to  100  psi  and  "bottled"  compressed  nitrogoi  to 
achieve  higher  pressures.  The  pressures  in  the  casings  (and  on  the  membranes)  are 
measured  by  pressure  transducers  attached  on  the  base  plate  of  each  casing  of  the  three 
pairs.  The  pressure  transducers  are  made  by  Data  Instruments,  Inc.  of  Lexington,  MA. 
They  are  model  AB  type  with  a  pressure  range  of  200  psi  and  a  maximum  signal  ou^ut 
of  100  mV. 

Pore  pressure  at  the  center  of  the  specimen  is  measured  by  custom  made  "needle" 
piezometer.  A  thin  needle  (0.820  mm  OD)  is  inserted  through  a  bottom  diagonal  port 
into  the  specimen.  At  the  tip  of  the  needle  several  small  holes  were  drilled  for  the  first 
centimeter,  and  the  tip  was  covered  with  a  No.  200  wire  mesh  to  prevent  clogging  of  the 
tube.  This  technique  worked  well  except  in  one  or  two  tests  out  of  about  a  hundred  tests 
the  needle  got  clogged.  The  piezometer  needle  is  connected  to  a  pressure  transducer  of 
model  AB.  The  pore  pressure  response  is  immediate  and  very  accurate.  The  needle  is 
cleaned  after  every  test  using  compressed  air  to  prevent  drying  of  kaolin  particles  in  the 
needle.  Before  starting  the  test,  the  needle  is  flushed  thoroughly  with  deionized  deaiied 
water,  saturated  and  then  connected  to  the  pressure  transducer. 

All  four  transducers  are  connected  to  a  multiplexer  MUXl  with  a  gain  of  50.0, 
and  then  to  the  analog-to-digital  (A/D)  convertor.  Data  received  from  the  A/D  convertor 
are  processed  by  the  microprocessor.  According  to  the  instructions  given  to  the 
computer  and  microprocessor,  measurements  from  the  LVDTs  and  the  pressure 
transducers  are  analyzed  and  decisions  made.  Based  on  these  decisions,  directions  are 
passed  to  the  relay  board  to  activate  or  inactivate  the  aj^rt^riate  relays.  This  enables 
opening  or  closing  the  solenoid  valves  and  ther^y  automatically  regulating  the  pressures 
in  all  three  directions  independently.  The  schematic  diagram  of  the  servo  control 
interfacing  is  shown  in  Fig.  3.10. 
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3.3.2  Deformation  Measurement 

The  deformations  of  the  specimen  are  measured  using  the  linear  variable 
differential  transformers.  The  LVDTs  are  fixed  on  the  front  plate  of  the  pressure  casing, 
contained  within  the  casing,  and  the  leads  are  taken  through  the  end  cap  to  the  A/D 
convertor,  the  DASH-8  board.  The  LVDTs  are  DC  (q)erated  and  have  a  nominal  linear 
range  of  ±0.250  inches.  The  LVDTs  used  are  manu^tuied  by  Schaevitz  Engineering 
(GCD/GPD-121-2S0  type),  and  are  hermetically  sealed  and  spring  loaded.  The 
sensitivity  of  the  LVDTs  is  1.6  volt/mm.  A  ±15V  DC  external  power  supply  provides 
power  to  all  LVDTs.  The  outputs  of  the  LVDTs  are  taken  to  the  multiplexer  MUXO. 
As  already  described  in  section  3.3.1  the  readings  from  LVDTs  are  analyzed  in  the 
microprocessor  and  the  proper  instructions  are  sent  to  the  solenoid  valves  to  regulate  the 
pressure.  By  regulating  these  pressures,  either  stress  controlled  tests  or  strain  controlled 
tests  can  be  performed. 

A  total  of  eight  LVDTs  are  used  in  the  deformation  measurement,  one  each  in  all 
six  sides  of  the  cube  and  two  more  on  a  lateral  side.  This  arrangemrat  is  to  check  and 
confirm  the  planar  deformation  of  each  face.  The  three  LVDTs  located  staggeredly  on 
one  side  yielded  essentially  the  same  readings  for  all  tests,  thus  conforming  to  planar 
displacement.  Teflon  disks  of  12.7  mm  diameter  with  a  slight  convexity  are  attached  to 
the  tips  of  the  LVDTs  in  order  to  reduce  the  contact  pressure  on  the  membrane  and  the 
sample. 

Since  the  LVDTs  have  a  limited  linear  range  (maximum  of  0.5  in),  the  initial 
position  of  the  LVDTs  is  adjusted  depending  on  the  strain  path  in  order  to  utilize  the 
LVDTs  maximum  linear  range.  For  example,  a  LVDT  would  be  placed  in  its  maximum 
compressed  level  where  a  compressive  deformation  would  be  measured  on  that  side.  On 
the  other  hand,  initial  position  of  an  LVDT  will  be  set  with  fully  extended  position  of 
rod  where  it  will  measure  a  maximum  extensive  deformation  of  the  sample. 


3.4  Properties  of  Soils  Used 


Two  different  soils  were  utilized  in  the  testing  program.  One  is  Georgia  kaolinite 
clay  in  dry  powdered  form,  obtained  from  Akrochem  Corporation,  Akron,  Ohio  under 
the  trade  name  Akrochem  SC-25  (soft  clay  #25).  The  other  soil  used  is  ground  silica 
obtained  from  U.S.  Silica  Company,  Ottawa,  Illinois  under  the  commercial  name  Sil-Co- 
Sil  #270. 

The  Atterberg  limits  and  specifrc  gravity  of  the  soils  are  given  in  Table  3.2,  and 
their  particle  size  distributions  are  shown  in  Fig.  3.11. 


Table  3.2  Properties  of  the  Kaolinite  and  Kaolin-Silica  Mix 


Soil 

Liquid  Limit 

Plastic  Limit 

Specific  Gravity 

Kaolinite 

63 

33 

2.60 

Kaolin-Silica 

37 

22 

2.65 

Additional  characteristics  of  the  kaolinite  and  ground  silica  are  given  in  Appendix 
D. 

3.5  Slurry  Preparation 

A  proper  slurry  preparation  method  is  an  important  factor  to  obtain  a  high  quality 
uniform  specimen.  This  includes  an  appropriate  water  content,  proper  mixer  and  a 
vacuuming  system.  According  to  Sheeran  and  Krizek  (1971),  samples  prepared  from  the 
higher  water  content  slurry  will  not  be  influenced  by  the  techniques  adopted  to  place  the 
slurry  into  the  slurry  consolidometer.  However,  higher  water  content  slurries  requires 
a  larger  equipment,  longer  piston  and  larger  loading  frame.  From  their  experiments, 
Sheeran  and  Krizek  found  out  that  a  water  content  of  2  to  2.5  times  the  liquid  limit  of 
the  soil  would  be  an  ideal  amount. 
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Fig.  3. 1 1  Particle  Size  Distribution  of  Kaoliniie  and  Ground  Silica 
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To  obtain  kaolin  specimens,  1700  g  dry  powdered  kaolin  was  added  to  2380  g 
(140%  water  content)  deionized  and  deaired  water,  and  mixed  thoroughly  by  a  hand  held 
mixer  for  about  10  minutes.  [Note:  If  the  clay  is  put  first  in  the  vessel  and  then  water 
is  added,  clay  sticks  to  the  bottom  and  wall  and  it  is  very  difficult  to  mix  because  of  clay 
lumps.]  Then  the  slurry  was  subjected  to  a  vacuum  of  about  6S0  mm  mercury  by  an 
aspirator  for  about  4  hours.  After  the  vacuuming,  soil  particles  were  separated  from 
water,  and  therefore  the  mix  was  stirred  by  a  spatula  slowly  and  gratly  but  thoroughly 
without  inducing  any  air  into  the  slurry,  until  it  became  a  smooth  homogeneous  mix. 
Before  the  slurry  was  poured  into  the  consolidometer  (as  shown  in  Fig.  3.1),  silicone 
oil  was  applied  to  the  inner  wall  in  order  to  reduce  the  wall  friction.  Wet  filter  papers 
were  placed  in  between  the  slurry  and  porous  stones  (on  the  base  and  under  the  piston). 
Gaskets  were  provided  in  between  the  upper  and  lower  chamber,  and  lower  chamber  and 
base  to  avoid  leaks.  A  rubber  seal  was  attached  around  the  sides  of  the  piston  to  capture 
any  squeezing  of  slurry  between  the  piston  and  the  walls. 

3.6  Consolidation  in  the  Slurry  Consolidometer 

Initially  a  small  pressure  of  3  psi  is  applied  to  the  piston  for  a  day.  This  prevents 
any  leak  between  the  compartments  or  between  the  piston  and  walls  and  the  slurry 
becomes  semi-solid.  After  the  semi-solid  cake  is  formed  closer  to  the  porous  stones, 
there  will  not  be  any  leak  even  at  high  pressures.  The  second  pressure  increment  of  18 
psi  is  applied  for  another  day  to  complete  the  slurry  consolidation. 

From  the  Cy  values  of  the  soils  tested  it  was  estimated  the  primary  consolidation 
would  be  over  within  a  day.  Thus,  for  both  soils,  two  days  were  allowed  to  complete 
the  consolidation.  The  16  holes  of  14"  diameter  drilled  through  the  piston  for  top 
drainage  also  help  to  remove  the  piston  without  suction.  The  lower  chamber  is 
reinforced  externally  to  obtain  exactly  a  four  inch  cube  sample. 


After  the  slurry  consolidation,  the  consolidometer  is  removed  from  the  loading 
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frame,  and  the  water  collected  on  the  top  is  poured  into  a  sink.  After  removing  the 
screws  fastening  the  upper  and  lower  chambers,  the  chambers  are  separated  slowly  using 
a  flat  screwdriver,  allowing  raough  time  to  equalize  the  internal  pressure  to  atmo^heric. 
The  excess  consolidated  cake  is  cut  with  a  wire  cutter  and  the  top  surface  of  the  sample 
is  trimmed  to  the  exact  size  with  a  sharp  straight  edge.  A  wet  4x4  filter  paper  was 
placed  on  the  top  surface.  The  lower  chamber  with  the  specimen  was  removed  from  the 
base  and  the  bottom  surface  of  the  sample  was  covered  with  a  filter  paper  and  transported 
to  an  extrusion  jack.  A  teflon  block  was  used  to  extrude  the  sample  and  after  extruding 
all  other  four  sides  were  covered  with  wet  filter  papers.  Then  the  sample  was  taken  to 
the  space  frame  with  the  teflon  block. 

The  wall  friction  plays  an  important  role  on  the  d^sition,  sedimentation  and 
consolidation  of  the  slurry.  Deschamps  (1991)  found  that  the  water  content  varied  about 
15%  from  top  to  bottom  of  a  sample  consolidated  from  a  montmorillonite  slurry,  and 
from  empirical  correlations  estimated  that  stress  differences  of  25%  existed  between  the 
top  and  bottom.  To  minimize  frictional  effect  in  the  lower  compartment  of  the  slurry 
consolidometer  where  the  test  specimen  is  collected,  the  walls  are  lined  with  teflon 
sheets.  In  addition,  silicone  oil  is  applied  to  the  teflon  surface  before  the  slurry  is  poured 
into.  To  check  the  efficiency  of  this  method  in  reducing  wall  friction,  a  four  inch  cube 
sample  from  a  consolidated  slurry  was  cut  into  64  one  inch  cubes,  numbered  as  shown 
in  Fig.  3.12,  and  the  water  contents  of  the  64  samples  were  determined  (Fig.  3. 13).  The 
water  content  distribution  is  almost  uniform  throughout  the  sample.  The  overall  variation 
is  only  of  ±0.8%  of  the  mean,  indicative  that  the  effect  of  wall  friction  is  almost 
negligible. 

3.7  Placing,  Flushing  and  Back  Pressure  Saturation 

Before  the  sample  which  is  covered  with  wet  filter  paper  is  transferred  to  the 
space  frame,  the  following  steps  are  followed.  The  tip  of  the  piezometer  needle  is 
covered  with  200  wire  mesh  cloth.  Then  the  needle  is  thoroughly  flushed  and  filled  with 
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Fig.  3.12  Cubical  Specimen  for  Water  Content  Determination 
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For  Kaolin  Soil  After  Slurry  Consolidation  of  17  psi 


Fig.  3.13 


Water  Content  Distribution  within  Cubical  Specimen 
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water,  and  the  pressure  transducer  is  attached  carefully  without  any  air  bubble  entrapped. 
To  minimize  friction  during  the  testing,  silicone  oil  is  applied  to  the  space  frame  and  to 
the  membranes.  A  very  small  pressure  is  applied  to  the  bottom  membrane,  just  enough 
to  carry  the  self  weight  of  the  sample  and  the  teflon  block.  The  sample  is  brought  to  the 
space  frame  with  the  teflon  block  and  placed  through  the  top  opening.  The  teflon  block 
is  carefully  slipped  out  and  the  sample  adjusted  exactly  at  the  center.  The  bottom 
spaghetti  tube  is  left  beneath  the  sample  at  the  center,  and  the  top  spaghetti  tube  is  placed 
on  the  top  center  of  the  sample.  These  spaghetti  tubes  are  used  for  flushing,  back 
saturation,  back  pressurizing  and  for  draining  during  the  consolidation  phase.  The 
piezometer  needle  (attached  with  the  pressure  transducer)  is  inserted  through  a  diagonal 
port  in  the  bottom  and  guided  to  penetrate  diagonally  through  the  sample  such  that  the 
tip  reaches  the  center  of  the  sample.  All  five  sides  of  the  specimen  are  covered  by  the 
membranes  and  the  pressure  casings  are  fixed  tightly  to  hold  the  membranes  in  place. 
To  measure  the  deformations,  LVDTs  had  already  been  placed  inside  the  casings  and 
their  positions  adjusted  according  to  the  type  of  test  to  be  run,  in  order  to  use  their 
maximum  range. 

The  "  stainless  steel  tubings  are  connected  to  the  opposite  pressure  casings,  and 
to  the  outlets  of  normally  closed  and  normally  opened  solenoid  valves.  After  all  the 
connections  were  tightened,  the  computer  program  "stress.bas"  (given  in  Appendix  E) 
is  loaded  and  executed,  to  measure  the  pressures  and  deformations  and  then  to  record  the 
data. 


In  order  to  expel  the  air  entrapped  in  between  the  sample  and  the  membranes  a 
small  pressure  of  5  psi  is  applied  in  all  directions.  The  outlet  tubing  for  flushing  is  kept 
in  a  beaker  with  water  so  that  the  expulsion  of  air  can  be  monitored.  After  the  air 
bubbles  have  stopped,  the  valve  for  back  pressure  saturation  is  opened  and  a  small  water 
pressure  is  applied  through  the  burette.  By  closing  and  opening  the  flushing  valve  most 
of  the  air  pockets  are  removed.  Remaining  air  is  dissolved  in  the  water  during  the  back 
pressure  saturation  stage.  The  all  around  pressures  and  the  back  pressure  are  increased 


slowly  and  simultaneously  by  5  psi.  The  sample  is  left  for  one  day  under  that  pressure 
(10  psi  cell  and  5  psi  back  pressure).  The  cell  pressure  and  the  back  pressure  are  then 
increased  to  30  psi  and  25  psi,  respectively,  and  left  for  another  day.  After  checking  the 
B-parameter,  all  pressures  are  increased  by  another  20  psi  and  kept  for  another  day.  The 
difference  between  the  all  around  cell  pressure  and  back  pressure  are  always  kept  at  5 
psi.  The  pressure  increments  are  applied  steadily  and  slowly  at  about  2  psi  per  minute, 
allowing  enough  time  for  pore  pressure  to  stabilize  and  not  allowing  to  prestress  the 
sample  during  the  saturation.  The  applied  back  pressure  is  maintained  throughout  the 
whole  testing  procedure.  At  the  end  of  saturation,  B-values  were  found  greater  than  0.98 
in  all  tests. 

3.8  1-D  Consolidation  in  the  Cuboidal  Shear  Device 

Although  the  data  in  Fig.  3.13  indicate  that  uniform  specimens  are  obtained 
(moisture  content  variation  of  less  than  0.8%  from  the  mean)  by  minimizing  friction  in 
the  lower  chamber  where  the  specimen  is  consolidated,  the  disturbance  during  the 
leveling,  extruding,  handling  during  transferring  and  placing  could  have  some  effect  on 
the  specimen.  To  overcome  this,  after  ensuring  saturation  the  specimen  was  consolidated 
one  dimensionally  using  the  flexible  membrane  boundaries  and  the  servo  controlled 
system  to  higher  stresses  (55  psi  vertical  and  about  35  psi  horizontal)  than  subjected 
during  slurry  consolidation.  It  has  been  shown  (Sivakugan,  1987)  that  after  such 
consolidation  in  the  CSD,  uniformity  of  water  content  throughout  the  specimen  had 
improved  to  ±0.5%  of  the  mean. 

One  dimensional  consolidation  in  the  cuboidal  shear  device  is  achieved  by 
increasing  simultaneously  to  100  psi  all  the  pressures  in  all  three  directions,  keeping  the 
same  back  pressure  of  45  psi  at  the  end  of  saturation  phase  with  the  drainage  valve 
closed.  Since  the  B-parameter  was  close  to  1.0,  the  pore  pressure  has  to  increase  by  the 
same  amount  (i.e.  by  50  psi,  to  95  psi).  The  pore  pressure  dropped  by  2  to  3  psi  and 
remained  constant  after  10  minutes.  After  the  pore  pressure  is  stabilized,  the  drainage 


valve  is  opened  and  the  sample  allowed  to  consolidate  one-dimensionally  by  decreasing 
the  lateral  stresses  through  the  servo-controlled  system  maintaining  very  little 
deformations  in  lateral  directions.  The  flow  chart  for  the  servo  controlled  one 
dimensional  (KJ  consolidation  is  given  in  Fig.  3.14.  A  tolerance  of  0.0005  inch  is 
permitted  for  the  zero  lateral  deformations.  Eighty  percrat  of  the  primary  consolidation 
was  basically  over  within  80  minutes  in  all  tests.  This  is  i^>parrat  from  the  typical  pore 
pressure  vs  time  curves  and  settlement  vs  time  curves  for  kaolin  and  kaolin-silica  mix  in 
Figs.  3.15  to  3.18.  The  primary  consolidation  was  completely  over  within  24  hours. 
The  main  reason  for  the  fast  consolidation  is  the  shorter  drainage  path  provided  by  the 
filter  papers  surrounded  all  six  sides. 

3.9  Stress  Controlled  and  Strain  Controlled  Tests 

One  of  the  main  advantages  of  the  cuboidal  shear  test  is  the  ability  of  achieving 
any  kind  of  stress/strain  path.  All  three  stresses  can  be  applied  independently,  and 
combined  with  the  servo-controlled  system  it  is  possible  to  perform  both  stress  controlled 
and  strain  controlled  tests.  The  one  dimensional  (KJ  consolidation  test  CSD  is  a  typical 
example  of  a  strain  controlled  test,  with  the  lateral  strains  maintained  at  zero. 

The  flexible  membrane  boundary  condition  is  ideal  for  applying  uniform  stress 
on  to  the  specimen.  In  the  CSD,  air  pressures  can  be  directly  applied  to  the  membranes 
and  controlled  through  the  regulators  and  thus  stress  controlled  tests  can  be  easily 
performed.  On  the  other  hand  strain  controlled  tests  can  be  executed  directly  by  rigid 
platen  at  the  expense  of  non  uniform  stress  distribution.  In  the  CSD,  strain  controlled 
condition  is  achieved  indirectly  with  the  aid  of  the  servo-controlled  system.  A  small 
pressure  increment  is  applied  and  the  deformation  is  measured  and  checked  against  the 
condition  given.  If  the  deformation  is  found  higher  than  the  required  value,  the  pressure 
is  decreased  using  normally  closed  solenoid  valve,  and  if  it  is  less  than  the  required,  an 
additional  pressure  increment  can  be  applied  through  the  normally  opened  solenoid  valve. 
The  response  is  immediately  checked  again  and  readjustments  made.  This  iteration 
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process  continues  until  the  required  condition  is  achieved,  and  can  be  independently  and 
simultaneously  performed  for  all  three  directions.  A  flow  chart  for  strain  controlled 
loading  is  shown  in  Fig.  3.19. 

For  example,  to  perform  tests  under  various  strain  rates,  the  deformation 
measured  is  divided  by  the  time  interval  and  checked  with  the  prescribed  strain  rate  and 
the  iteration  continues  until  both  strain  rates  are  equal.  To  achieve  plane  strain  condition 
and  Ko  consolidation  the  strain  is  maintained  to  zero  using  the  same  above  procedure. 

The  speed  of  the  data  acquisition  was  increased  by  several  fold  by  using  Turbo 
Basic  compiler  and  the  compiled  version  of  the  coupled  program. 

3.10  Summary 

The  servo-controlled  cuboidal  device  is  an  inexpensive  but  well  versatile 
equipment  in  which  almost  any  stress  path  can  be  achieved.  With  the  flexible  boundary 
loading  and  minimized  friction  between  the  membrane  and  fhune  using  silicone  oil,  a 
uniform  principal  stress  can  be  applied.  The  boundary  conditions  are  well  defined  and 
the  pressure  and  deformation  measurements  are  obtained  accurately  and  quickly.  Using 
the  servo-control  system,  accurate  K,  consolidation,  strain  controlled  loading  and  any 
type  of  stress  path  tests  can  be  performed  in  the  cuboidal  shear  device. 

Exactly  same  procedures  were  adopted  to  prepare  the  test  specimen,  to  make  sure 
that  the  different  results  obtained  were  only  due  to  the  variation  in  the  test  pix)ceduies. 
The  uniformity  of  the  specimen  was  voified  by  the  moisture  content  profile  (indirect 
measurement  of  void  ratio)  after  the  slurry  consolidation  and  before  performing 
pressuremeter  test.  The  variation  of  the  moisture  content  improved  from  ±0.8%  of  the 
mean  after  the  slurry  consolidation  to  ±0.5%  of  the  mean  after  reconsolidation  in  the 
CSD.  This  moisture  content  variation  is  so  small  and  it  can  be  assumed  that  the 
specimen  preparation  method  is  satisfactory  and  it  yields  uniform  samples. 
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Fig.  3.19  Flow  Chart  for  Strain  Controlled  Test  (After  Sivakugan  et  al.,  1988) 
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CHAPTER  4 

STRAIN  RATE  EFFECTS  ON  PRESSUREMETER  TESTING 

4.1  Introduction 

Even  though  the  shear  strength  of  a  soil  can  be  interpreted  from  a  pressuremeter 
test,  one  can  not  expect  the  same  value  from  a  triaxial  test  because  the  tests  follow 
different  stress  paths.  The  pressuremeter  consistently  yields  high  undrained  shear 
strength,  as  high  as  100%  more  than  the  triaxial  test.  This  difference  was  attributed  to 
the  speed  of  the  test,  because  the  strain  rate  in  the  pressuremeter  (about  1  %/min)  is  about 
100  times  higher  than  the  strain  rate  in  conventional  triaxial  tests  (about  0.01  %/min). 
However,  stress  path  effects  also  contribute  to  these  differences. 

The  influence  of  the  strain  rate  in  triaxial  testing  has  been  thoroughly  studird, 
however  only  a  few  experiments  have  been  performed  to  study  strain  rate  effects  in 
pressuremeter  tests,  and  they  were  nonconclusive.  Therefore,  this  study  was  undertaken 
to  quantify  strain  rate  effects  in  pressuremeter  testing. 


4.2  Effect  of  Rate  of  Strain  in  Triaxial  Tests 


Shear  strength  and  deformation  behavior  of  clays  are  time  dependent.  Several 
^  researchers  (e.g. ,  Casagrande  and  Wilson,  195 1 ;  Crawford,  1959;  Perloff  and  Osterberg, 

1963;  Richardson  and  Whitman,  1963)  have  shown  that  the  undrained  shear  strength 
determined  from  laboratory  tests  depends  on  the  speed  of  testing.  From  the  undrained 
•  triaxial  compression  tests  on  the  plastic  clay  from  Drummen  (Norway),  Bjerrum  (1972) 
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concluded  that  the  shear  strength  was  increasing  about  10%  for  every  ten  fold  increase 
of  strain  rate. 

Vaid  and  Campanella  (1977)  performed  conventional  constant  rate  of  strain  shear, 
constant  stress  creep  and  various  other  tests  on  a  local  undisturbed  clay  (called  Haney 
clay).  Their  tests  showed  a  linear  increase  in  undrained  shear  strength  with  the  log  of 
strain  rate  in  the  higher  strain  rate  regions.  However,  in  the  lower  strain  rate  domain, 
the  undrained  strength  reached  a  limit  (called  upper  yield),  and  a  further  reduction  in  rate 
did  not  result  in  additional  loss  of  strength.  On  the  other  hand,  Bjemim’s  test  results  did 
not  show  this  upper  yield  limit  although  his  lowest  strain  rate  is  less  than  that  of  Vaid 
and  Campanella. 

Nakase  and  Kamei  (1986)  investigated  the  influence  of  strain  rate  on  undrained 
shear  characteristics  of  Ko-consolidated  cohesive  soils  by  performing  triaxial  compression 
and  extension  tests.  They  used  Kawasaki  clay  (PI  =  30)  and  two  reconstituted  soils  of 
Toyura  sand  and  Kawasaki  clay  with  PI  of  15  and  10.  Prapaharan  et  al.  (i989) 
combined  all  the  above  results  and  plotted  them  together  using  shear  stress  at  a  strain  rate 
0.01%  as  the  normalized  value  (Fig.  4.1).  The  upper  yield  strength  was  assumed  to 
occur  at  a  strain  rate  of  0.001  %/min  and  after  that  the  shear  strength  increases  by  8-10% 
for  a  tenfold  increase  in  strain  rate.  All  these  results  were  based  on  conventional  triaxial 
tests  performed  on  various  types  of  soils. 

4.3  Effect  of  Rate  of  Strain  in  Pressuremeter  Test 

The  strain  rate  used  in  the  pressuremeter  test  is  commonly  one  to  two  orders  of 
magnitude  larger  than  that  used  in  laboratory  tests.  Furthermore  the  strain  rate  varies 
inversely  with  the  square  of  the  radius  in  the  soil  mass  around  the  probe  during  the  test 
(Wroth,  1975). 


Fig.  4.1  Normalized  Shear  Strength  Versus  Strain  Rate  from  Triaxial  Tests 
(After  Prapaharan  et  al.,  1989) 


There  are  basically  three  types  of  pressuremeters  available:  (1)  pressuremeters 
used  in  pre-drilled  borehole,  e.g.,  Menard  pressuremeter,  Lateral  Load  Tester  G-LT), 
Elastmeter  100,  and  TEXAM  pressuremeter  (Briaud  et  al.,  1986);  (2)  self  boring 
pressuremeters,  e.g.,  PAF  (Jazequel,  1982),  Camkometer  (Wroth  and  Hughes,  1973), 
Offshore  pressuremeters  such  as  Push  in  Pressuremeter  (PIP)  (Fyffe  et  al.,  1986),  and 
PAM  (Brucy  and  LeTirant,  1986);  and  (3)  displacement  type  such  as  full-displacement 
pressuremeter  (Hughes  and  Robertson,  1985)  and  cone  pressuremeter  (Withers  et  al., 
1986  and  1989).  Generally  the  pre-drilled  pressuremeter  tests  are  stress  controlled  tests 
where  the  pressure  is  applied  in  increments,  each  increment  being  held  for  a  specified 
period  of  time.  During  this  holding  time  some  consolidation  and  creep  will  inevitably 
occur  around  the  expanding  probe.  The  assumption  of  undrained  behavior,  therefore, 
is  questionable. 

Winter  (1982)  suggested  a  procedure  for  pre-drilled  pressuremeter  testing  in  both 
granular  and  cohesive  soils.  It  recommends  to  apply  pressure  in  equal  steps  until  the 
expansion  of  the  probe  during  one  load  increment  exceeds  about  one  fourth  of  the 
original  probe  volume  (typically  200  cm*  for  a  800  cm*  probe).  The  load  increments 
should  be  selected  in  such  a  manner  that  about  ten  load  increments  are  required  to  reach 
the  maximum  loading.  Generally,  20-,  50-,  100-,  or  200-kPa  pressure  steps  are  used. 
The  volumetric  readings  are  taken  15,  30,  and  60  seconds  after  the  load  is  applied.  This 
practice  is  followed  in  order  to  take  care  of  creep  effects  and  to  dissipate  the  pore 
pressures  generated  due  to  the  probe  expansion. 

Briaud  et  al.  (1986)  suggest  both  pressure  and  volume  controlled  tests,  either 
increasing  the  probe  pressure  or  its  volume  in  equal  increments.  From  field  tests,  they 
found  that  the  recommended  pressure  and  volume  increment  procedures  could  be  used 
interchangeably.  The  test  should  be  performed  by  doubling  the  initial  size  of  the  probe 
in  10  minutes.  There  are  merits  and  drawbacks  to  both  the  pressure  and  volume 
increment  procedures.  The  disadvantage  of  the  pressure  increment  procedure  is  that  the 
limit  pressure  p^  must  be  estimated  before  running  the  test,  on  the  other  hand  this  is  the 
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advantage  of  the  volume  increment  procedure  where  no  estimate  of  Pi,  is  required.  The 
disadvantage  of  the  volume  increment  procedure  is  that,  if  the  volume  increments  are  not 
small  enough,  the  modulus  part  of  the  curve  may  not  be  defmed  by  enough  points.  This 
is  the  merit  of  the  pressure  increment  procedure  where  the  modulus  is  well  deftned. 
Considering  these  aspects,  Briaud  et  al.  (1986)  recommended  ten  one  minute  increments 
equal  to  Pl/10  for  the  pressure  incremoit  procedure,  and  forty  15  seconds  increments 
equal  to  WJAO  for  the  volume  increment  procedure  where  V„  is  the  deflated  volume  of 
the  probe. 

The  effects  of  consolidation  and  creep  on  stress  controlled  pressuremeter  tests 
have  been  studied  experimentally  and  numerically  by  Pyrah  et  ai.  (1985)  and  Anderson 
et  al.  (1987).  They  performed  undrained  pressuremeter  tests  in  a  modified  triaxial  cell 
so  that  hollow  cylindrical  specimens  150  mm  OD  and  150  mm  high  with  a  25  mm 
diameter  cylindrical  cavity  (Fig.  4.2)  could  be  subjected  to  internal  pressure  increments 
simulating  expansion  of  a  borehole  during  a  pressuremeter  test.  The  numerical  study 
used  the  finite  element  program  CAMFE  (Carter,  1978),  based  on  a  Biot  type 
consolidation  and  the  modified  cam  clay  model  with  solutions  obtained  from  an 
incremental,  time-marching  technique  that  can  deal  with  both  material  and  finite 
deformation  nonlinearity. 

These  studies  showed  that  the  effect  of  both  consolidation  and  creep  is  to  reduce 
the  deduced  modulus  values,  but  that  the  consolidation  around  the  probe  tends  to  produce 
higher  deduced  undrained  shear  strengths,  while  creep  tends  to  have  the  opposite  effect, 
i.e.  give  a  lower  deduced  strength.  Soil  parameters  derived  from  a  stress-controlled  test 
are  thus  dependent  on  the  relative  effect  of  consolidation  and  creep,  and  for  stress- 
controlled  tests  the  effect  of  creep  appears  to  be  more  critical. 

Generally  the  self  boring  type  pressuremeter  tests  are  performed  by  increasing  the 
probe  volume  at  a  constant  rate  of  volume  increase,  thus  they  are  strain  controlled  tests. 
Pyrah  et  al.  (1988)  and  Huang  et  al.  (1991)  have  studied  strain  controlled  pressuremeter 
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Cross-Section  of  Hollow  Cylinder  Triaxial  Apparatus  (After  Anderson 
etal.,  1987) 


69 


tests.  Similar  to  the  stress  controlled  tests,  Pyrah  and  Anderson  simulated  the  strain 
controlled  pressuremeter  test  in  the  modified  triaxial  apparatus  using  hollow  cylindrical 
samples  (Fig,  4.2).  They  carried  out  series  of  tests  on  two  clays  (kaolin  and  pottery 
clay)  with  radial  strain  rates  ranging  from  0.2%  to  4.0%  per  minute.  The  creep  effects 
were  studied  incorporating  Singh  and  Mitchell  (1968)  creep  model  into  the  CAMFE 
program.  The  strain  controlled  test  was  simulated  by  adopting  displacement  defined  type 
of  loading  at  the  inner  boundary.  Constant  stress  was  assumed  at  the  outer  boundary 
simulating  the  in  situ  horizontal  stress  in  a  pressuremeter  test.  The  numerical  simulations 
using  CAMFE  were  performed  for  strain  rates  ranging  from  0.2%  to  4.0%  per  minute, 
and  maximum  expansion  of  the  cavity  was  taken  as  20%  radial  strain. 

By  considering  consolidation  as  the  sole  time-dependent  phenomenon  in  the 
analysis  (consolidation  only),  the  simulated  expansion  curves  for  tests  at  different  strain 
rates  are  shown  in  Fig,  4.3.  For  comparison  purposes,  Pyrah  et  al.  (1988)  simulated  an 
undrained  test  using  a  very  fast  strain  rate  (40%  per  minute)  with  little  time  allowed  for 
consolidation.  It  can  be  noticed  from  the  curves  that  the  slower  strain  rate  tests  show 
higher  limit  pressure  which  could  be  the  result  of  consolidation  taken  place  and 
subsequent  strength  gain  of  the  soil  around  the  cavity.  Pyrah  et  al.  used  Palmer  (1972) 
method  to  generate  stress-strain  curves  from  the  expansion  curves.  The  derived  stress- 
strain  curves  for  different  strain  rates  are  shown  in  Fig.  4.4.  As  expected  the  faster  rate 
tests  gave  higher  strength  at  the  beginning  but  the  slower  rate  tests  gained  strength  with 
strain  (i.e.  with  duration)  due  to  consolidation.  Faster  strain  rate  tests  show  significant 
amount  of  strain  softening. 

With  creep  also  included  in  the  analysis  (consolidation  and  creep),  the  simulated 
expansion  curves  and  the  interpreted  stress-strain  curves  are  shown  in  Figs.  4.5  and  4.6. 
These  curves  clearly  show  that  the  strength  values  increase  with  the  strain  rate.  From 
this  comparison  (consolidation  alone  and  consolidation  with  creep),  creep  phenomenon 
is  the  main  factor  to  yield  smaller  strength  values  at  lower  strain  rates. 
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Average  radial  strain  iV«) 

Fig.  4.3  Simulated  Expansion  Curves  for  Different  Rates  of  Strain 
(Consolidation  Only)  (After  Pyrah  et  al.,  1988) 
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Fig.  4.5  Simulated  Expansion  Curves  for  Different  Rates  of  Strain 
(Consolidation  and  Creep)  (After  Pyrah  et  al.,  1988) 


40%  &  4%  Min.  cur 


73 


Fig.  4.6  Derived  Stress-Strain  Curves  for  Different  Rates  of 
(Consolidation  and  Creep)  (After  Pyrah  et  al.,  1988) 
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Unfortunately,  the  laboratory  tests  (Fig.  4.7)  conducted  by  Pyrah  et  al.  (1988)  did 
not  support  the  results  of  numerical  simulation.  With  the  exception  of  the  fastest  test 
(4%  per  minute),  which  gives  the  highest  strength  at  the  beginning,  the  other  tests  show 
a  very  irregular  pattern.  This  erratic  behavior  can  be  due  to  the  experimental  procedure 
such  as  the  possibility  of  partial  drainage.  Pyrah  et  al.  concluded  that  creep  is  the  more 
significant  time-dependent  phenomenon  in  stress  controlled  test,  while  it  is  of  minor 
importance  in  a  strain  controlled  test  creep  is  where  consolidation  is  predominant. 
However,  the  numerical  studies  and  the  triaxial  test  results  (Section  4.2)  contradict  their 
conclusion. 

The  only  other  experimental  study  on  strain  rate  effect  on  pressuremeter  test  by 
strain  controlled  method  was  carried  out  by  Huang  (1986)  and  Huang  et  al.  (1991). 
Huang  performed  a  series  of  strain-controlled  model  pressuremeter  tests  in  cohesive  soils 
inside  a  calibration  chamber  to  study  the  strain  rate  effects  on  pressuremeter  expansion 
curve  and  on  derived  shear  strength.  The  strain-controlled  pressuremeter  tests  were 
performed  with  radial  strain  rates  of  0.1%  to  4.4%  per  minute.  Two  types  of  soils 
namely,  kaolinite  (KlOO)  and  kaolinite  and  silt  mixture  (K50)  were  used  in  the  calibration 
chamber  test.  A  typical  comparison  of  the  expansion  curves  (i.e.  probe  pressure,  P„ 
versus  radial  strain,  e,)  and  the  deduced  principal  stress  differences  ((r,-a«)  for  normally 
consolidated  K50  soil  are  shown  in  Fig.  4.8  and  4.9,  respectively.  A  summary  of  the 
derived  parameters  for  all  the  strain  rate  tests  is  given  in  Table  4. 1 .  The  results  indicate 
that  for  normally  consolidated  clays  (Fig.  4.8)  the  initial  shear  modulus,  Gj,  increased 
approximately  five  times  as  strain  rates  varied  from  0.75%  to  4.4%  per  minute.  The 
change  in  Gj  for  tests  in  over-  consolidated  clays  was  insignificant  (Table  4. 1).  The  data 
further  indicated  the  limit  pressure  Pj,  is  relatively  insensitive  to  the  strain  rate. 

The  numerical  studies  performed  by  Anderson,  Pyrah  and  their  co-workers  and 
the  influence  of  strain  rate  in  strength  measurement  by  triaxial  tests  (Section  4.2)  suggest 
that  the  shear  strength  (peak)  increases  with  the  strain  rate.  However,  similar  to  Pyrah 
et  al.’s  (1988)  experiments,  Huang’s  (1986)  experiments  did  not  show  an  increase  in 
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undrained  shear  strength  with  strain  rate.  In  fact,  the  shear  strength  decreased  with  the 
strain  rate  in  most  of  the  tests  performed  with  the  exception  of  two  tests  (CP16  and 
CP20,  Table  4. 1).  To  clarify  this  contradiction  and  to  quantify  the  effect  of  strain  rate 
on  pressuremeter  testing  an  extensive  laboratory  study  was  undertaken  as  described  in 
Section  4.5. 

4.4  Simulation  of  Pressuremeter  Test  in  Laboratory 

Since  the  pressuremeter  is  a  large  in  situ  device  (about  40  in  long),  it  is  difficult 
to  perform  full  scale  testing  to  study  the  influence  of  various  factors  that  affect  its  results. 
Another  problem  in  field  tests  is  that  the  influence  of  natural  soil  variation  is  unknown. 
To  overcome  these  problems  the  PMT  is  simulated  and  studied  in  the  laboratory  in 
known  soils  by  several  ways  such  as  expansion  of  hollow  cylindrical  sample,  using  model 
(small  scale)  pressuremeter,  and  using  a  true  triaxial  apparatus  to  simulate  stresses  on  an 
element  around  the  pressuremeter  probe. 

4.4. 1  Simulation  of  PMT  in  Thick  Hollow  Cylinder  Test 

The  cylindrical  cavity  expansion  has  been  studied  theoretically,  numerically  and 
experimentally.  Some  of  the  geotechnical  applications  include  bearing  capacity  of  piles 
(e.g.,  Vesic,  1977,  and  Sayed  and  Hamed,  1987),  pullout  resistance  of  anchors  (Vesic, 
1971)  and  interpretation  of  pressuremeter  tests  (e.g.,  Baguelin  et  al.,  1972,  Ladanyi, 
1972,  Palmer,  1972  and  Vesic,  1972). 

Pressuremeter  tests  were  simulated  in  the  laboratory  by  expanding  cylindrical 
cavities  of  thick  hollow  samples  by  Anderson  and  his  co-worters  at  the  University  of 
Sheffield  (Anderson  and  Pyrah,  1986,  Anderson  et  al.,  1987,  Pyrah  et  al.,  1988  and 
Anderson  and  Pyrah,  1989).  A  triaxial  cell  was  modified,  as  shown  in  Fig.  4.2,  so  that 
hollow  cylindrical  specimens  with  a  150  mm  OD  and  £q)proximately  150  mm  high,  and 
a  25  mm  diameter  cylindrical  cavity  could  be  subjected  to  internal  pressure  increments 
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simulating  expansion  of  a  borehole  during  a  pressuremeter  test.  Stress  controlled  tests 
were  performed  on  three  different  clays,  pottery  clay,  Bamsley-clay  and  kaolin  by 
applying  pressures  in  increments  and  using  different  holding  times  (30  sec,  1  min,  2  min 
and  until  creep  ceased).  Strain  controlled  pressuremeter  tests  were  simulated  to  study 
the  strain  rate  effects  on  two  clay  soils  (kaolin  and  pottery  clay)  with  radial  strain  rates 
ranging  from  0.2%  to  4%  per  minute.  Numerical  simulation  using  finite  element 
technique  was  also  performed  to  simulate  the  above  tests  and  effects  of  consolidation  and 
creep  were  studied  both  separately  and  combined.  The  results  were  presented  and 
discussed  in  the  previous  section. 

4.4.2  Simulation  of  PMT  by  Model  Pressuremeter 

Gangopadhyay  and  Nasr  (1986),  and  Nasr  and  Gangopadhyay  (1988)  used  a 
model  pressuremeter  of  l.S  cm  diameter  in  a  IS  cm  diameter  and  10  cm  high  artificially 
sedimented  kaolin  specimen  (Figs.  4.10  and  4.11).  To  avoid  the  disturbance  during  the 
insertion  of  the  probe,  the  probe  was  placed  first  and  kaolin  slurry  was  consolidated 
around  the  probe.  A  series  of  17  pressuremeter  model  tests  were  performed  using 
different  maximum  vertical  pressures  and  overconsolidation  ratios  of  1,  2,  4  and  8.  The 
undrained  shear  strength  of  the  laboratory  pressuremeter  model  samples  were  predicted 
for  each  test  by  the  inteipretation  methods  proposed  by  Menard  in  1957  (Baguelin  et  al., 
1978),  Gibson  and  Anderson  (1961),  Baguelin  et  al.  (1972),  Ladanyi  (1972),  Palmer 
(1972),  Prevost  and  Hoeg  (1975),  and  Denby  and  Clough  (1980).  A  wide  range  in  shear 
strength  values  obtained  with  the  different  inteipretation  methods  was  observed  for  the 
same  test,  and  the  difference  increased  with  the  OCR.  Prevost’s  method  gave  the  best 
agreement  with  the  results  of  triaxial  test. 

Huang  (1986),  Huang  et  al.  (1988)  and  Huang  et  al.  (1991)  performed  a  series 
of  strain  controlled  model  pressuremeter  tests  in  cohesive  soils  [kaolin  (KlOO)  and  kaolin- 
silt  mix  (K50)]  inside  a  calibration  chamber  (Figs.  4.12  and  4.13).  Their  study  showed 
that  the  initial  shear  modulus  Gj  and  undrained  shear  strength  are  sensitive  to  both  strain 


Sample  Preparation  for  Model  PMT  (A)  &  (B)  Sedimentation  Stage  (C) 
Consolidation  Stage  (After  Gangopadhyay  and  Nasr,  1986) 


Shear  Stage  of  PMT  (Field  Test  Simulation)  (After  Gangopadhyay  and 
Nasr,  1986) 


lower  compartment  upper  compartment 


Fig.  4.12  Cross-Section  of  Chamber  Slurry  Consolidation  (After  Huang  et  al., 
1988) 
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Fig.  4.13  Cross-Section  of  Double  Wall  Calibration  Chamber  (After  Huang  et 
al.,  1988) 
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rates  and  soil  disturbance.  For  the  N.C.  KlOO  and  KSO  clays,  Gj  increases  by  four  to 
five  times  as  strain  rates  vary  from  0.73%  to  4.4%  per  min.  But  for  the  O.C.  clays  the 
strain  rate  effects  on  Gj  are  less  dramatic.  As  Gangopadhyay  and  Nasr  (1986)  obtained, 
Huang  et  al.  (1991)  also  found  out  that  the  undrained  shear  strength  values  intopreted 
from  the  model  pressuremeter  test  by  Prevost’s  (1979)  method  were  very  close  to  the 
triaxial  values. 

4.4.3  Simulation  of  PMT  in  True  Triaxial  Test 

Wood  and  Wroth  (1977)  simulated  pressuremeter  expansion  tests  in  the  laboratory 
by  testing  single  elements  of  soil  in  a  true  triaxial  device.  The  strain  path  of  the  PMT 
on  an  element  involves  one  dimensional  consolidation  (to  simulate  the  field  condition) 
followed  by  shearing  under  conditions  of  plane  strain,  at  constant  volume,  in  the  plane 
perpendicular  to  the  direction  of  consolidation.  True  triaxial  devices  are  the  only  types 
of  apparatus  that  can  apply  this  complete  strain  path  in  one  continuous  operation  without 
the  need  for  unloading,  trimming  and  reorientating  the  sample  (Wroth,  1984).  Wood 
used  the  Cambridge  true  triaxial  apparatus  (Hambly,  1969),  i.e.  the  rigid  platen  type 
where  the  strain  controlled  tests  can  be  performed  easily.  Large  deformations  can  be 
achieved  because  each  of  the  sides  of  the  cuboidal  sample  can  be  indq)endently  varied 
between  70  and  130  mm  (a  maximum  strain  of  30%). 

The  clay  used  for  the  laboratory  tests  by  Wood  and  Wroth  was  Spestone  kaolin 
(LL  =  72,  PL  =  40).  The  kaolin  was  mixed  as  a  slurry  at  a  water  content  of  160%, 
then  pumped  into  the  true  triaxial  apparatus,  and  given  the  required  history  of  1-D 
consolidation  and  unloading.  The  value  of  obtained  from  the  1-D  consolidation  was 
about  0.7.  By  consolidating  kaolin  in  an  oedometer  with  measuremrat  of  horizontal 
stress,  Nadarajah  (1973)  obtained  a  value  of  0.64  for  K„.  Sketchley  (1973)  using  the 
biaxial  plane  strain  apparatus  found  a  value  of  0.63.  Skandarajah  et  al.  (1991)  obtained 
K„  values  of  0.62  after  1-D  consolidation  in  the  cuboidal  shear  device  for  Georgia  kaolin 
(LL  =  63  and  PL  =  33). 


Wood  and  Wroth  (1977)  presented  stress-strain  curves  and  total  and  effective 
stress  paths  for  pressuremeter  tests  simulated  in  the  true  triaxial  apparatus  on  three 
different  clays  at  different  overconsolidation  ratios.  Patterns  of  b^avior  observed  in  the 
ground  and  in  the  laboratory  are  compared  and  the  validity  of  some  of  the  assumptions 
made  in  interpreting  the  pressuremeter  test  results  is  assessed  in  the  light  of  these 
comparisons. 

4.4.4  Full  Size  PMT  in  the  Laboratory  (in  Calibration  Chamber) 

Bellotti  et  al.  (1989)  used  a  SBPM  and  reported  extensive  data  obtained  from  47 
tests  performed  in  a  large  calibration  chamber  using  pluvially-dqposited  silica  sand,  and 
from  25  tests  performed  in  situ  in  a  natural  deposit  of  relatively  clean  silica  sand  of  the 
River  Po,  Italy.  Based  on  the  results  of  these  experiments  a  method  to  correct  the 
measured  unload-reload  shear  modulus  from  SBPM  tests  in  sands  was  proposed  to 
account  for  the  variations  in  stress  and  strain  levels.  The  ENEL-CRIS  calibration 
chamber  (Fig.  4.14)  was  designed  to  calibrate  and  evaluate  different  in  situ  testing 
devices  in  sands.  The  equipment  consists  of  a  double-wall  chamber,  a  loading  frame, 
a  mass  sand  spreader  for  sand  deposition  and  a  saturation  system.  The  chamber  can  test 
a  cylindrical  sample  of  sand  1.2  m  in  diameter  and  1.5  m  in  height.  The  cylindrical 
probe  of  the  SBPM  was  the  Camkometer  Mark  Vm  type  82  mm  in  diameter  and  490 
mm  long  (L/D  is  about  6). 

Anderson  and  Pyrah  (1989)  carried  out  full  scale  80  mm  dia.  SBPM  tests  in  clay 
in  a  triaxial  calibration  chamber  (Fig.  4.15).  A  rigid  cylinder  1.7  m  high  and  0.8  m  in 
diameter  was  filled  with  kaolin  slurry  mixed  at  1.5  times  LL,  which  was  consolidated 
one  dimensionally  using  top  and  bottom  drainage.  When  the  height  of  the  slurry  had 
reduced  to  about  1.0  m,  the  rigid  cylinder  was  removed  and  an  outer  membrane  was 
fitted.  The  SBPM  was  bored  into  the  soil  and  the  clay  consolidated  again  around  it  using 
equal  vertical  and  horizontal  stresses.  This  attempted  to  simulate  perfect  insertion  of  the 
SBPM  with  minimal  disturbance.  They  found  a  time  lag  between  the  pore  pressure 
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transducers  placed  in  the  soil  and  the  SBPM  pore  pressure  cells,  casting  doubt  on  the 
validity  of  the  pore  pressures  measured  by  current  commercially  available 
pressuremeters.  Furthermore,  the  pressuremeter  tests  carried  out  in  the  calibration 
chamber  confirmed  the  influence  of  time  dependent  phoiomena  such  as  consolidation  and 
creep,  and  also  suggested  the  possibility  of  vertical  consolidation  away  from  the  mid¬ 
height  of  the  pressuremeter. 

4.4.5  Full  Scale  Self-Boring  Pressuremeter  Testings  in  Soft  Clay 

Benoit  and  Clough  (1986)  conducted  32  SBPM  tests  in  soft  San  Francisco  Bay 
Mud  wherein  key  parameters  were  varied  one  at  a  time  to  isolate  their  influence.  The 
parameters  considered  included:  expansion  rate,  cutter  position,  cutting  rate,  and  cutting 
shoe  dimension.  Furthermore,  they  determined  lateral  earth  pressure,  shear  strength, 
stress-strain  behavior,  and  in  some  cases,  coefficient  of  consolidation,  from  the  probe 
expansion  data. 

The  strain  rate  or  rate  of  membrane  inflation  is  one  of  the  major  factors 
influencing  undrained  shear  strength  from  SBPM  tests.  It  could  be  due  to  the  possible 
drainage  during  the  test  and  to  the  effects  of  rheologic  parameters.  In  the  pressuremeter 
test,  these  factors  are  bound  together  and  are  difficult  to  isolate.  Benoit  and  Clough 
(1986)  carried  out  seven  SBPM  tests  on  Bay  Mud  to  study  the  influence  of  strain  rate  in 
the  shear  strength.  Young  Bay  Mud  exhibits  the  typical  rheological  behavior  for  clays 
in  undrained  conditions,  wherein  the  strength  increases  with  rate  of  loading  (Lacerda, 
1976).  On  the  other  hand,  Jain  (1985)  determined  from  finite  element  analysis  and 
laboratory  tests  on  Bay  Mud  that  if  drainage  occurs  in  the  pressuremeter  test  for  a 
contractive  soil  like  Bay  Mud  and  the  results  are  interpreted  assuming  undrained 
conditions,  the  strength  will  be  underestimated. 

Membrane  expansion  rates  used  in  the  SBPM  tests  were  varied  from  1.7  to  47.6 
kN/mVmin  (0.24-6.90  psi/min)  as  opposed  to  the  normal  rate  of  6.9  kN/m^/min  (1 
psi/min).  The  result  of  their  tests  is  shown  in  Fig.  4.16  where  the  shear  strengths 
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Fig.  4.16  Compajison  of  Shear  Strengths  Obtained  at  Various  InflaUon  Rates 
(After  Benoit  and  Clough,  1986) 
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obtained  from  the  tests  are  normalized  by  dividing  each  strength  by  the  corresponding 
value  of  strength  determined  in  the  normal  tests  at  that  depth,  and  the  membrane  inflation 
rate  is  normalized  by  the  normal  rate  (1  psi/min).  From  the  normalized  stmigth  versus 
inflation  rate  ratio  plot  (Fig.  4.16),  it  is  observed  that  the  shear  strength  increases 
steadily  by  about  10-15%  per  seven  fold  strain  rate  increase. 

4.5  Simulation  of  Pressuremeter  Test  in  the  Cuboidal  Shear  Device 

As  described  in  the  previous  section,  the  PMT  can  be  simulated  several  ways  in 
the  laboratory,  as  expansion  of  cylindrical  cavity  of  a  thick  hollow  cylindrical  sample  in 
a  modified  triaxial  cell,  small  scale  or  full  scale  PMT  in  a  calibration  chamber  or  in  a 
modified  triaxial  cell,  or  by  testing  a  single  element  of  a  soil  in  a  true  triaxial  apparatus 
subjecting  to  PM  strain  path.  Out  of  all  these  techniques,  testing  an  element  in  a  true 
triaxial  apparatus  is  the  easiest  and  most  versatile  method. 

As  shown  in  Fig.  4. 17,  a  small  cylindrical  element  around  the  cavity  is  simulated 
as  a  cube  element.  Before  the  cavity  expansion,  the  vertical  stress,  is  the  major 
principal  stress,  a\,  and  the  radial  and  hoop  stresses,  which  are  equal,  a\  =  a'l,  are  the 
minor  principal  stresses,  a\  =  tr’s.  When  the  cavity  is  expanded  by  applying  a  pressure 
on  the  cavity  wall  (by  probe  pressure)  the  radial  pressure  increases  and  hoop  pressure 
decreases  and  the  total  vertical  pressure  remains  constant.  At  one  point  the  radial 
pressure  becomes  more  than  the  vertical  pressure  and  thus  rotation  of  principal  stress 
direction  occurs. 

The  pressuremeter  test  is  simulated  in  the  cuboidal  shear  device  (Chapter  3)  in 
the  following  sequence.  Test  specimens  are  prepared  in  a  slurry  consolidometer  from 
the  slurry  at  about  2.5  times  the  LL  water  content.  The  four  inch  cube  specimen  is 
transferred  to  the  CSD  and  reconsolidated  one  dimensionally  to  a  higher  pressure  to 
overcome  disturbance  effects.  Up  to  this  consolidation  phase  exactly  the  same  procedure 
is  followed  for  all  tests,  ^ince  the  pressuremeter  probe  is  very  long,  plane  strain 
condition  can  be  assumed  in  the  plane  perpendicular  to  the  direction  of  consolidation. 
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Fig.  4. 17  Simulation  of  Curvilinear  Element  to  Cubic  Element 
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This  condition  is  implemented  in  the  CSD  by  maintaining  no  deformation  in  the  vertical 
direction,  i.e.  e,  =  0.  The  pressuremeter  test  is  usually  performed  very  fast,  so 
undrained  condition  can  be  assumed  in  clays.  Therefore,  during  the  pressuremeter 
expansion,  the  sharing  occurs  at  constant  volume.  This  condition  is  enforced  in  the  CSD 
by  closing  the  drainage  valves  and  maintaining  =  -£,-  Since  there  is  no  volume  change 
and  the  axial  strain  is  zero,  equating  the  volumes: 

(1  +  £,)(1  +  £.)  =  ! 

yields  =  -€„  neglecting  higher  order  terms.  A  computer  program  was  written  to 
incorporate  the  above  conditions  (c,  =  0  and  e,  =  -e,)  by  adjusting  remotely  (r„  ffy  and 
a*  through  solenoid  valves. 

4.6  Strain  Rate  Tests  Varying  from  0.01%/min  to  5.00%/min 

As  described  in  Section  3.9  various  strain  rate  tests  were  performed  using  a  strain 
controlled  method.  The  strain  rate  can  be  prescribed  for  every  test  and  a  tolerance  of 
0.05  times  the  strain  rate  was  used.  The  tolerance  for  no  deformation  (e^  =  0)  was  given 
as  0,0005  inch. 

After  the  1-D  consolidation  in  the  CSD,  the  stress  path  of  the  pressuremeter 
testing  was  simulated  in  the  CSD  and  the  samples  were  sheared  at  different  strain  rates. 
The  pressures  were  servo-controlled  to  achieve  the  required  strain  path  and  strain  rates. 
At  least  three  tests  were  performed  at  each  strain  rate  (at  0.01  % ,  0.05 % ,  0. 10% ,  0.50% , 
1.00%  and  5.00%  per  min).  In  general  the  undrained  shearing  was  carried  out  until  the 
applied  pressures  reached  200  psi  (maximum  capacity  of  the  pressure  transducers)  or 
until  the  strains  reached  15%.  Very  good  agreement  was  observed  between  tests 
performed  at  the  same  strain  rate.  Hence,  under  controlled  conditions  the  CSD  testing 
yields  very  repeatable  results.  It  was  noticed  after  the  test  that  the  edges  of  the  specimen 
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remained  in  position,  showing  that  the  floating  type  boundary  conditions  minimize  the 
edge  effects. 

4.7  Results  and  Discussion 

The  frequency  of  data  recording  was  selected  deprading  on  the  speed  of  the  test 
(i.e.,  strain  rate  of  shearing).  For  example,  for  the  &stest  test  (strain  rate  = 
S.00%/min),  readings  were  recorded  every  3  seconds  (every  reading),  and  for  the 
slowest  test  (strain  rate  =  0.01  %  per  min)  readings  were  taken  every  IS  seconds  at  the 
beginning  (one  in  S  readings),  and  at  about  S  minutes  interval  (one  in  100)  after  the 
stress  reached  the  plateau.  This  procedure  was  adopted  in  order  to  avoid  storing  too 
much  data  without  losing  accuracy. 

In  the  flexible  membrane  CSD,  strain  controlled  tests  are  performed  indirectly. 
As  described  in  section  3.9,  strain  controlled  loading  is  a  trial  and  error  iteration  process, 
thus  there  is  a  time  lag  in  the  reaction  of  the  soil  and  the  measurement.  Even  if  the 
speed  of  data  acquisition  system  is  increased,  the  solmoid  valves  take  some  time  to  react 
(delay  in  the  mechanical  opening  and  closing  system)  and  thus  slower  response  compared 
to  the  measurement  speed.  So,  there  is  noise  in  the  readings,  particularly  in  the  shear 
stress  data,  because  it  is  the  difference  of  a,  and  Oy  readings,  which  are  not  always 
completely  synchronized. 

To  obtain  a  smoother  curve,  a  regression  program  was  written  to  fit  the  data 
points  into  a  3rd,  4th  and  Sth  order  polynomials.  The  best  polynomial  was  sdected  for 
each  test.  For  the  same  data  set  (e.g.  Test  49),  third,  fourth  and  fifth  order  polynomials 
are  shown  in  Figs.  4. 18-4.20,  respectively.  Out  of  these  three  polynomials,  a  fifth  order 
polynomial  At  was  selected  because  it  had  the  lowest  variation  and  it  is  consistent  with 
other  similar  tests,  i.e.  strain  softoiing  behavior.  For  another  test  (Test  58),  the  best 
suited  fifth  order  polynomial  regression  curve  is  shown  in  Fig.  4.21.  Using  the  above 
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Shear  Strain  (%) 

Fig.  4.20  Curve  Fitting  by  Regression  Analysis  -  5th  Order  Polynomial  -  Test  49 


Shear  Strain  (%) 

Fitting  by  Regression  Analysis  -  5th  Order  Polynomial 
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regression  method  all  experimental  curves  were  smoothened  and  are  given  in  Appendix 
B. 


Some  typical  results  are  shown  in  Fig.  4.22  and  Fig.  4.23  for  kaolin  clay  and 
kaolin-silica  mix  soil,  respectively.  In  order  to  compare  the  test  results  with  various 
strain  rates  they  are  plotted  together,  but  only  three  curves  are  shown  in  each  figure  to 
avoid  congestion.  The  horizontal  axis  represents  the  compression  strain  in  the  x  direction 
of  the  CSD  which  simulates  the  radial  strain  in  the  pressuremeter  test.  The  vertical 
stress  represents  the  shear  stress  which  is  obtained  from  the  half  of  the  difference 
between  the  radial  and  hoop  stresses  C/i 

The  stress-strain  curves  are  very  consistent.  They  show  clearly  the  increase  of 
tangent  modulus  values  and  shear  strength  values  with  the  increase  of  strain  rate.  The 
kaolin  specimens  (Fig.  4.22)  show  a  slight  strain-softening  while  the  kaolin-silica 
specimens  (Fig.  4.23)  exhibit  strain-hardening  effects. 

The  summary  of  the  peak  shear  strength  versus  strain  rate  test  results  are  given 
in  Tables  4.2  and  4.3  for  kaolin  and  kaolin-silica,  respectively.  The  corresponding 
stress-strain  curves  are  attached  in  Appendix  B.  The  shear  strength  values  at  various 
strain  rates  were  normalized  by  the  average  shear  strength  at  0.01  %  per  minute,  in  order 
to  compare  the  test  results  (for  pressuremeter)  with  the  previously  published  results  (e.g. , 
Prapaharan  et  al.,  1989).  The  normalized  shear  strength  values  are  plotted  with 
logarithmic  of  strain  rates  and  shown  in  Figs.  4.24  and  4.25  for  kaolin  and  kaolin-silica, 
respectively. 

From  the  plots  (Figs.  4.24  and  4.25)  and  the  regression  analysis,  the  slopes  of  the 
curves  were  found  as  14.3%  and  15.3%  for  kaolin  and  kaolin-silica,  respectively.  For 
the  tests  performed  in  kaolin  clay,  the  strength  values  fall  very  closely  to  the  straight  line 
(with  r^  =  0.98,  standard  error  of  estimation  only  0.02  and  standard  error  of  coefficient 
of  0.0048),  whereas  for  kaolin-silica  tests  the  variation  is  higher  (r^  =  0.64,  standard 


Table  4.2  Shear  Strength  Normalized  with  Respect  to 
0.01  %  per  min  Versus  Strain  Rate  for 
Kaolinite  Clay 


Test  No. 

Strain  Rate 
(%  per  min) 

t  - 

^  2 

(psi) 

35 

0.01 

11.3 

47 

0.01 

11.5 

48 

0.01 

11.4 

52 

0.01 

11.1 

49 

51 

12.2 

13.0 

32 

0.10 

13.2 

33 

0.10 

13.2 

34 

0.10 

13.15 

44 

0.10 

13.0 

50 

0.50 

14.2 

36 

1.00 

14.8 

37 

1.00 

14.4 

38 

1.00 

14.5 

39 

1.00 

14.9 

40 

5.00 

15.9 

41 

5.00 

15.65 

43 

5.00 

15.60 

0.01%/Bln 


l998 
.015 
1.007 
0.980 

1.071 

1.148 

1.166 

1.166 

1.161 

1.148 

1.254 

1.307 

1.271 

1.280 

1.316 

1.404 

1.382 

1.377 


Table  4.3  Shear  Stroigth  Normalized  with  Respect 
to  0.01%  per  min  Versus  Strain  Rate 


for  Kaolin-Silica  Mixture 

Test  No. 

Strain  Rate 
(%  per  min) 

-  . 

2 

(psi) 

0.01%/inln 

55 

0.01 

11.7 

0.946 

56 

0.01 

12.1 

0.979 

64 

0.01 

13.15 

1.064 

78 

0.01 

12.5 

1.011 

61 

16.0 

1.294 

71 

14.5 

1.173 

77 

0.05 

12.8 

1.035 

58 

0.10 

14.4 

1.165 

59 

0.10 

17.0 

1.375 

60 

0.10 

17.0 

1.375 

74 

0.10 

16.0 

1.294  j 

66 

0.50 

16.7 

1.351 

67 

0.50 

17.5 

1.416 

72 

0.50 

16.5 

1.335 

75 

0.50 

14.0 

1.133 

53  1.00  15.9  1.286 

54  1.00  15.7  1.270 

62  1.00  16.4  1.327 

79  1.00  15.5  1.254 


57 

68 


5.00 

5.00 


18.5 

18.4 


1.497 

1.488 


Normalized  Shear  Strength  with  Respect  to  0.01  %/min  Versus  Strain 
Rate  Plot  for  Kaolin  Clay 


2.00 


Fig.  4,25  Normalized  Shear  Strength  with  Respect  to  0.01  %/min  Versus  Strain 
Rate  Plot  for  Kaolin-Silica  Mixture 
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error  of  estimation  =  0.10  and  standard  error  of  coefficient  =  0.0262).  The  reason  for 
this  difference  could  be  related  to  the  heterogeneous  property  of  the  mixture  of  two 
entirely  different  soils  and/or  the  difficulty  to  pin  point  exactly  the  peak  strength  for 
kaolin-silica  curves  because  they  were  strain  hardening  type  (Fig.  4.23).  On  the  other 
hand,  kaolin  clay  is  a  very  uniform  homogeneous  soil  with  a  flaky  structure,  which  can 
be  rearranged  along  the  slip  plane  and  be  the  cause  for  yielding  strain  softening  behavior 
(Fig.  4.22). 

The  straight  line  relationship  obtained  for  kaolin  (Fig.  4.24)  and  kaolin-silica  (Fig. 
4.25)  are  superimposed  on  the  plots  (Fig.  4.1)  developed  by  Prapaharan  et  al.  (1989)  and 
shown  in  Fig.  4.26.  Since  the  curves  in  Fig.  4.1  based  on  triaxial  tests  are  also 
normalized  using  the  shear  strength  value  at  0.01  %  per  min  strain  rate,  it  is  possible  to 
superimpose  and  compare  the  pressuremeter  stress  path  tests  with  the  triaxial  tests.  As 
can  be  noticed  from  Fig.  4.26,  the  undrained  shear  strength  predicted  by  the 
pressuremeter  tests  are  more  sensitive  to  the  strain  rate  than  the  triaxial  tests.  On  the 
average,  it  can  be  concluded  that  the  shear  strength  increases  about  15%  for  every  ten 
fold  of  strain  rate  in  pressuremeter  tests  and  about  10%  in  triaxial  tests.  This  result  is 
consistent  with  the  full  scale  SBPM  tests  (Benoit  and  Clough,  1986)  conducted  on  San 
Francisco  Bay  Mud  where  the  normalized  strength  increases  about  15%  (Fig.  4.16). 

It  is  encouraging  to  notice  that  the  test  results  from  the  CSD  (simulating  PMT) 
are  very  consistent,  and  the  undrained  shear  strength  increases  with  strain  rate,  similar 
to  the  triaxial  tests  (section  4.2).  These  results  are  consistent  with  the  numerical  studies 
performed  by  Pyrah  et  al.  (1988),  which  included  consolidation  and  creep  effects  (Fig. 
4.6).  But  the  laboratory  tests  performed  in  thick  hollow  cylinder  by  Pyrah  et  al.  (1988) 
and  in  model  pressuremeter  in  a  calibration  chamber  by  Huang  et  al.  (1991)  did  not  give 
any  conclusive  results  (already  discussed  in  section  4.3). 

Even  though  the  model  pressuremeter  test  in  a  calibration  chamber  and  the  hollow 
cylinder  test  may  be  expected  to  physically  represent  the  actual  pressuremeter  test  better 
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than  the  simulation  of  pressuremeter  strain  path  in  a  cuboidal  shear  device,  the  latter  test 
had  better  control  in  the  simulation  of  undrained  condition  and  better  measurement 
system.  From  the  erratic  nature  and  odd  shapes  of  some  curves  in  Fig.  4.7,  it  is  very 
clear  that  testing  procedures  must  have  beat  a  problem  in  the  thick  hollow  cylindrical 
tests  performed  by  Pyrah  et  al.  (1988).  As  Huang  et  al.  (1991)  indicates  in  their  paper, 
the  partial  drainage  between  the  model  pressuremeter  probe  and  surrounding  soil  was  an 
inevitable  problem  in  the  calibration  chamber  testing.  Therefore,  the  results  from  the 
simulation  in  cuboidal  shear  device  and  the  numerical  studies  appear  more  definite  and 
consistent  than  that  of  previous  studies. 

4.8  Conclusion 

The  shear  strength  and  deformation  behavior  of  clays  are  time  dependent. 
Usually  the  faster  the  load  applied  the  higher  the  strength  of  soil.  This  phenomenon  is 
mainly  due  to  the  creep  in  between  the  soil  particles.  When  a  load  is  applied  slowly,  it 
has  more  time  to  sustain  the  load  hence  it  creq^s  more  and  looses  its  strength.  On  the 
other  hand,  when  the  load  is  applied  quickly,  the  soil  has  less  time  to  creq>  and  it  shows 
higher  strength. 

From  the  compression  and  extension  tests  performed  in  the  conventional  triaxial 
tests  on  various  soils,  it  has  been  found  that  the  undrained  shear  strength  increases  about 
10%  for  every  tenfold  of  strain  rate.  There  is  no  influence  in  shear  strength  below  a 
strain  rate  of  0.001%  per  min.  Since  the  pressuremeter  test  follows  a  different  stress 
path  than  a  triaxial  test,  it  was  anticipated  that  the  strain  rate  would  have  a  different 
influence  in  the  pressuremeter  test. 

Anderson,  Pyrah  and  their  co-workers  studied  the  effects  of  consolidation  and 
creep  in  pressuremeter  testings  by  both  experimentally  using  thick  hollow  cylindrical  soil 
specimens  and  numerically  using  a  finite  element  program  CAMFE  (Carter,  1978)  for 
both  stress  controlled  and  strain  controlled  conditions.  Unfortunately,  the  experimental 
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tests  did  not  yield  any  conclusive  results,  but  the  numerical  studies  shed  some  light. 
These  studies  showed  that  the  effect  of  both  consolidation  and  creq>  is  to  reduce  the 
deduced  modulus  values,  but  that  the  consolidation  around  the  probe  tends  to  produce 
higher  shear  strengths,  while  creep  toids  to  have  the  opposite  effect,  i.e.,  give  a  lower 
deduced  strength.  Soil  parameters  derived  from  a  stress-controlled  test  are  thus 
dependent  on  the  relative  effects  of  consolidation  and  creq>,  and  the  effect  of  creq) 
appears  to  be  more  influential. 

The  same  conclusions  were  reached  for  the  strain  controlled  pressuremeter  tests. 
When  only  the  consolidation  is  considered,  slower  strain  rate  tests  show  higher  limit 
pressure  which  could  be  the  result  of  consolidation  taking  place  and  subsequent  strength 
gain  of  the  soil  around  the  probe.  The  faster  strain  rate  tests  gave  higher  strength  at  the 
beginning  but  the  slower  rate  tests  gained  strength  with  duration  due  to  consolidation. 
Faster  strain  rate  tests  show  significant  amount  of  strain  softening.  With  creep  also 
included  in  the  numerical  analysis  (consolidation  and  creep),  the  strength  values  increased 
with  the  strain  rate.  From  this  comparison  (consolidation  alone  and  consolidation  with 
creep),  creep  phenomenon  is  the  main  factor  to  yield  smaller  strength  values  at  lower 
strain  rates.  These  numerical  studies  show  trends  similar  to  that  of  triaxial  tests,  with 
the  deduced  undrained  shear  strength  increasing  with  strain  rate,  however,  the  increase 
could  not  be  quantified  because  the  range  of  strain  rate  analyzed  was  only  one  magnitude 
wide  (from  0.2%  to  4%  per  min). 

Similar  to  Pyrah  et  al.  ’s  (1988)  experiments,  Huang’s  (1986)  model  pressuremeter 
tests  in  a  calibration  chamber  also  did  not  show  an  increase  in  shear  strength  with  strain 
rate.  His  tests  suffered  by  the  partial  drainage  along  the  probe-soil  interface. 

The  pressuremeter  strain  path  was  simulated  in  a  cuboidal  shear  device, 
conforming  to  plane  strain  condition  in  vertical  direction  and  undrained  condition  was 
implemented  by  closing  the  drainage  valve  and  maintaining  the  expansion  in  the  y- 
direction  equal  to  the  compression  in  the  x-direction  to  satisfy  no  volume  change.  From 


no 


more  than  forty  well  controlled  tests,  it  has  been  found  that  the  undrained  shear  strength 
increases  about  IS  %  for  every  ten  fold  of  strain  rate  increase.  From  the  seven  full  scale 
SBPM  tests  conducted  on  San  Francisco  Bay  Mud,  Boioit  and  Clough  (1986)  found  out 
that  the  undrained  shear  strength  steadily  increased  by  about  10-15%  for  an  increase  of 
sevra  times  faster  than  the  normal  inflation  rate  of  membrane.  The  consistency  of  the 
results  from  the  laboratory  tests  and  fleld  tests  enhance  the  conclusion  that  the  shear 
strength  increases  about  15%  per  log  cycle. 

The  next  chapter  discusses  how  to  incorporate  the  strain  rate  variation  in  the 
radial  direction  in  the  cavity  expansion  theory  and  quantitatively  estimate  the 
overprediction  of  shear  strength  in  the  pressuremeter  test. 
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CHAPTERS 

ERROR  IN  PRESSUREMETER  TEST  INTERPRETATION 
DUE  TO  STRAIN  RATE  EFFECT 


5.1  Introduction 

The  effects  of  strain  rate  on  triaxial  tests  and  pressuremeter  tests  have  beoi 
discussed  in  Chapter  4.  The  undrained  shear  strength  increases  with  strain  rate,  about 
10%  and  15%  per  every  ten  fold  of  strain  rate  for  triaxial  and  pressuremeter  tests, 
respectively.  During  the  probe  expansion  in  the  PMT,  the  strain  rate  varies  radially 
across  the  soil  mass  surrounding  the  probe.  Theoretically,  the  strain  rate  in  the  probe- 
soil  boundary  is  the  same  as  the  strain  rate  of  probe  expansion  and  the  strain  rate  is  zero 
at  infinity.  Prapaharan  et  al.  (1989)  showed  that  at  a  distance  of  about  ten  radius  of  the 
probe  the  strain  rate  is  very  small  and  can  be  neglected. 

Since  the  strain  rate  varies  hyperbolically  (^own  later),  the  total  effect  of  the 
variation  of  strain  rate  has  to  be  obtained  by  integrating  the  effect  in  each  small  element. 
The  purpose  of  this  chs^ter  is  to  evaluate  quantitatively  the  combined  effect  of  higher 
strain  rate  adopted  in  the  PMT  and  the  variation  of  strain  rate  within  the  soil  mass  on  the 
undrained  shear  stroigth  derived  from  pressuremeter  tests. 

5.2  Expansion  of  Cylindrical  Cavity  -  Brief  Review 

Initially  the  cylindrical  cavity  expansion  was  studied  in  metals  because  of  its 
application  to  pressure  vessels  and  gun  barrels  (Hill,  1950).  The  expansion  of  a 
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cylindrical  cavity  in  soil  mass  was  theoretically  studied  by  Gibson  and  Anderson  (1961) 
assuming  the  soil  as  elastic-perfectly  plastic  with  a  Tresca  yield  criterion,  to  derive  the 
pressuremeter  expansion  curve  for  undrained  tests  in  clays. .  Displacements  were 
calculated  assuming  zero  volume  change  during  plastic  flow,  i.e.  the  soil  was  assumed 
to  behave  as  a  rigid-plastic  incompressible  solid  in  a  plastic  zone  surrounding  the  cavity, 
and  as  a  linear  elastic  solid  beyond  that  zone.  The  effect  of  volume  change  in  the  plastic 
region  was  not  considered.  Using  experimoitally  determined  stress-strain  and  volume 
change-strain  relationships,  Ladanyi  (1963)  introduced  the  volume  change  effects  in  the 
cavity  expansion  analysis. 

In  1972,  Baguelin  et  al.,  Ladanyi  and  Palmer  independently  presented  simple 
interpretation  methods  for  pressuremeter  test,  removing  the  rheological  restrictions 
inherent  in  the  elastic-perfectly  plastic  analysis,  implicitly  imposed  by  Gibson  and 
Anderson  (1961).  The  only  restrictive  assumption  was  that  of  deformation  under 
undrained  conditions. 

Vesic  (1972)  presented  general  solutions  for  the  problems  of  expansion  of 
cylindrical  and  spherical  cavities  in  soil  having  both  cohesion  and  friction.  He  used  a 
linear  elastic-plastic  model  and  considered  compressive  volumetric  strains  during  the 
plastic  phase.  However,  the  elastic-plastic  model  is  unable  to  take  into  account 
volumetric  strains  due  to  shear  (Baguelin  et  al.,  1978).  Nor  can  the  model  consider 
decreases  in  shear  strength  with  strain.  In  other  words  dilatancy  and  sensitivity  are 
ignored. 

Prevost  and  Hoeg  (1975)  proposed  differait  equations  for  strain  hardening  and 
strain  softening  soils.  Using  those  equations  in  association  with  equilibrium 
considerations  and  plasticity  theory,  they  derived  expressions  for  pressuremeter  expansion 
curves. 


Prapaharan  (1987)  found  that  even  though  the  modified  Prevost  and  Hoeg  method 


(Ladd  et  al. ,  1980)  was  the  most  promising  ^ce  it  could  closely  fit  a  curve  while  having 
a  theoretical  background,  it  did  not  fit  theoretical  pressuremeter  expansion  curves  with 
sufficient  accuracy.  Therefore,  Pr^paharan  proposed  an  alternative  equation  which  gave 
a  bettmr  fit  for  theoretical  pressuremeter  expansion  curves. 

Whra  more  sc^histicated  constitutive  models  ate  used  to  rq)resait  the  stress-strain 
relationships,  cavity  expansion  problems  cannot  be  solved  analytically.  In  such  instances 
numerical  methods  are  used  to  solve  the  complicated  equations  (e.g.  Carter  et  al.,  1979). 
Most  classical  solutions  for  stress  and  stnun  distribution  around  an  ^panding  cavity 
which  is  so  far  assumed  in  pressuremeter  test  interpretation,  have  been  based  essentially 
on  a  small  strain  assumption.  Soulie  et  al.  (1986)  presented  a  finite  element  method 
solution  for  large  strain  problems,  based  on  an  incremental  formulation  for  the  problem 
of  an  expanding  cylindrical  cavity. 

5.3  Formulation  of  Cavity  Expansion  with  Strain  Rate 

The  following  section  uses  the  formulation  and  assumptions  from  Prapaharan  et 
al.  (1989): 

1.  Cavity  expansion  is  taking  place  under  undtained  condition. 

2.  Cavity  expansion  is  axisymmetric  and  taking  place  under  plane  strain 
conditions. 

3.  Tensile  strains  are  positive. 

Fig.  S.l(a)  shows  the  initial  in  situ  stress  state  of  a  cavity,  with  center  O  and 
initial  radius  a«.  The  point  B  rq)resents  a  generic  material  point  in  the  soil  mass,  located 
at  a  distance  (r-u).  During  the  cavity  expansion,  the  cavity  radius  increases  from  ao  to 
ao  +  Uo,  and  the  generic  point  moves  from  (r-u)  to  r  (point  B’  in  Fig.  5. 1(b)).  Since  the 
deformation  takes  place  under  undrained  conditions,  the  volume  of  soil  displaced  in  the 
annulus  between  A  and  A’  should  be  equal  to  the  volume  of  material  in  the  annulus 


Fig.  S.l  Expansion  of  Cylindrical  Cavity 


115 


between  B  and  B’.  By  equating  these  volumes: 


»  [(a,  +  =  T  [r*  -  (r  -  u)^ 


5.1 


After  simplification: 


u  =  r  -  [r*  -  (2aUo  +  u,*)]'* 


5.2 


The  positive  root  is  neglected  because  u  can  not  be  larger  than  r.  The  circumferential 
strain  C|  is  tensile  everywhere  and  defined  as: 


= 


r-u 


5.3 


Substituting  for  u  from  Eq.  5.2: 


6,  =  -1  +  [1  -  Uo(2ao  +  uj/r^-'' 


5.4 


The  equilibrium  equation  in  the  cylindrical  system  requires: 


da,  a,  -  o* 

— ^  +  -L - •  =  0 

dr  r 


5.5 


The  difference  between  the  radial  and  circumferoitial  effective  stresses  is  given 
as  a  function  of  the  circumferential  strain  6«,  and  the  strain  rate  e, .  Therefore: 


-  o  •  *  Or  ■  *  9  («•»  ^e) 


5.6 


Substituting  Eq.  5.6  in  Eq.  5.5: 


da,  1 


5.7 


As  r  tends  to  infinity,  the  radial  stress  tends  to  the  in  situ  horizontal  stress  a^, 
which  is  independent  of  time.  At  the  cavity  boundary,  <r,  is  equal  to  the  applied  pressure 
which  is  a  measured  quantity  P(Cg,  a  function  of  the  circumferential  strain  cq  at  the 


cavity  wall  and  strain  rate  ^  for  strain  controlled  expansion.  By  integrating  Eq.  5.7 
from  infinity  to  the  cavity  boundary  (i.e.,  at  r  =  a^  +  uj: 


f  Jo,  «  J  -  -  q(6,  k)dr 

i*(€o»  ^o)  ■ 


To  simplify  the  notation,  the  drcumferaitial  strain  e,  will  be  Fq)resented  by  6 
from  here  onwards.  Using  Eq.  5.4  the  integration  variable  r  is  transformed  to  e: 


1  - 


*  “>)  _  1 


(1  +  6)2 


1  - 


1 


(1  +  €)' 


1 


+  «o) 


5.9 


Differentiating  Eq.  5.9: 


^  _  <fe _ 1 

r’  (1  +  €)’  +  «„) 

Dividing  Eq.  5.10  by  Eq.  5.9: 


^  _ _ Je _ 

r  '  €(1  +  eK2  +  e) 

Substituting  Eq.  5.11  in  Eq.  5.8: 


w 

-  O*  •  / 


e(l  ♦  €)(2  ♦  €) 


~d€ 


5.10 


5.11 


5.12 


If  the  function  q(e,  k)  is  known,  the  Eq.  5.12  can  be  numerically  integrated 
to  obtain  the  strain  rate  dependent  pressuremeter  curve.  In  conventional  interpretation 
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methods,  the  strain  rate  effect  is  not  included,  i.e.,  ^(e,  i)  »  ^e).  In  this  case, 
differentiating  Eq.  5.12  gives: 

*  0(2  * 

0 

which  is  the  same  as  the  equation  derived  by  others  (e.g.  Bagudin  et  al.,  1972,  Ladanyi, 
1972,  and  Palmer  1972). 

5.4  Variation  of  Strain  Rate  within  the  Soil  Mass 


In  Section  5.3,  the  generic  point  in  the  initial  state  (B)  was  taken  at  distance  (r-u) 
and  after  deformation  that  point  (B’)  moved  to  r.  This  step  was  taken  in  order  to 
compare  the  results  with  Palmer’s  (1972)  results.  In  order  to  make  the  differentiation 
possible,  to  get  an  expression  for  strain  rate  in  terms  of  distance,  r  should  be  referred 
to  the  initial  state.  See  Fig.  5.2  for  the  new  distances  of  the  generic  point,  before  and 
after  the  deformation.  When  the  cavity  expands  from  the  initial  radius  a„  to  a«  +  u^,  the 
material  generic  point  B  moves  from  r  to  point  B’  at  r+u.  Following  the  same 
procedure  as  in  Section  5.3,  equating  the  volumes  before  and  after  the  deformation: 

T  [(ao  +  -  a*J  =  T  [(r  +  u)^  -  5.14 


which  leads  to 

u  = -r  +  [r^  +  Uo(2a,  +  ujr  5.15 


The  negative  root  is  neglected  because  u  is  positive  during  probe  expansion.  The  strain 
6,  results  from  Eq.  5.15: 


^  =  -1  + 
r 


1  + 


|V2 


5.16 
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Substituting  €„  =  —  in  Eq.  5.16: 


e  =  -1  +  1  +  (^)2  c^(2  +  6^ 


Differentiating  strain  with  respect  to  time: 


1 


Substituting  Eq.  5.17  in  Eq.  5.18: 


{€,(£^+(2+€^€,} 


*  =  ^  7^  -  2*.0  ♦  O 

2  1  +  €  I  r  j 


To  eliminate  r  from  Eq.  5,19,  I — J  is  obtained  from  Eq.  5.17: 


f£  I  =  (1  c)^  -  1  ^  e(2  £) 

r  €^(2  + 


Substituting  Eq.  5.20  in  Eq,  5.19: 


€  =  €  e(2  f  6) 

^  e,(2  +  (1  +  €) 


The  equation  5.21  gives  the  variation  of  strain  rate  within  the  soil  mass  for  a 
given  strain  €„,  and  strain  rate  at  the  cavity  wall.  Figures  5.3  and  5.4  show  typical 
variations  of  strain  and  strain  rate,  respectively,  with  distance  from  the  c^ter  of  the 
cavity  that  were  obtained  using  equations  5.17  and  5.21.  The  strain  and  strain  rate 


aRCUMFERENTIAL  STRAIN  (%) 


STRAIN  RATE  (%  per  min  ) 


2 


Fig.  5.4 


decrease  with  distance  and  are  essentially  negligible  at  a  distance  equal  to  ten  times  the 
radius  of  the  cavity. 

S.S  Stress-Strain  Modeling  Including  Strain  Rate 

The  variation  of  undrained  shear  strength  with  strain  rate  has  been  discussed  in 
Chapter  4.  Based  on  triaxial  tests  conducted  on  various  soils  with  strain  rate  varying 
from  0.0001%  per  min  to  10%  per  min,  Prapaharan  et  al.  (1989)  proposed  a  bilinear 
model  to  describe  the  variation  of  undrained  shear  straigth  with  strain  rate  (Fig.  4. 1). 
From  the  tests  performed  in  CSD  simulating  pressuremeter  strain  path,  the  influence  of 
strain  rate  in  pressuremeter  tests  was  studied  and  the  outcome  was  included  in  the  same 
plot  and  shown  in  Fig.  4.26. 

Prapaharan  et  al.  (1989)  proposed  the  following  relationship  between  undrained 
strength  and  strain  rate: 

where  qu  is  the  shear  strength  at  a  strain  rate  k,  q,.  is  the  strength  at  a  reference  strain 
rate  a  and  is  slope  of  strength  against  logarithm  of  strain  rate  curve. 

From  Fig.  4.26,  it  can  be  obtained  as  =0.10  for  triaxial  tests  and  =0.15 
for  pressuremeter  tests. 

Assume  a  hyperbolic  model  to  express  the  stress-strain  behavior  of  a  true  strain 
hardening  soil: 


q(€)  = 


6 

D  +  e 


9u 


5.23 


where  q„  is  the  ultimate  strength  and  D  is  a  constant. 


Combining  Eqs.  5.22  and  5.23: 


9(€,  €) 


qjl  *  filogn/^«))  ^ 
(D*€) 


Substituting  Eq.  5.24  in  Eq.  5.12: 


^o)  - 


}  plogipC^g))  ^ 

I  (1  *  e)C2  *  e)(D  *  e) 


5.24 


5.25 


The  strain  rate  £  in  Eq.  5.25  is  a  function  of  strain,  e,  and  can  be  calculated 
from  Eq.  5.21.  Eq.  5.25  can  be  integrated  numerically  to  obtain  the  piessuremeter 
expansion  curve  which  includes  the  strain  rate  effect.  Eq.  5. 13  can  be  used  to  develop 
the  stress-strain  curve  from  the  pressuremeter  curve  which  ignores  the  strain  rate 
variation  in  the  soil  mass.  By  comparing  the  derived  stress-strain  curve  q(6o)  -  which 
does  not  include  the  strain  rate,  against  the  "true”  stress-strain  curve  q(€^  -  which 

includes  the  strain  rate,  error  involved  in  neglecting  the  strain  rate  in  soil  mass  can  be 
estimated. 


The  "true"  pressuremeter  expansion  curves  can  be  developed  from  Eq.  5.25. 
The  parameters  needed  are  q^,  a,  jS  and  the  constant  D  which  describes  the  stress- 
strain  curves  for  strain  hardening  soils.  The  values  chosen  to  these  parameters  and 
constants  are: 


q^:  The  derived  stress-strain  curve  is  normalized  with  respect  to  the  reference 
strength  q„.  Therefore,  for  strain  hardening  soils,  in  Eq.  5.25,  the  term  q„  will 
be  cancelled  out  and  the  absolute  value  is  not  needed. 
a:  The  reference  strain  rate  a  is  tak«i  as  0.01  %  per  min  as  it  is  the  standard  strain 

rate  used  in  laboratory  tests. 

:  Two  strain  rates  are  used:  (i)  1.0%  per  min,  the  strain  rate  recommended  by 
Baguelin  et  al.  (1978)  to  be  used  in  the  PMT,  and  (ii)  0.1%  per  min  for 
comparison  purposes. 


jS:  The  slope  of  the  shear  strength  versus  strain  rate  curve.  From  Fig.  4.26,  the 

slope  for  the  PMT  (from  CSD  tests)  is  about  0.15. 

* 

The  constant  D  which  describe  the  stress-strain  curves  for  strain  hard^iing  soils  can  be 
obtained  from  empirical  relationships  or  from  experimental  data. 

(i)  From  empirical  relationship 

The  assumed  hyperbolic  relationship  to  rq)resent  the  strain  hardening  soil  is  (Eq. 

5.23): 


^  ^  5.26 

d€  (D  +  6)* 

As  €  tends  to  zero,  the  slope  tends  to  q,/D  (Eq.  5.26).  That  is,  the  slope  of  the  stress- 
strain  curve  at  zero  strain  is  equal  to  q/D.  This  slope  can  be  related  to  the  shear 
modulus.  The  shear  modulus  is  defined  as: 

G  =  - 
Y 

where  r  =  q/2  and  7  =  e,  -  c,  =  2c  (for  small  strains). 

Therefore: 

g  =  -SL 
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Hence: 


Initial  slope  =  ^ 

6 


9. 

D 


=  4G 


Therefore: 
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The  shear  modulus  G  is  related  to  the  elastic  modulus: 


2(1  +  V) 

For  undrained  conditions  in  clays,  it  can  be  assumed  that  the  Poisson’s  Ratio  p  is  equal 
to  0.5.  Therefore: 


G 


4£ 


Since  E/q„  values  for  most  clays  fell  between  250  and  750  (Bjemim,  1972),  1/D  will 
range  from  300  to  1000.  For  example,  a  value  of  1/500  was  selected  for  by  Prapaharan 
et  al.  (1989)  in  their  analysis. 

(ii)  Determination  of  D  from  the  experimoital  curves 

From  the  simulated  pressuremeter  tests  ir.  the  CSD  (Chapter  4),  it  was  found  that 
the  kaolin-ground  silica  soil  behaved  like  strain  hardening  soils  (Fig.  4.23).  The 
behavior  of  strain  hardening  soils  can  be  simulated  with  an  hyperbolic  relationship  (Eq. 
5.23): 


q(€)  = 


D  *  € 


9. 


Oi  -  o,  = 


D  +  e 


5.27 


Rearranging  the  terms: 


126 


From  the  experimental  data  (Table  4.3),  -  versus  e  curves  are  plotted  (Fig. 

S.S-S.8,  a  few  examples  in  each  strain  rate)  and  usmg  the  linear  regression  the  slope  and 
the  intercept  are  determined  and  the  calculated  D  is  presrated  in  Table  S.  1.  An  average 
value  of  1/D  ==  160  for  kaolin-silica  was  used  in  the  present  numerical  analysis  for  strain 
hardening  soil. 

S.6  Procedure  for  Numerical  Simulation 

The  relationship  derived  for  the  applied  pressure  on  the  cavity  wall  as  a  function 
of  strain  and  strain  rate  for  strain  hardening  soil  is  (Eq.  S.2S): 

^  ?  1  piog(4/tt)  ^  5.29 

(1  +  €)(2  +  €)(D  +  €) 

Since  q«,  is  the  strength  at  a  reference  strain  rate  a,  it  is  a  known  quantity  and  a 
constant.  Therefore,  it  can  be  separated  from  the  integral  and  brought  to  the  left  hand 
side  of  the  equation,  which  is  in  the  normalized  form. 

Method  I:  Numerical  Integration  by  Incremental  Approach 

Numerical  integration  for  a  complex  function  can  be  performed  using  an 
incremental  method  such  as  Simpson’s  Rule.  Any  function  can  be  tq>proximated  by  a 
higher  order  of  polynomial  and  the  increment  width  a^roximated  to  the  variation  of  the 
function.  For  Simpson’s  Rule  the  range  [a,b]  is  divided  into  n  equal  parts  where  n  is  an 
even  integer.  Then,  according  to  Simpson’s  Rule,  the  integral  can  be  approximated  by: 

h 

[  fix)  dx  =  ^  4y,  +  2y^  +  ...  +  2y..j  +  4y,.,  +  y.) 


5.30 


Linear  Conversion  of  Hypeibolic  Curve  for  0.01%/n)in 


0.0035 


Fig.  5.6  Linear  Conversion  of  Hyperbolic  Curve  for  0.0S%/min 


where  Ax  =  (b  -  a)/n. 


When  n  =  2,  the  Simpson’s  Rule  becomes: 

{fi,x)dx-  Wa)  ♦  *  X»)l 

Now  the  whole  integral  can  be  obtained  by  adding  the  divided  areas  as  shown  in  Fig. 
5.9: 

»*7Ax  b 

/  =  j  fix)  dx  *  f  fx)  dx  *  ...  *  I  fx)  dx  5.32 

a  a*Ax  c*{fm-l)Ax 

Step  1  Step  2  Step  m 

In  the  analysis  reported  herein  each  step  is  repeated  201  times  (i.e.,  n  -  201)  and  1000 
steps  were  taken  (i.e.,  m  =  1000).  These  values  of  n  and  m  were  found  adequate  for 
the  function  considered.  The  integral  in  equation  5.29  is  a  function  of  strain  and  strain 
rate: 


/,  =  /  M 
0 


The  strain  rate  is  related  to  strain  by  Eq.  5.21: 

4  .  4  ^  ‘O  *  (5.34) 

'  4.0  ♦  4.)  (1  ♦  4) 

where  c  varies  with  cavity  strain  e,.  A  variation  of  0%  to  10%  cavity  strain  is 
considered  in  the  numerical  analysis. 

In  Step  1,  a  =  0,  a  +  Ax  =  0.1/201  and  J  is  flitted  into  m  =  1(XX)  points 
to  evaluate  the  integral  numerically.  The  strain  rate  t  is  evaluated  at  every  step  0, 


Aco/lOOO,  2A6o/10(X)  ...  Ae„  and/(€,  e)  is  also  evaluated  at  those  steps.  Thra,  using 
Simpson’s  Rule  (Eq.  5.31): 


4±t 

1000 


...  and  so  on 


can  be  calculated  and  the  sum  of  the  above  integrals  will  give  the  integration  in  Eq.  S.29, 
and  (P  (c^,  6^  -  oj/qo.  can  be  calculated  for  all  varies  from  0  to  0.1.  This  gives  a 
pressuremeter  expansion  curve  which  includes  the  variation  of  strain  rate  in  the  soil  mass 
around  the  probe.  Then  the  pressuremeter  curve  is  differentiated  numerically  using 
Gauss  functions  with  Eq.  S.13  in  order  to  obtain  qfe^). 


The  above  method  is  called  the  "incremental  approach".  A  computer  program 
was  written  in  Quick  Basic  and  the  source  code  is  annexed  in  Appendix  V.  The  results 
are  discussed  in  Section  5.7. 


Method  II:  Numerical  Integration  Consistent  with  the  Derivation  of  the  Formulation 

In  Section  5.3,  the  applied  pressure  e^on  the  cavity  wall  was  found  by 
integrating  the  equilibrium  equation  in  cylindrical  coordinate  system  (Eq.  5.5)  and  using 
the  boundary  conditions,  as  r  tends  to  infinity,  the  radial  stress  tends  to  the  horizontal 
stress  (Th  and  at  the  cavity  boundary  r  =  a„  +  Uo,  a,  is  equal  to  P(€^  the  applied 
pressure  is  related  to  the  shear  stress  of  the  soil  (Eq.  5.8).  When  the  domain  of 
integration  was  switched  from  dr  to  de  Q)y  using  Eq.  5. 1 1),  the  boundary  conditions  were 
also  externally  forced  such  that  at  r  =  a,  +  u^,  6  =  €„  and  at  r  =  oo,  e  =  0,  however 
this  was  not  implicitly  included  or  enforced  anywhere  in  the  formulation.  This  will 
induce  errors  in  Method  I  as  shown  later.  In  Method  I,  without  implying  the 
transformation  of  boundary  condition,  the  integral  was  numerically  calculated  by  simply 
finding  Ae„  by  incremental  approach  and  adding  the  area  under  the  function. 
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Method  n  is  performed  as  follows: 

For  cavity  strain  changing  from  0%  to  10%,  i.e.  =  0.0  to  0.1  (as  earlier), 
using  201  steps  Ae^  =  (0.1  -  0)/201.  Then, 

A*. 

/,  *  /  * 

0 

split  this  integral  into  n  steps  (say  n  =  10(X)  as  earlier)  and  evaluate  e  at  each  of 
Aejn  points  and  calculate  /(e,  £)  at  these  points.  Thai,  calculate  I,  =  edi,  using 
Simpson’s  Rule.  This  gives  one  point  in  the  P(Aeo)  vs  Aeo  curve.  Similarly, 

2A*, 

can  be  calculated  using  rqieat  of  this  process  up  to  201  times  and 

0 

Onally, 

201Ac, 

4)1  =  /  can  be  calculated.  Each  numerical  integral  is  performed  using 

0 

n  =  1000  points. 

The  major  difference  between  this  method  and  the  previous  method  (incremental 
approach)  is  that  the  integration  is  carried  out  every  time  starting  from  0  and  ending  up 
with  cavity  strain,  i.e.  the  entire  soil  mass  is  considered  in  every  small  increment  of 
strain.  It  could  be  argued  that  if  the  same  number  of  integration  points  (n  =  10(X))  is 

Ac.  201  Ac. 

used  for  smaller  increment  such  as  j  and  the  larger  increment  such  as  j  ,  it 

0  0 

might  affect  the  accuracy  of  the  integration.  To  verify  this  aspect,  n  was  increased  by 
100  for  every  iteration  up  to  20,100  points  and  it  was  found  that  the  results  were  not 
improved  by  increasing  the  number  of  integration  points  (Fig.  5.10).  Hence,  n  =  KXX) 
(points)  was  found  to  be  a  reasonable  number  of  integration  points  and  was  used  in  the 
subsequent  numerical  integration. 

A  computer  program  was  written  in  Quick  Basic  to  perform  the  numerical 
integration  and  differentiation  (using  Gauss  functions)  in  order  to  calculate  q  (O- 


Verification  of  Number  of  Integration  Points  in  the  Numerical 
Integration  to  Simulate  Pressuremeter  Test 


137 


'!my 


source  code  is  given  in  Appendix  E.  The  results  are  discussed  in  the  next  section. 

S.7  Results  and  Discussion 

The  following  values  are  used  to  check  the  numerical  integration  with  previously 
published  results  (Prapaharan  et  al.,  1989): 

a  =  Reference  strain  rate  =  0.01%/min 

j8  =  Parameter  for  strain  rate  effect  =  0.10 
D  =  1/500  and 

=  0.1%/min  and  1.0%/min 

Fig.  5,11  shows  the  effect  of  strain  rate  on  undrained  shear  strength  derived  from 
pressuremeter  test  for  strain  hardening  soils  (Prapaharan  et  al.,  1989).  The  solid  lines 
represent  the  so-called  true  material  stress-strain  curves  obtained  from  the  assumed  strain 
hardening  relationship  (Eq.  5.23)  for  different  strain  rates  (0.001,  0.01,  0.1  and 
1.0%/min).  The  dashed  lines  show  the  stress-strain  curves  derived  from  the 

pressuremeter  exansion  curves  obtained  for  two  different  expansion  rates  (0. 1  %/min,  and 
1.0%/min).  These  stress-strain  curves  were  derived  using  Eq.  5.13  from  the 

pressuremeter  curves  obtained  with  Eq.  5.26.  The  difference  between  the  derived  curve 
and  the  true  material  curve  for  the  same  strain  rate  is  an  indication  of  strain  rate  effect 
induced  by  the  pressuremeter  test  condition. 

Using  the  computer  program  developed  for  Method  I,  the  pressuremeter  curves 
are  obtained  by  numerical  integration  procedure  described  in  Section  5.9  using  Eq.  5,25 
and  the  stress-strain  curves  are  derived  from  Eq.  5.13  and  shown  in  Fig.  5.12.  As  can 
be  noticed,  the  stress-strain  curve  derived  using  the  Method  I  (Incremental  Approach) 
coincided  exactly  with  the  "true"  material  stress-strain  curve  as  in  Fig.  5.11.  This  means 
that  the  incremental  approach  did  not  reflect  the  pressuremeter  test  condition  and  the 
algorithm  or  to  be  exact,  the  interpretation  of  the  boundary  condition  is  not  correct. 
Even  though  mathematically  there  is  nothing  wrong  in  the  numerical  integration  by  the 
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incremental  approach,  the  transformed  boundary  conditions  are  not  properly  implemented 
in  the  numerical  integration.  This  is  an  important  point  to  be  noticed  whra  equations 
alone  are  used  in  the  PMT  interpretation  without  relevance  to  the  actual  physical  aspects 
of  the  test. 

Using  Method  n  which  is  consistoit  with  the  transformation  of  the  boundary 
conditions,  the  stress-strain  curve  are  derived  from  the  pressuremeter  expansion  curve 
and  shown  in  Fig.  5.13  with  the  "true"  material  stress-strain  curve.  The  curves  shown 
in  the  figure  exactly  match  the  published  results  (Fig.  5.11)  for  the  strain  rate  0.1%  per 
min. 

Influence  of  Upper  Yield 

In  Chapter  4,  the  influence  of  strain  rate  on  undrained  shear  strength  determined 
by  triaxial  tests  was  discussed.  Bjerrum  (1972)  found  that  the  shear  strength  was 
increasing  about  10%  for  every  ten  fold  increase  of  strain  rate  without  any  upper  yield. 
But  others  (Vaid  and  Campanella  (1977),  and  Nakase  and  Kamei  (1986))  observed 
different  levels  of  upper  yield  values. 

To  study  the  influence  of  upper  yield  in  the  derived  stress-strain  curves  three 
cases  were  considered: 

1.  Upper  yield  at  0.(X)1  %/min  (Priq)aharan  et  al.,  1989) 

2.  Upper  yield  at  0.(XX)1  %/min 

3.  No  upper  yield  (Bjerrum,  1972) 

The  derived  stress-strain  curves  for  the  above  three  cases  are  shown  in  Fig.  5.14. 
As  the  upper  yield  is  decreasing  the  strain  softening  behavior  is  little  more  pronounced. 
In  fact,  even  for  the  no  upper  yield  case  (Bjerrum,  1972),  i.e.  the  strength  2q>proaching 
zero  as  the  strain  rate  tends  to  zero  which  is  an  unrealistic  behavior,  the  increase  of 
strain  softening  is  very  little.  Therefore,  for  all  practical  purposes  it  can  be  assumed  that 
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the  levd  of  upper  yield  is  insignificant  in  the  prediction  of  undrained  shear  strength  from 
the  pressuremeter  test.  This  ascertains  our  previous  notion  that  beyond  a  distance  of 
about  ten  radius  from  the  cavity  wall  the  strain  and  strain  rate  effects  could  be  neglected 
(Figs.  5.3  and  S.4)  because  the  very  small  strain  rates  (0.001%  per  min  and  below)  do 
not  affect  the  strength  predictions. 

Influence  of 

From  Fig.  4.26,  the  slope  of  the  normalized  undrained  shear  strength  versus 
logarithm  of  strain  rate  line,  /3,  was  determined  as  0.10  for  triaxial  tests  and  0.15  for 
pressuremeter  tests.  To  study  the  influence  of  j8  several  parametric  studies  were 
performed.  Fig.  5.15  shows  the  normalized  stress-strain  curves  for  two  strain  rates 
0.1%/min  and  1.0%/min,  D  =  1/500,  reference  strain  rate,  a  =  0.01%/min  and  0  = 
0.10  (the  same  as  Fig.  5.11  of  Prapaharan,  et  al.,  1989).  This  plot  is  rq)roduced  here 
in  order  to  compare  the  influence  of  different  parameters  that  affect  the  stress-strain 
curves.  For  the  higher  value  of  /S  =  0.15,  the  actual  ("true"  material)  and  the  derived 
stress-strain  curves  are  shown  in  Figs.  5.16  and  5.17  along  with  j8  =0.10  curves  for 
cavity  expansion  rate  of  0.1%/min  and  1%/min.  The  curves  are  basically  the  same  at 
the  beginning  (because  they  have  the  same  D)  up  to  a  cavity  strain  of  0.5  % .  Beyond  this 
strain  has  a  significant  influence  on  the  actual  and  derived  curves.  For  the  higher 
value  of  (PMT  condition),  the  peak  strength  shows  a  25%  decrease,  and  as  |3  increases 
the  softening  behavior  also  increases  significantly.  In  fact,  for  =  0.10  (triaxial 
condition)  the  actual  and  PMT  curves  are  strain  hardening  and  for  =  0.15,  they  are 
strain  softening  curves.  As  the  cavity  expansion  rate  increases  from  0.1%/min  to 
1.0%/min,  for  the  same  0  =  0.15,  the  softening  behavior  also  increases,  as  opposed  to 
an  increase  in  strain  hardening  for  0=0. 10.  For  a  0. 1  %/min  cavity  expansion  rate,  the 
peak  strength  is  almost  the  same  for  0  =0.10  and  0.15,  but  for  1.0%/min  expansion 
rate  jS  =  0. 15  produces  slightly  larger  strength  than  j8  =  0. 10.  For  the  higher  strain  rate 
(1.0%/min),  the  PMT  curve  for  =  0.15  is  strain  softening  and  jS  =  0.10  is  strain 
hardening  type.  So,  they  reach  almost  same  (residual)  strength  at  higher  strains. 


Cavity  Strain  % 

Fig.  5.  IS  Effect  of  Strain  Rate  on  Undrained  Shear  Strength  Derived  from  PMT 
(D  =  1/500,  Beta  =  0.1) 


Cavity  Strain  % 

Influence  of  Beta  on  Undrained  Shear  Strength  Derived  from  PMT  at 


Cavity  Strain  % 

Fig.  5. 17  Infl  lence  of  Beta  on  Undrained  Shear  Strength  Derived  from  PMT  at 
Strain  Rate  =  1.0%/min  (D  =  1/500) 
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In  general,  for  the  PMT  condition  =  0.15),  the  derived  curves  from 
pressuremeter  expansion  curves  show  higher  strain  softening  compared  to  the  triaxial  test 
conditions  (j3  =  0.10),  whereas  the  actual  curves  show  higher  strain  hardening  and 
higher  strength  as  increases  from  0.10  to  0.15.  That  is,  the  pressuremeter  test 
condition  induces  more  strain  softening  than  the  triaxial  test  condition  which  very  well 
agrees  with  Prevost  (1976)  who  used  a  strain  rate  sensitive  rheological  model,  and 
Jamiolkowsld  and  Lancellota  (1977)  results. 

Influence  of  D 

The  constant  D  in  the  hyperbolic  equation  which  is  an  inverse  measure  of  initial 
tangent  modulus,  also  has  a  significant  influence  on  the  actual  and  derived  (PMT)  curves 
(Figs.  5.18-5.23).  As  expected,  for  the  smaller  D  value  (which  Prapaharan  et  al.  (1989) 
arbitrarily  selected,  D  =  1/500),  the  actual  and  PMT  curves  have  steeper  slopes  than  for 
the  larger  D  value  (obtained  from  the  CSD  tests  performed  on  kaolin-silica  soil  in 
laboratory;  D  =  1/160)  for  both  cavity  expansion  rates.  The  constant  D  affects  not  only 
the  results  initially  but  also  throughout  the  curve.  For  the  larger  values  of  D,  the  actual 
as  well  as  derived  stress-strain  curves  show  higher  strain  hardening  behavior  irrespective 
of  the  cavity  expansion  rate.  Even  though  at  smaller  strain  levels,  the  larger  D  curves 
show  lesser  strength  than  the  stress-strain  curves  with  smaller  D  values,  they  reach 
essentially  the  same  strength  at  larger  strain  levels  (see  Figs.  5. 18-5.21).  For  botli  cavity 
expansion  rate  (0. 1  and  1.0%/min),  the  derived  curves  from  PM  expansion  yield  exactly 
the  same  ultimate  strength,  irrespective  of  different  D  values.  Similar  to  previous 
results,  higher  strain  rate  tests  and  higher  ^  value  curves  give  higher  shear  strength 
(Figs.  5.18-5.21). 

Combined  Influence  of  D  and  0 

For  the  values  obtained  from  triaxial  tests  (0  =  0.10)  and  from  the  empirical 
correlations  (D  =  1/500),  Prapaharan  et  al.  (1989)  studied  the  effect  of  strain  rate  in  the 
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Influence  of  Constant  D  on  Undrained  Shear  Strength  Derived  from 
PMT  at  Strain  Rate  =  0.1%/min  (Beta  =  0.10) 


Strain  Rate=1.0%/min,  Beta=0.10,  a=0.01 
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Influence  of  Constant  D  on  Undrained  Shear  Strength  Derived  from 
PMT  at  Strain  Rate  =  1.0%/min  (Beta  =  0.10) 
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Influence  of  Constant  D  on  Undrained  Shear  Strength  Derived  from 
PMT  at  Strain  Rate  =  0.1%/min  (Beta  =  0.15) 


Cavity  Strain  % 

Influence  of  Constant  D  on  Undrained  Shear  Strength  Derived  from 
PMT  at  Strain  Rate  =  1.0%/min  (Beta  =  0.15) 
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Fig.  5.22  Influence  of  D  and  Beta  on  Undrained  Shear  Strength  Derived  from 
PMT  at  Strain  Rate  =  0.1%/min 


Strain  Rate  =  l.C%/min,  Alpha =0.01%/min 
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Fig.  5.23  Influence  of  D  and  Beta  on  Undrained  Shear  Strength  Derived  from 
PMT  at  Strain  Rate  =  1.0%/min 
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pressuremeter  expansion  curves  and  derived  stress-strain  curves  (Figs.  5.11  and  5.15). 
From  the  laboratory  tests  performed  on  kaolin-silica  soil  simulating  pressuremeter  test 
condition  in  CSD  (Chapter  4),  0  was  found  as  0.15  and  D  was  determined  as  1/160. 
Using  these  pairs  of  results,  stress-strain  curves  are  produced  for  cavity  expansion  rate 
of  0.1%/min  and  1.0%/min  (Figs.  5.22  and  5.23).  Surprisingly,  for  the  strain  rate 
0. 1  %/min  case  the  ultimate  strength  are  the  same  for  the  actual  and  derived  curves  (see 
Fig.  5.22).  It  was  found  earlier  that  increasing  0  values  induced  more  strain  softening, 
while  increasing  D  values  caused  for  more  strain  hardening.  Hence,  when  combined  in 
this  case,  the  influence  of  0  and  D  compensated  or  nullified  their  softening/hardening 
behavior  and  produced  the  same  ultimate  strengths  (for  D  =  1/500,  0  =  0.10  and  D  = 
1/160,  0  =  0.15  in  both  actual  and  PMT  cases).  However,  in  the  case  of  larger  cavity 
expansion  rate  (1  %/min),  the  strain  hardening  effect  is  more  than  the  softening  effect, 
and  the  curves  produced  from  the  laboratory  PM  test  data  show  a  slight  increase  in  the 
ultimate  strength  (Fig.  5.23).  It  can  be  concluded  that  the  influence  of  D  and  0  which 
cause  hardening  and  softening,  respectively,  compensated  their  effects  and  the  ultimate 
strengths  were  essentially  the  same  from  the  triaxial  tests/empirical  correlations  and  from 
the  PMT  performed  in  CSD  in  kaolin-silica  soil.  Although  this  conclusion  shed  light  on 
the  reason  for  difference  between  triaxial  and  PMT  based  results,  its  generalization  is  not 
suggested  at  this  stage  as  different  soils  will  have  different  D  and  0  values,  and  their 
effects  may  not  always  nullify  each  other. 

5.11  Conclusions 

The  numerical  integration  performed  to  simulate  the  pressuremeter  expansion 
curve  and  to  derive  the  stress-strain  curve  was  used  to  study  the  influences  of  cavity 
expansion  (strain)  rate,  0  and  D.  Furthermore,  the  effect  of  upper  yield  in  the  prediction 
of  shear  strength  was  also  analyzed,  and  it  was  concluded  that  the  level  of  upper  yield 
has  no  significant  influence  in  the  derived  stress-strain  curve. 

As  the  cavity  expansion  rate  increases  the  shear  strength  also  increases  in  both 
actual  and  PMT  condition  (Fig.  5.15).  Furthermore,  for  strain  hardening  soils  the  PMT 
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give  strain  softening  type  stress-strain  curves.  As  increases  (for  pressuiemeter  test 
condition),  even  though  the  actual  material  curves  show  higher  strength  and  higher 
hardening,  the  stress-strain  curves  derived  from  PM  expansion  curves  show  higher 
softening.  This  parametric  study  which  was  carried  out  using  the  results  for  PMT 
condition,  confirms  the  findings  of  Prevost  (1976)  and  Jamiolkowski  and  Lancellotta 
(1977),  and  demonstrates  that  the  stress-strain  curves  derived  from  the  PMT  will  exhibit 
more  strain-softening  characteristics  under  constant  strain  rate  conditions  even  if  the  soil 
is  strain  hardening. 

The  studies  based  on  the  laboratory  test  data  shown  that,  as  D  increases,  although 
the  curves  have  smaller  initial  slopes,  they  show  strong  strain  hardening  behavior. 
However,  at  large  strains,  the  ultimate  strength  reaches  the  same  level  irrespective  of  the 
values  of  D.  As  far  as  the  ultimate  strength  is  concerned,  the  constant  D  does  not  play 
a  major  role  except  for  the  initial  portion  of  the  curve.  On  the  other  hand,  the  degree 
of  strain  softening  slightly  increases  with  decreasing  values  of  D. 

For  the  soils  studied,  the  combined  effect  of  P  and  D  compensate  their  strain 
softening/hardening  behavior  and  the  ultimate  strength  are  essentially  equal  for  both  the 
triaxial  test  and  pressuremeter  test  values. 


CHAPTER  6 


STRESS  RELIEF  AND  STRESS  RELAXATION  EFFECTS 
IN  PRESSUREMETER  TESTING 


6.1  Introduction 

The  knowledge  of  time  dependent  deformation  behavior  of  soils  is  important  to 
understand  the  stress-strain-time  effects  in  soils  such  as  strain  rate  effect,  creep,  creep 
rupture,  strength  after  creep,  stress  relaxation,  strain  hardening,  strain  softening  and  long 
term  strength.  The  following  time  dependent  behaviors  related  to  pressuremeter  tests 
will  be  discussed  in  this  report;  strain  rate  effect  (already  introduced  in  Chapters  4  and 
5),  stress  relief,  stress  relaxation  time  (or  normalizing  period),  creep,  and  stress 
relaxation. 

In  the  Menard  type  pressuremeter  or  even  the  self-boring  pressuremeter  the 
borehole  is  slightly  larger  than  the  probe.  Hence,  before  the  membrane  is  inflated  the 
soil  around  the  probe  tends  to  move  inward  and  All  the  gap.  During  this  process  of 
stress  relief  the  lateral  stresses  are  released  and  the  in  situ  horizontal  stress  is  reduced. 
In  the  SBPM,  the  influence  of  the  size  of  cutting  shoe  on  the  stress  relief  is  also 
important. 

The  stress  relaxation  time  is  the  period  usually  allowed  between  the  end  of 
drilling  and  beginning  of  probe  expansion.  This  time  is  provided  to  allow  for  dissipation 
of  the  excess  pore  pressure  developed  during  the  drilling  operation.  The  relaxation  time 
depends  on  the  type  of  soil  and  the  drilling  method.  In  the  laboratory  simulation  of  PMT 


in  the  CSD,  relaxation  tests  were  performed  with  various  strain  levels  to  study  the 
influence  of  initial  strains  on  the  relaxation  time. 

Most  material  masses,  when  subjected  to  a  sustained  loading,  deform  or  creep 
continually,  i.e.  their  dimensions  change  with  time.  Conversely,  if  a  material  is 
deformed  to  a  certain  amount  and  then  its  dimensions  are  maintained  unchanged,  in 
general,  a  time-dependent  decreasing  stress,  acting  in  the  direction  whose  dimension  is 
maintained,  will  be  required.  This  phraomenon  is  commonly  referred  to  as  "stress 
relaxation".  "Creep"  can  be  defined  as  the  deformation  of  a  soil  mass  under  a  constant 
stress  state  and  constant  volume  or  pore  pressure.  Stress  relaxation  can  be  considered 
as  the  reverse  of  creep,  i.e.  stress  change  while  keeping  the  same  dimensions.  In  simple 
terms,  creep  is  pressure  controlled  and  relaxation  is  volume  controlled. 

6.2  Stress  Relief  in  Pressuremeter  Test 

One  of  the  major  advantages  of  the  in  situ  testing  over  the  conventional  laboratory 
tests  is  that  there  is  an  opportunity  to  prevent  or  limit  sample  disturbance  due  to  stress 
relief.  However,  in  the  case  of  Menard  pressuremeter  test,  the  probe  has  to  be  inserted 
in  a  predrilled  hole.  During  the  time  between  drilling  and  probe  expansion,  the  borehole 
wall  is  subjected  to  stress  relief.  Even  in  the  SBPM,  the  cutting  shoe  is  designed  such 
that  it  makes  a  hole  slightly  larger  than  the  membrane  in  order  to  reduce  friction  between 
the  probe  and  wall  to  decrease  the  disturbance  to  the  wall  surface  and  to  avoid  damage 
to  the  membrane.  If  the  borehole  is  larger  than  the  probe,  then  a  gap  is  created  between 
the  soil  and  the  membrane,  thus  causing  stress  release  (relief)  and  reduction  of  the  in  situ 
horizontal  stress.  This  unloading  is  represented  by  path  AB  in  Fig.  6.1,  where  point  A 
corresponds  to  the  in  situ  condition  (Law  and  Eden,  1980).  Reloading  by  inflating  the 
membrane  is  indicated  by  path  BC  where  C  represents  the  condition  when  the  inflated 
membrane  diameter  is  exactly  equal  to  that  of  the  cutting  shoe.  That  is,  the  clay  is  now 
pushed  back  to  the  original  location  before  the  boring.  Further  increasing  the  probe 
pressure  leads  to  point  D  beyond  which  the  path  will  be  identical  to  the  ideal  test.  If  the 
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cutting  shoe  is  smaller  than  the  membrane,  then  the  clay  has  to  be  pushed  laterally  to 
make  room  for  the  membrane,  thus  introducing  a  load  on  the  membrane  which  is 
represented  by  the  path  A£  in  the  figure.  If  the  test  is  allowed  for  a  stress  relaxation 
process,  the  pressure  will  drop  to  F  (path  £F),  where  the  inflation  of  the  membrane 
begins  and  joins  the  ideal  loading  path  at  G.  The  corresponding  interpreted  shear 
stresses  for  ideal  and  oversized  cutting  shoe  tests  are  given  in  Fig.  6.2.  The  oversized 
cutting  shoe  leads  to  stress  relief  and  overestimation  of  the  undrained  shear  strength. 
Law  and  Eden  (1980)  determined  that  stress  relief  caused  by  an  oversized  cutting  shoe 
in  a  SBPM  overestimated  the  deduced  shear  strength  by  IS  or  80%  depending  on  the 
choice  of  the  zero  reference  strain  (point  C  or  B  in  Fig.  6.1).  Furthermore,  the  derived 
modulus  was  also  overestimated  by  30%  due  to  unloading  prior  to  shear. 

Benoit  and  Clough  (1986)  conducted  full  scale  SBPM  test  in  soft  San  Francisco 
Bay  Mud  to  study  the  effect  of  oversized  cutting  shoe  on  various  parameters.  It  was 
observed  that  the  lateral  pressures  obtained  were  consistently  around  20%  lower  for  the 
1.1%  oversized  shoe  than  for  the  perfectly  fitting  shoe.  The  reason  for  the  lower  lateral 
pressure  is  obviously  the  stress  relief  caused  by  the  oversized  cutting  shoe.  This  is 
supported  by  the  data  of  Hughes  et  al.  (1980),  where  the  in  situ  horizontal  stresses  were 
reduced  by  60  to  65%  in  the  Coode  Island  silty  clay  when  a  3%  oversized  cutting  shoe 
was  used.  Similarly,  in  stiff  Leda  clay.  Law  and  Eden  (1980)  found  the  lateral  pressure 
to  be  underestimated  by  about  40-45%  with  a  1.1%  oversized  cutting  shoe.  Benoit  and 
Clough  (1986)  also  determined  from  the  SBPM  tests  that  the  1 . 1  %  oversized  cutting  shoe 
tests  led  to  an  overestimation  of  60-100%  in  undrained  shear  strength,  which  is  consistent 
with  the  80%  increase  in  shear  strength  reported  by  Law  and  Eden  (1980). 

All  these  results  show  that  the  stress  relief  occurrs  in  predrilled  boreholes  and  in 
holes  drilled  by  self  boring  pressuremeter  with  an  oversized  cutting  shoe,  and  causes 
unloading  of  the  soil  around  the  borehole  and  leads  to  underestimation  of  the  in  situ 
horizontal  pressure.  From  these  experimental  evidence,  the  importance  of  stress  relief 
in  PMT  is  evident  and  it  should  be  taken  care  in  experimental  and  analytical  studies. 
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6.3  Stress  Relaxation  Time  in  Pressuremeter  Test 

The  stress  relaxation  time,  also  called  the  normalization  period,  is  the  time 
interval  between  the  end  of  self-boring  and  the  beginning  of  the  probe  expansion.  During 
this  period,  the  excess  pore  pressure  developed  during  the  drilling  operation  must  be 
allowed  to  dissipate,  and  the  disturbed  stress  state  around  the  probe  will  reach 
equilibrium,  i.e.  return  to  its  original  state.  This  time  period  is  dependent  upon  the 
drilling  technique  used  and  the  type  of  soil  tested,  and  it  varies  from  a  few  minutes  to 
several  hours. 

Jamiolkowski  and  Lancellotta  (1979)  showed  that  for  Porto  Tolle  clay,  a 
relaxation  time  of  about  100  minutes  was  enough  to  determine  the  proper  total  in-situ 
horizontal  stress.  For  two  Norwegian  clays,  Lacasse  et  al.  (1981)  i^x>rted  a  relaxation 
time  between  2  and  22  hours.  Lacasse  and  Lunne  (1982)  indicated  that  for  a  relaxation 
period  varying  from  90  to  1300  minutes  for  Drummen  clay  there  was  little  effect  on  the 
total  in-situ  lateral  stress  determination.  Denby  (1978)  showed  that  higher  lift-off 
pressures  were  obtained  when  a  test  was  performed  with  a  relaxation  time  of  30  minutes 
compared  to  results  from  tests  using  a  period  of  120  minutes  or  more.  However,  he 
obtained  identical  shear  strengths,  irrespective  of  the  different  relaxation  times  used. 
Denby  used  a  normalization  period  of  at  least  120  minutes  in  the  San  Francisco  Bay  Mud 
tests.  In  a  Tokyo  Bay  soft  clay,  Mori  (,*981)  indicated  that  a  120  minute  period  was 
necessary  to  bring  down  the  excess  pore  pressures  to  less  than  0.75  psi.  The  SBPM  tests 
carried  out  in  Boston  Blue  clay  after  permitting  only  a  ten  to  thirty  minute  equilibrium 
period  (Ladd  et  al.,  1980)  led  to  inconsistent  values  of  lateral  pressures.  From  the 
SBPM  tests  conducted  in  San  Francisco  Bay  Mud,  Benoit  (1983)  found  that  normalization 
periods  of  90  minutes  or  more  did  not  influence  the  results.  A  summary  of  several  tests 
which  used  various  relaxation  times  is  shown  in  Table  6.1.  From  these  results  it  can  be 
concluded  that  most  of  the  excess  pore  pressures  induced  from  drilling  are  dissipated  in 
about  two  hours  for  typical  soft  clays. 


Table  6. 1  Relaxation  Time  for  Several  Field  Pressuremeter  Tests 

(After  Benoit,  1983) 


TEST  SITE 

NORMALIZATION 

EXPANSION  RATE 

OR 

SOIL  TYPE 

DCO  Tnn 

rtKiUU 

(minutes) 

psi/min 

S/min 

REFERENCES 

Porto  Tolle  clay 

15-120 

0.83 

Jamiol kowski 

Guasticce  clay 

15-120 

Bjfl 

0.83 

and  Lancellotta 
(1977) 

urammen  clay 

90-1300 

Lacasse  and 

Onsoy  clay 

90-1300 

— 

- — 

Lunne 

(1982) 

Hamilton  Air 

Force  Base: 

■Oenby 

180 

1.0 

Denby  (1978) 

Benoi t 

60-8460 

0.24-6.9 

0.04-0.99 

Benoit  (1983) 

Boston  Blue  clay 

generally  10-30 

1.0 

Ladd  et  al . 

others  42-720 

(1980) 

South  Gloucester 

60-1080 

1.4 

Law  and  Eden 

Matagami 

60-1080 

1.4 

— 

(1980) 

Coode  Island 

60-960 

WWW 

Hughes  et  al . 

(1980) 

Tokyo  Bay 

120 

— 

Mori  (1981) 

162 


For  clays,  the  total  horizontal  stress  during  self-boring  is  higher  than  the  in  situ 
value.  During  relaxation,  it  decreases  gradually  and  becomes  asymptotic  to  the  in  situ 
horizontal  stress  (PJ.  In  the  case  of  loose  sands  or  silts,  the  lateral  stress  decreases  and 
again  increases  to  the  in  situ  horizontal  stress  (Baguelin  et  al.,  1974).  Too  short  a 
relaxation  time  usually  leads  to  an  incorrect  estimate  of  (Dmby  1978,  and 
Jamiolkowski  and  Lancellotta  1979),  and  the  effect  seems  to  be  larger  as  depth  increases. 
This  is  confirmed  by  results  from  tests  in  Boston  Blue  clay  that  showed  a  total  in-situ 
lateral  stress  which  was  consistently  underestimated  using  normalization  period  of  10-30 
minutes  (Ladd  et  al.,  1980).  In  other  cases,  the  values  were  overestimated  for  short 
normalization  periods  in  the  Porto  Tolle  clay  and  for  long  normalization  periods  in  the 
Onsoy  clay  (Jamiolkowski  and  Lancellotta  1977,  and  Lacasse  and  Lunne  1982).  Benoit 
(1983)  observed  relaxation  periods  varying  from  1  to  22  hours  with  Bay  Mud  with  a 
larger  cutting  shoe  having  no  effect  on  the  measured  lateral  pressure. 

Selection  of  P,,,  the  datum,  has  significant  influence  on  the  stress-strain  curve  and 
the  undrained  strength  derived  from  the  pressuremeter  tests  (Ladd  et  al.,  1980). 
Therefore,  a  sufficient  relaxation  time  should  be  allowed  for  the  correct  determination 
of  Po.  Baguelin  et  al.  (1978)  recommend  that  the  test  be  started  if  the  change  in  pressure 
over  a  period  of  10  minutes  is  less  than  O.IS  psi  (1  kPa). 

6.4  Pressuremeter  Creep  Test 

Many  soils  -  clays,  silts,  and  all  frozen  soils  -  have  rheological  properties,  i.e., 
the  ability  to  develop  creep  deformations  and  reduce  strength  under  sustained  stresses. 
The  time-dependent  deformation  behavior  of  clays  under  continued  loading  depends  upon 
several  variables  such  as  soil  type,  soil  structure,  stress  history,  effective  stress, 
temperature,  etc. 

There  has  been  extensive  research  performed  to  characterize  the  creep  and 
relaxation  behavior  of  clays  using  rheological  models  composed  of  linear  springs  in 
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combination  with  linear  or  nonlinear  dashpots.  Singh  and  Mitchell  (1968)  using  a  rate 
process  theory  and  published  experimental  results,  formulated  a  stress-strain-time 
relationship  in  which  the  creep  rate  varied  exponentially  with  stress  level,  and  decreased 
nonlinearly  with  time.  Haefeli  (1965)  assumed  a  creep  law  for  snow,  ice,  rock  and  soil 
in  which  the  creep  rate  varied  nonlinearly  with  stress  but  did  not  vary  with  time  or 
strain.  Prevost  (1976)  developed  a  phenomenological  equation  to  describe  the  stress- 
strain-time  behavior  of  normally  consolidated  clays  loaded  under  undrained  conditions. 

Most  of  the  above  research  is  based  on  triaxial  testt.  Very  few  investigators  used 
pressuremeters  to  study  the  creep  behavior  of  soil.  Ladanyi  (1982)  and  his  collaborators 
(for  example,  Ladanyi  and  Johnston,  1978)  have  used  the  pressuremeter  to  study  the 
creep  potential  of  frozen  soils.  The  volume  of  the  pressuremeter  cell  is  a  limitation  of 
pressuremeter  creep  tests.  Therefore,  PM  creep  test  is  feasible  only  for  short  creep  times 
and  a  medium  range  of  stress.  To  obtain  the  long-term  creep  parameters,  Ladanyi  et  al. 
(1991)  used  borehole  relaxation  tests  as  the  alternative.  The  main  advantage  of  the 
relaxation  test  is  that  the  strain  is  controlled  and  the  stress  variation  is  observed  so  that 
there  is  no  danger  of  exceeding  the  volume  capacity  of  the  PM  cell.  As  a  consequence, 
borehole  relaxation  tests  can  cover  easily  the  whole  range  of  stresses  and  long  periods 
of  time. 

Ladanyi  (1982)  derived  the  following  expressions  to  determine  creep  parameters 
from  Menard  PMT.  The  total  strain  attained  after  a  given  time  under  constant  stress  is: 

€*  =  fe  +  e'e  6.1 

where  6^  is  the  instantaneous,  not  necessarily  elastic,  portion  of  the  total  strain,  and  €% 
is  the  creep  strain  given  by: 


<  ■  uyirtoyo/  »* 


6.2 
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where  the  subscript  e  denotes  the  von  Mises  equivalent  stress  and  strain,  is  the  creep 
modulus  corresponding  to  an  arbitrary  strain  rate  t  is  the  time,  and  b  and  n  are 
creep  exponents. 

When  Eq.  6.2  is  applied  to  the  problem  of  an  expanding  cylindrical  cavity  in  an 
infinite  medium,  originally  acted  upon  by  an  isotropic  lateral  stress  Po,  the  radial  creq) 
strain  rate  under  a  constant  stress  p,  >  Po  can  be  derived  as: 

drj/dt  =  ri  F  b  t^‘  6.3 

where  F  is  a  function  of  (p,-Po),  given  by: 

F  =  (M/2)[(prPo)/aJ“ 
and 

M "  2(j3nr' 

In  a  stage-loaded  creep  test,  if  p;  is  the  stress  applied  in  the  borehole  during  the  stage  k, 
following  a  smaller  stress  in  the  previous  stage  (k-1),  the  resulting  radius  increase  with 
time  is  given  by: 

ln(r/ri^-,)  =  Ft*-  6.6 

or,  in  terms  of  the  borehole  volume: 

In  (V/Vk.,)  =  2  Ft**  6.7 

where  V  =  rr^iL  is  the  current  volume  of  the  cavity  of  length  L. 


6.4 

6.5 


Finally,  the  value  of  creep  modulus  is  given  by: 


<r.  =  (PrPjN  [M/(2F)h]''“ 


6.8 


where  M  is  given  by  Eq.  6.S  and  (2F)n  denotes  the  value  of  (2F)  at  an  arbitrary  point 
(Pi'PJnj  located  on  a  (2F)  versus  (pj-pj  straightline  segment  in  a  log-log  plot  (Fig.  6.3). 

In  section  6.6  this  approach  will  be  combined  to  relaxation  tests  to  determine 
creep  parameters. 

6.5  Stress  Relaxation  Parameters  from  Triaxial  Tests 

The  creep  behavior  of  soils  has  been  studied  primarily  from  the  point  of  view  of 
deformation.  There  are,  however,  problems  in  which  the  deformations  are  negligible, 
and  the  prediction  of  the  stresses  acting  on  a  structure  due  to  interaction  with  a  soil  mass 
(for  example  a  retaining  wall)  is  of  primary  interest.  This  type  of  test,  which  studies  the 
variation  of  stresses  while  the  deformation  is  kept  constant  is  called  stress  relaxation  test. 

Relatively  few  researchers  have  studied  stress  relaxation  in  soils.  Murayama  and 
Shibata  (1961)  that  the  decay  of  deviator  stress  with  the  logarithm  of  time  was  linear,  up 
to  a  certain  limit,  and  then  remained  constant.  Vialov  and  Skibitsky  (1961)  also  noticed 
the  linear  decrease  of  stress  with  logarithm  of  time  but  they  did  not  find  the  existence  of 
a  final  relaxed  level  of  stress.  Probably  their  relaxation  time  was  not  sufficient  as  the 
tests  lasted  only  4  hours.  Saada  (1962)  obtained  a  linear  decay  stress  with  logarithm  of 
time,  up  to  about  SO  days,  and  then  the  deviatoric  stress  abruptly  fell  to  zero.  From  the 
stress  relaxation  tests  carried  out  on  undisturbed  Sault  St.  Marie  clay,  Christensen  and 
Wu  (1964)  also  obtained  linear  relation  between  stresses  and  logarithm  of  time. 
Murayama  et  al.  (1974)  and  Akai  et  al.  (1975)  presented  similar  experimental  results. 

Lacerda  and  Houston  (1973)  and  Lacerda  (1976)  have  thoroughly  studied  the 
stress  relaxation  and  creep  effects  on  soil  deformation  by  performing  several  tests  on 
undisturbed  San  Francisco  Bay  soft  clay,  remolded  kaolinite,  clean  quartz  sand  and 
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6.3  Cr^  Pai^eter  Determination  from  Stage-Loaded  Pressuremeter  Test 
(After  Ladar«'n,  1982) 
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compacted  clay.  A  summary  of  the  stress  relaxation  tests  on  various  soils  is  shown  in 
Fig.  6.4.  These  results  show  that  when  the  stress  relaxation  process  is  started  by 
confining  the  deformations,  the  shear  stress  starts  to  relax  significantly  after  a  ^te  time 
elapses  and  the  stress  relaxation  is  solely  dq)endent  on  the  strain  rate  and  the  initial  stress 
level.  Consistently  with  previous  data,  the  stress  relaxation  was  found  to  have  an  inverse 
linear  dependence  with  the  logarithm  of  time.  The  ratio  of  the  deviator  stress  at  any  time 
to  the  stress  at  the  beginning  of  stress  relaxation  is  linear  with  the  logarithm  of  time,  and 
the  slope  is  independent  of  the  confining  pressure,  strain,  and  strain  rate  but  dependent 
on  initial  stress  level.  It  is  observed  that  the  variation  of  pore  pressure  was  very  small 
in  the  undrained  stress  relaxation  tests.  Similarly,  the  volume  change  during  relaxation 
in  the  drained  tests  was  also  negligible.  Furthermore,  there  was  no  difference  between 
the  stress  relaxation  of  anisotropically  and  isotropically  consolidated  soils. 

Lacerda  (1976)  derived  the  following  equation  for  stress  relaxation  from  Singh- 
Mitchell’s  creep  equation  (1968): 

«  =  =  1  -  J  log(-i),  for  t  >  t,  6.9 

O.  D,  K 

where  s  =  slope  of  stress  relaxation  curve  =  (shown  in  Fig.  6.5),  and  4>  is  defined 
as: 


♦  = 


2.3(1 -w) 
a 


6.10 


D  =  deviator  stress 
D  =  deviator  stress  level  =  D/D^ 

Do  =  deviator  stress  at  time  to 

Do  =  deviator  stress  level  at  time  to  =  Do/Do  „„ 

to  =  time  at  beginning  of  stress  relaxation 


The  value  of  4  is  obtained  from  the  intercept  of  the  straight  line  portion  of  the 


relaxation  curve  on  the  time  axis  of  the  stress  level-log  time  plot.  It  is  dependent 
exclusively  on  soil  type  and  strain  rate  prior  to  the  relaxation  test.  The  higher  the  strain 
rate,  the  smaller  is  the  value  of  tg. 

The  above  equations  are  valid  for  m  <  1  and  they  express  the  relationship 
between  q,  the  relaxed  stress  relative  to  the  initial  stress,  and  time,  t,  which  may  be 
derived  by  inverting  Singh  and  Mitchell’s  (1968)  creep  equation; 

c  =  i4  exp(a5)  (tjty  6.11 

where  €  =  strain  rate  (96/min),  D  =  D/Dg,,,  (stress  level),  a  =  slope  of  the  linear 

portion  of  the  logarithm  strain  rate  versus  deviator  stress  plot  obtained  from  creep  test 
(in  units  1/pressure),  5  =  a  (dimensionless  parameter),  m  =  absolute  value  of 
the  slope  of  the  straight  line  relationship  between  the  logarithm  of  strain  rate  and  the 
logarithm  of  time,  t,  =  unit  time,  e.g.,  1  minute,  and  A  =  extrapolated  value  of  strain 
rate  for  zero  stress  level  in  the  linear  plot  of  log  strain  rate  versus  deviator  stress,  for 
time  ti,  in  %/min. 

In  Eq.  6.9  tg  can  be  related  to  strain  rate  by: 

t,  =  hj€  6.12 

where  hg  is  the  numerical  value  of  the  strain  rate  (in  units  of  strain)  necessary  to  yield 
a  "delay  time"  of  t„  before  stresses  begin  to  relax  (shown  in  Fig.  6.6).  The  values  of  s, 
the  slopes  of  the  stress  relaxation  curves,  were  found  to  be  independent  of  confining  or 
consolidation  pressure  and  of  initial  relative  density  of  sand. 


The  relationship  between  A  (creq>)  and  hg  (relaxation)  was  derived  from  Eq. 

6.11: 


hg  =  13.2  A“ 


6.13 


From  several  published  and  his  own  experimental  data,  Lacerda  (1976)  obtained 


Log  Strain  Rate  Versus  Log  Time  Plot  (After  Lacerda  and  Houston, 
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the  following  correlations  between  the  stress  relaxation  parameters  and  plasticity  index 
(PD: 

^ - EL -  +  0.022  6.14 

4.4  P/+280 

where  PI  is  in  %,  and 

log  h,  =  0.0285  (PD  -  2.82  6.15 


These  results  can  be  used  to  obtain  the  stress  relaxation  and  creep  parameters  and 
incorporated  in  the  time  dependent  plasticity  models  to  predict  the  stress-strain-time 
behavior  of  soils. 


6.6  Pressuremeter  Stress  Relaxation  Test 


Unlike  the  stress  relaxation  studies  carried  out  by  many  researchers  using  the 
triaxial  apparatus,  only  Ladanyi  and  his  collaborators  (Ladanyi  et  al.,  1991)  have  used 
the  pressuremeter  to  obtain  relaxation/creep  parameters.  Because  of  the  volume 
limitations  of  the  pressuremeter  cell,  borehole  relaxation  tests  are  performed  as  an 
alternative  to  borehole  creep  test  to  obtain  the  long-term  creep  information.  In  related 
experiments,  Ladanyi  (1982)  used  the  dilatometer  to  perform  stress  relaxation  tests,  and 
Ladanyi  and  Huneault  (1989)  obtained  creep  parameters  of  frozen  soil  with  the  cone 
penetrometer. 

Using  the  aging  (time  hardening)  theory  of  creep,  which  assumes  that  there  is  a 
unique  and  continuous  surface  in  space  relating  stress  with  strain  and  time,  Ladanyi 
(1982)  derived  expressions  for  stress  relaxation  from  which  creep  parameters  could  be 
inferred.  In  the  aging  theory,  it  is  assumed  that  creq>  and  relaxation  are  closely  related, 
so  that  a  relaxation  curve  is  nothing  else  but  a  creep  curve  under  a  continuously  decaying 
stress,  resulting  in  a  constant  value  of  strain.  In  other  words,  according  to  this 
assumption,  any  constitutive  creq)  relation  can  be  directly  transformed  into  a  relaxation 
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relation  by  making  the  creep  strain  constant  and  equal  to  the  applied  initial  strain. 


According  to  the  aging  creep  theory  and  using  an  approximate  method,  Ladanyi 
(1982)  presented  the  equation  for  a  family  of  borehole  relaxation  curves  in  an  expanding 
cylindrical  cavity  problem: 


wy;,) 

Af(f'+r)* 


where  M  is  given  by  £q.  6.S,  Pj  is  the  internal  pressure  applied  in  the  borehole,  p^  is  the 
total  lateral  ground  stress  around  the  borehole,  and  t’  is  a  very  short  time  interval.  If  t’ 
is  neglected,  and  the  relaxation  curves  are  plotted  in  a  log  (p.-pj  versus  log  t  plot  (Fig. 
6.7),  with  the  strain  In  (V/VJ  as  the  parameter,  then  their  slope  at  the  end  of  interval 
gives  the  ratio: 

b  ^  _  Alog(p,-/>o)  ^  V  617 

n  Alogr  h 

with  V  and  h  given  in  Fig.  6.7.  For  the  same  time  interval  (for  the  same  relaxation 
period),  where  t  =  tj  =  constant,  from  the  log  nn(V/Vo)]  versus  log  (pj-pJ  plot,  the 
slope  of  the  curve  is: 

„  .  =  h  6.18 

Alog(p,-p^  V, 

with  V,  and  h,  shown  in  Fig.  6.5 


When  b  and  n  are  known,  the  value  of  the  creep  modulus  for  a  given  can 
be  calculated  for  any  point  k  on  the  line  In  (V/VJ  versus  (Pj-pJ  by  using  Eq.  6.16: 


=  (PrPJk 


When  Eq.  6.16  is  applied  to  two  consecutive  points  in  a  relaxation  line  with 


coordinates  (ti,  p^)  and  (t2,  Pjq))  the  value  of  in  situ  horizontal  pressure  Po  can  be 
determined  from: 


6.20 


«  5  21 

It  is  possible  to  perform  the  stress  relaxation  tests  with  the  model  pressuremeter 
in  Calibration  Chamber  and  can  obtain  the  creep  parameters.  The  in  situ  horizontal 
stress  could  be  computed  using  the  Eq.  6.21  and  be  compared  with  the  Po  obtained  from 
the  pressuremeter  expansion  curve. 

6.7  Pressuremeter  Stress  Relaxation  Test  in  CSD 

To  the  author’s  knowledge  there  has  not  been  any  systematic  studies  (field  tests 
or  laboratory  simulations)  conducted  to  evaluate  the  relaxation  effects  (relaxation  time 
and  stress  relaxation)  in  pressuremeter  testing.  So  far,  only  Ladanyi  et  al.  (1991) 
performed  field  pressuremeter  stress  relaxation  tests  to  determine  creep  parameters  of 
frozen  soils. 

The  effect  of  stress  relief  due  to  oversize  cutting  shoe  on  in  situ  horizontal  stress 
and  undrained  shear  strength  prediction  was  already  discussed  in  section  6.2.  From  the 
field  SBPM  tests,  it  was  estimated  that  lateral  pressures  were  20-45%  lower  and 
undrained  shear  strengths  60-100%  higher  for  a  1. 1  %  oversized  borehole  (Law  and  Eden 
1980,  and  Benoit  and  Clough  1986),  and  60-65%  lower  in  situ  horizontal  stress  was 
predicted  for  a  3%  oversized  borehole  (Hughes  et  al.,  1980).  However,  there  was  no 
attempt  made  to  correlate  the  stress  relief  (reduction  of  in  situ  horizontal  stress)  and  the 
stress  relaxation  or  relaxation  time  (normalizing  period).  Table  6.1  gave  the  relaxation 
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times  adopted  by  various  researchers  but  they  were  not  correlated  to  the  initial  strain  or 
to  the  initial  strain  rate  which  can  be  approximated  to  the  cutting  rate. 

Therefore,  a  laboratory  simulation  study  was  undertaken  to  attempt  to  correlate 
the  stress  relief,  stress  relaxation  time  and  stress  relaxation  with  the  initial  amount  of 
strain  and  strain  rate  in  a  pressuremeter  test.  The  simulation  of  pressuremeter  test  in  the 
cuboidal  shear  device  has  been  already  described  in  Chapter  4.  Samples  were  prepared 
from  artificially  sedimented  kaolin  clay  as  described  in  Chapter  3.  All  the  specimens 
were  prepared  and  one  dimensionally  consolidated  in  exactly  the  same  way  as  for  the 
strain  rate  tests. 

Pressuremeter  stress  relaxation  tests  were  simulated  in  the  cuboidal  shear  device 
as  follows:  After  the  one-dimensional  consolidation,  the  drainage  valve  was  closed  and 
was  increased  such  that  the  specimen  is  deformed  at  the  specified  strain  rate,  varying 
from  0.005  %/min  to  0.1%/min  -  depending  on  the  test,  up  to  the  required  amount  of 
strain.  The  strain  rates  and  strains  before  relaxation  are  tabulated  in  Table  6.2.  Since 
the  sample  was  undergoing  a  pressuremeter  stress/strain  path,  a,  was  adjusted  such  that 
no  vertical  deformation  was  permitted  (plane  strain  condition  in  vertical  direction).  The 
value  of  Oy  was  adjusted  such  that  6,  =  -£y,  to  assume  no  volume  change  (undrained 
condition).  Having  achieved  the  required  amount  of  strain,  the  sample  was  permitted  to 
relax.  The  relaxation  condition  was  imposed  on  the  specimen  by  adjusting  all  three 
pressures  to  maintain  zero  deformation  in  all  three  directions.  A  tolerance  of  0.0(X)1" 
was  permitted  in  the  computer  control. 

In  actual  practice,  when  the  soil  around  the  PM  probe  is  relaxing  there  is  a 
possibility  for  drainage  because  the  soil  is  permitted  to  relax  for  a  normalization  period 
of  120  min  to  several  hours.  Therefore,  during  the  relaxation  period  the  drainage 
condition  is  neither  fiilly  undrained  nor  fully  drained.  In  order  to  cover  both  extremes, 
relaxation  tests  were  performed  in  both  undrained  and  drained  conditions. 


Table  6.2  Relaxation  Time  for  Differait  Stram  Levels  in  Pressuremeter  Test 


Test  No. 

Strain  Rate 

(%  per  min.) 

Strain  Before 
Relaxation 
(%) 

Relaxation 

Time 

(min.) 

87-U 

0.1 

1.0 

200 

89-U 

0.1 

1.0 

200 

90-U 

0.1 

1.0 

250 

93-U 

0.05 

0.5 

250 

94-U 

0.05 

0.5 

200 

95-U 

0.05 

0.5 

200 

88-U 

0.01 

0.1 

225 

91-U 

0.01 

0.1 

200 

‘99-D 

0.01 

0.1 

250 

100-D 

0.005 

0.1 

200 

U  -  Undrained  Test 


D  -  Drained  Test 
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6.8  Results  and  Discussion 

Typical  plots  of  effective  stresses  versus  time  are  shown  in  Figs.  6.8-6.11,  for 
each  strain  level  of  1.0,  O.S,  0.1%,  and  for  undrained  condition  and  drained  conditions. 
Other  plots  are  presented  in  Appendix  C.  Table  6.2  summarizes  the  relaxation  times 
obtained  from  all  the  relaxation  tests  reported  herein  and  in  Appendix  C. 

Irrespective  of  the  initial  strain  level,  the  strain  rate,  or  the  drainage  condition, 
the  relaxation  time  needed  to  take  the  stresses  to  their  initial  (steady)  state  is  in  between 
2(X)  and  250  minutes.  The  relaxation  time  may  depend  on  the  type  of  soil  and  the  type 
of  drilling  technique  used  because  the  pore  pressure  dissipation  time  depends  on  the 
permeability  of  the  soil  and  the  disturbance/pore  pressure  generation  dq)ends  on  the  type 
of  drilling  equipment  used.  However,  for  kaolin  clay,  the  relaxation  time  is  about  four 
hours  regardless  of  cutting  shoe  size  or  the  rate  of  cutting. 

From  the  relaxation  curves,  the  coefficient  of  lateral  earth  pressure  (K,,)  was 
estimated  before  and  after  relaxation,  and  the  percent  of  change  in  was  also  calculated 
and  tabulated  in  Table  6.3.  Again,  the  AK^  is  more  or  less  constant  (varies  between 
6.3%  and  6.5%)  within  the  range  of  strain  tested  (0.1%  to  1.0%  strain  and  0.005%/min 
to  0. 1  %/min  strain  rate).  Therefore,  it  can  be  concluded  that  for  an  oversize  boring  of 
0.1%  to  1.0%  the  stress  relief  causes  an  underprediction  of  lateral  earth  pressure  by 
6.4%  only. 

6.9  Conclusion 

There  are  several  time  related  factors  which  influence  the  measurement  and 
prediction  of  pressuremeter  tests.  Most  importantly,  the  following  time  dependent 
behaviors  are  of  concern  such  as  loading  rate  or  strain  rate,  stress  relief,  creep,  stress 
relaxation,  strain  hardening  and  strain  softening.  The  strain  rate  effect  in  pressuremeter 
testing  was  found  to  be  about  50%  more  than  the  effect  in  triaxial  testing. 


Fig.  6.9  Relaxation  Curves  for  Kaolin  After  0.5%  Strain-Undrained  Test  (No.  94) 


Table  6.3  Change  in  Coefficient  of  Lateral  Pressure  During  Relaxation 
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Stress  relief  and  stress  relaxation  in  PMT  are  important  at  least  when  compared 
to  triaxial  testing.  It  is  commonly  expected  that  the  larger  the  size  of  overcut  larger  the 
amount  of  stress  relief.  However,  the  pressuremeter  stress  paths  simulated  in  the 
cuboidal  shear  device  revealed  that  the  change  in  the  coefficient  of  lateral  earth  pressure 
was  not  affected  by  the  amount  of  strain  (i.e.,  size  of  over  cut).  The  change  in  Ko  was 
found  to  be  constant,  about  6.4%  for  0.1%,  0.5%  and  1.0%  strain. 

Usually  the  stress  relaxation  time  and  the  probe  expansion  rate  are  rarely 
quantified  in  the  literature.  But  these  are  important  parameters  and  little  knowledge  is 
available  to  guide  the  practitioner  as  to  the  proper  value.  Commonly  they  are  referred 
to  in  a  descriptive  manner  such  as  the  stress  relaxation  time  must  be  long  enough  to 
permit  the  excess  pore  pressure  to  dissipate,  and  the  probe  expansion  rate  fast  enough 
to  prevent  drainage.  Unfortunately,  these  parameters  are  not  only  soil  dependent  but  also 
drilling  technique  and  operator  dependent.  To  shed  more  light  into  this  qualitative  nature 
of  defining  stress  relaxation  time,  pressuremeter  tests  were  simulated  in  the  cuboidal 
shear  device  and  the  influence  of  strain  and  strain  rate  on  the  stress  relaxation  was 
studied.  For  kaolin  clay,  in  the  limited  strain  range  of  0.1%  to  1.0%  of  strain  (which 
reflects  the  amount  of  overcutting)  and  for  the  strain  rate  of  0.005  %/min  to  0.1%/min 
(which  is  a  measure  of  drilling  rate),  the  relaxation  time  was  about  200  to  250  minutes. 


CHAPTER? 


CONCLUSIONS 

This  chapter  summarizes  the  main  conclusions  drawn  from  this  investigation. 
Other  results  and  conclusions  are  presented  at  the  end  of  each  chapter. 

1.  The  automated  flexible  wall  calibration  chamber  model  pressuremeter  tests 
performed  on  two  soils,  kaolin  and  kaolin-silica  mixture,  show  a  satisfactory  performance 
of  the  improved  pressure  control  system  (electro-pneumatic  control)  and  the  data 
acquisition  system.  The  newly  designed  piezometer  also  performed  very  well.  It  was 
a  major  problem  in  previous  studies. 

2.  The  systematic  procedure  adopted  in  sample  preparation  yielded  a  very 
consistent  uniform  samples.  The  cuboidal  shear  device  is  a  very  versatile  equipment 
which  can  simulate  any  kind  of  stress  path  with  the  proper  servo-controlled  hardware  and 
software.  The  mechanical  solenoid  valves  hindered  the  speed  of  the  experiment, 
otherwise  all  components  of  the  simulation  system  worked  satisfactorily. 

3.  Pressuremeter  tests  were  simulated  in  the  cuboidal  shear  device  by  applying 
a  shear  load  on  a  one  dimensionally  consolidated  sample  such  that  no  deformation  was 
allowed  in  the  vertical  direction  (plane  strain  condition)  and  allowing  the  expansion  in 
y-  direction  to  be  equal  to  the  compression  in  the  x-direction  (no  volume  change  in  the 
undrained  condition).  Three  or  more  tests  were  performed  at  strain  rates  of  0.01%, 
0.05%,  0.10%,  0.50%,  1.00%,  and  5.00%  per  min.  Excellent  agreement  was  observed 
for  the  tests  performed  at  the  same  strain  rate.  It  was  found  from  the  experiments  that 


the  normalized  shear  strength  (with  respect  to  0.01%/min)  increases  linearly  with  the 
logarithm  of  strain  rate.  The  increase  in  undrained  shear  strength  in  the  pressuremeter 
stress  path  is  about  14.3%  per  log  cycle  for  the  kaolin  clay  and  15.3%  for  the  kaolin- 
silica  mixture.  The  undrained  shear  strength  in  the  conventional  triaxial  test  was  found 
to  increase  about  8  to  10%  for  a  tenfold  increase  in  strain  rate.  Therefore,  it  can  be 
concluded  that  the  undrained  shear  strength  increases  about  40-50%  more  in 
pressuremeter  stress  path  tests  than  in  the  triaxial  stress  path. 

From  the  numerical  studies,  the  influence  of  consolidation  and  creep  were  studied. 
Consolidation  tends  to  increase  the  strength  around  the  probe  for  slower  strain  rate  cases, 
however,  creep  tends  to  decrease  the  strength.  The  creep  effect  is  more  than  the 
consolidation  influence,  and  the  combined  effect  of  consolidation  and  creep  gives  lesser 
shear  strength  at  smaller  strain  rates  and  higher  strength  at  higher  strain  rates. 

4.  The  model  developed  based  on  cavity  expansion  theory  is  able  to  incorporate 
the  influence  of  decreasing  strain  rate  along  the  surrounding  soil  mass.  The  findings 
from  the  strain  rate  test  for  the  pressuremeter  was  used  in  the  model  to  estimate  the 
difference  in  the  interpretation  of  the  undrained  shear  strength.  The  influence  of  the 
other  parameters  also  studied  parametrically.  The  level  of  upper  yield  did  not  affect  the 
strength  prediction.  The  slope  of  normalized  shear  strength  versus  logarithmic  strain  rate 
line,  0  was  found  as  0. 15  and  0. 10  for  pressuremeter  tests  and  triaxial  tests,  respectively. 
For  the  higher  values  of  0  (PMT  condition),  the  peak  strength  shows  about  a  25% 
decrease,  and  the  strain  softening  is  also  higher  than  for  the  triaxial  case.  In  fact,  for 
0  =  0.10,  the  actual  material  curve  and  the  PMT  interpreted  curves  are  of  strain 
hardening  type,  and  the  0  =0.15  curves  are  strain  softening.  Thus,  the  different  test 
condition  influences  the  type  of  results  obtained.  The  important  conclusions  from  this 
study  are  that  the  PMT  condition  induces  more  strain  softening  than  the  triaxial  test 
condition,  which  agrees  well  with  previous  studies,  and  higher  strain  rate  tests  (PMT) 
show  higher  shear  strength. 


S.  The  stress  relief  causes  unloading  of  the  soil  around  the  borehole  and  leads 
to  underestimation  of  the  in  situ  horizontal  pressure  which  in  turn  influences  the 
interpreted  shear  strength  from  the  pressuremeter  expansion  curve.  It  was  found  from 
the  full  scale  tests  that  1.1  %  oversized  cutting  shoe  tests  led  to  an  overestimation  of  60- 
100%  in  undrained  shear  strength. 

The  relaxation  time  or  normalization  period  is  the  time  period  which  should  be 
allowed  to  dissipate  the  excess  pore  pressure  developed  during  the  drilling  operation. 
This  period  is  dependent  upon  the  drilling  technique  used,  the  operator  and  the  type  of 
soil  tested.  In  the  literature,  the  relaxation  time  reported  varies  from  a  few  minutes  to 
several  hours.  The  pressuremeter  test  was  simulated  in  CSD  and  the  specimen  was 
strained  to  different  level  of  strains  using  different  strain  rates  and  the  relaxation  time 
was  monitored.  Irrespective  of  the  initial  strain  level,  for  kaolin  clay  the  relaxation  time 
was  about  200  to  250  minutes  in  all  cases. 

Ladanyi  (1982)  developed  equations  to  obtain  creep  parameters  from  the 
pressuremeter  test  in  a  borehole.  Lacerda  (1976)  proposed  several  correlations  to  obtain 
relaxation  parameters  from  Singh-Mitchell  creep  equation  and  empirically  from  plasticity 
index.  Using  the  aging  theory  of  creep,  Ladanyi  et  al.  (1991)  derived  expressions  for 
stress  relaxation  from  which  creep  parameters  could  be  inferred.  These  creep/relaxation 
parameters  are  required  to  predict  the  time  dependent  behavior  of  soil. 
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Fig.  B.l  Shear  Stress  Vs  Radial  Strain  for  Kaolin  Clay 
(Test  No.  35,  Radial  Strain  Rate  =  0.01%/min) 
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Fig.  B.2  Shear  Stress  Vs  Radial  Strain  for  Kaolin  Clay 
(Test  No.  47,  Radial  Strain  Rate  =  0.01%/min) 
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Test  38 


Shear  Stress  Vs  Radial  Strain  for  Kaolin  Clay 
(Test  No.  38,  Radial  Strain  Rate  =  l.(X)%/niin) 


Test  39 


Shear  Stress  Vs  Radial  Strain  for  Kaolin  Clay 
(Test  No,  39,  Radial  Strain  Rate  =  1.00%/min) 


Shear  Stress  Vs  Radial  Strain  for  Kaolin  Clay 
(Test  No.  40,  Radial  Strain  Rate  =  5.00%/min) 
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Fig.  B.20  Shear  Stress  Vs  Radial  Strain  for  Kaolin  -  Silica  Mix 
(Test  No.  56,  Radial  Strain  Rate  =  0.01%/min) 
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Shear  Stress  Vs  Radial  Strain  for  Kaolin  -  Silica  Mix 
(Test  No.  77,  Radial  Strain  Rate  =  0.05%/min) 


Shear  Stress  Vs  Radial  Strain  for  Kaolin  -  Silica  Mix 
(Test  No.  60,  Radial  Strain  Rate  =  0.10%/ntin) 


Shear  Stress  Vs  Radial  Strain  for  Kaolin  -  Silica  Mix 
(Test  No.  70,  Radial  Strain  Rate  =  0.10%/min) 


Test  7 


iolin  -  Silica  Mix 
0.50%/min) 


Shear  Stress  Vs  Radial  Strain  for  Kaolin  -  Silica  Mix 
(Test  No.  67,  Radial  Strain  Rate  =  0.50%/min) 
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STRESS  RELAXATION  PLOTS  FOR  PMT 


Relaxation  Curves  After  0.1%  Strain  (Undrained 


Relaxation  Curves  After  0.5%  Strain  (Undrained 


Relaxation  Curves  After  0.5%  Strain  (Undrained 


Stress  Relaxation  Curves  After  1.0%  Strain  (Undrained  -  Test  87) 


Relaxation  Curves  After  1.0%  Strain  (Undrained 


Relaxation  Curves  After  0.1%  Strain  (Drained  -  Test  99) 
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APPENDIX  D 


Characteristics  of  Kaolinite 

Trade  Name 
Chemical  Name 
Chemical  Family 
Molecular  Formula 

Color 

Fineness  (thru  325  Mesh) 

Free  Moisture 
Particle  Size 

Less  than  2  Microns 
Greater  than  5  Microns 
Alumina  -  % 

Ignition  Loss  -  5  % 

Manganese  in  deleterious  form 
Silica  -  % 

GE  Brightness 

Producer 

Address 

Phone 

Fax 


AKROCHEM  SC-25  Clay 
Kaolin  or  Kaolinite 
Aluminum  Silicate 
AljOj,  2Si02,  2H2O 

Light  Cream 
99.7% 

1%  Max. 

61% 

20% 

38.5- 39.5 

13.6- 14.0 
Nil 

43.2-44.5 

79 

Akrochem  Corporation 

255  Fountain  Street 

Akron,  Ohio  44304 

(216)  535-2108  (800)  321-2260 

(216)  535-8947 
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• 

Characteristics  of  Ground  Silica 

Trade  Name 

sil-co-su  mo 

Chemical  Name 

Silica 

Chemical  Family 

Silicate 

Molecular  Formula 

SiOj 

Mine 

Ottawa,  IL 

Sieve  Analysis  (ASTM  C  371-56) 

Cum.  %  +  100  Mesh 

11.0 

Cum.  %  +  200  Mesh 

35.4 

Cum.  %  +  325  Mesh 

50.1 

Particle  Size 

Median  (/x)  (ASTM  C-958) 

45.0 

Average  (/x)  (ASTM  B-330) 

10.2 

Specific  Surface  Area  (cmVg) 

2226 

(ASTM  B-330) 

Oil  Absorption  (lbs/ 100  lbs) 

15.0 

(ASTM  D-1483) 

Moisture  (Max) 

0.5% 

Bulking  value  gal/lb 

0.4536 

Ibs/gal 

22.05 

Refractive  Index 

1.54 

GE  Brightness 

76.3 

Yellowness 

0.040 

pH  Value,  37.5%  solids 

8.1 

Hardness  (Moh’s  scale) 

7 

Chemical  Analysis  (%) 

(ASTM  C- 146-72) 

• 

SiOz  99.8  TiOz 

0.013  AI2O3  0.047 

FljOj  0.015  CaO 

<0.01  MgO  <0.01 

0.09 
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10 . * . 

20 

30  ‘  *  STRESS.BAS  -  To  Road  and  Raooid  Data  Iftaoturaa  and  Oioplacamanttl  And 

40'*  To  Pacfofm  1-0  Conaoddatton 

46  '• 

60 . 

66  ‘OCCIARE  SUB  DASH8 
00  CtSiKEY  OFF 

as  DIM  DISPmmi0t.DISP|10l,CAUVDT(10|,TEMF«|11(M.PmmAL(6|.P(6l 
aa  DIM  CALTRANS(6I.  TRANMdl,  0I0%(1I.  LT»(1| 

70  KEYOl  ON:  ON  KEY(1|  OOSUB  1660 
80  KEYI2I  ON:  ON  ICEY(2I  GOSUB  1690 
90  KEYI3I  ON:  ON  KEY(3I  OOSUB  laSO 
100  ICEY(4|  ON:ON  KEY(4|  OOSUB  2260 
110  ICEYI6I  ON:ON  ICEY(6>  OOSUB  2360 
120  KEYIB)  ON:ON  KEYIB)  OOSUB  2380 
130  KEYI7)  ON:ON  KEYI7)  OOSUB  3180 
140  OPEN  *b:data1'  AS  #2  LEN  -  60 

160  FIELD  #2,  4  AS  XF*.  4  AS  YF«,  4  AS  ZF«,  4  AS  BF*.  4  AS  X1F».  4  AS  X2FI.  4  AS  Y1  Ft,  4  AS  Y2Ft.  4  AS  ZIFt,  4  AS  Z2Ft,4 
AS  Y21Ft,  4  AS  Y22Ft.8  AS  TIMEFt 
160  M  -  LOF|2)/6a 

1 80  REM  A  routina  for  loadina  daih8.bin  outiida  basic  woric  space. 

1 90  REM  May  bo  margod  at  the  beginning  of  a  program. 

200  'DEF  SEO  >  SH4000  ‘Change  this  toad  addraas  to  atdt  your  memory 
210  REM  Loads  at  92k.  A  zero  added  automatically  at  right  of  Sh1700  making  it 
220  REM  Shi  7000. 

230  'BLOAD  ‘dashB.bin'.O 

240  OPEN  ’dashS.adr*  FOR  INPUT  AS  #1 

250  INPUT  f1.  BASAOR% 

260  CLOSE  «1 
270  DASH8  -  0 
280  FLAQ%  >  0 
290  M0%  ■  0 

300  CALL  DASH8  (MOIL.  BASAOR%.  FLAO«) 

310  IF  FLAO%  <  >0  THEN  PRINT  'tostolation  arror* 

320  M0%o10:  DiU%(0|o2:  DIO%l1|o2  'Rata  gonorotor  vvith  counter  2 
330  CALL  0ASH8  (M0«.  DIO%(OI.  FLAO%| 

340  M0%«11:  0l0%|1|ir2000 

350  CALL  DASH8  (MOIL.  DIO%IOI.  FLAO%l 

360  ' 

370  CALLVOT(OI»41.91:CALLVOT11|«41.56:CALLVOTi2|-41.24:CALLVDTI3l-41.78: 

CALLVDT|4|>4O.91:CALLVOT(S|>41.13:CALLVOTiai>41.06:CALLVDTI7l«4O.87  'Calibration  constants.  Volts/inch 

380  CALTRANSIOI  -  2.027:  CALTRANSd )  -  2.02:  CALTRANSI2I > 2.028:  CALTRANSI3I  >  2.014  'Calibration  Constsntsin  PSI/m.Volt 

390  ' 

400  ' 

410  LOCATE  1,1:  PRINT  ‘miTTAL  READINGS  OF  LVDTs:*:  PRINT 

420  LOCATE  13.1:  PRINT  'INITIAL  READINGS  OF  PRESSURE  TRANSDUCERS:’: PRINT 

430  LOCATE  15.1:  PRINT  'Ch.  No.':TA8dOI:'lnitial  prassursipa'I* 

440  LOCATE  22,1:PRINT  'Press  s  to  record  data' 

450  LOCATE  23.1: PRINT  'Press  F4  to  check  LVOTs/Transduears' 

460  LOCATE  3.1;  PRINT  'Ch.  No.';TABdOt;'Volts' 

470  ' 

480  'INITIAL  KAOINCS  OF  THE  LVDTS 
490  FOR  1%  -  0  TO  7 
500  OOSUB  1360 

610  DlSPmrridtl  -DISPIdt)  'in  Inchos 

620  LOCATE  l«  4  4.1 :  PRINT  USINO  '  ##  «««.F««':l«:VOLT 

630  NEXT  l« 

640  PRMT:  PRINT 
650  ' 

660  'mrriAL  readbios  of  the  transducers 

670  PINITIAL(0|-3.6:  PmiTIALdl- 1.2:PMmALI2l ■•.3:PB«niAll3) -  II 

680  FOR  l«  •  0  TO  3 

690  GOSUB  1170 

600  PIIHI  «  P1I%I  -  PINITIALIdtl 

610  LOCATE  1 6  ♦  1%.  1 :  PRINT  USMO  *  ##  ««# .«f «':I%:P1I%| 

620  NEXT  1% 

630  ' 
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640  ' 

650  IF  INKEYI  «  ’•*  GOTO  670 

660  GOTO  490 

670  N  -  0;  L-6 

680  GOSUB  720 

690  GOSUB  850 

700  END 

710' 

720  CIS 

730  LOCATE  1.1:PRMT  ’DISPLACEMENTS  M  MCHES:* 

740  LOCATE  3.1:PB»IT  TAB16I;  ’OXI’;  TAB(13I:  •DX2’:  TAB(21>:  ’DYI’;  TAB(29I:  •DY2’;  TAB(37I;  ’DZI’:  TAB(45);  ’0Z2*; 
TABI53):  *DY21’:  TABIBII;  ’DY22’ 

760  LOCATE  9.1;PRMT  'PRESSUflES  M  PSI:' 

760  LOCATE  11.1:PfUNT  ’SigiiiaX*:TAB(12l:‘SigmaY*:TAB(23l:'SiainaZ‘;TAB(31l:’PBra  Ptanura’ 

770  LOCATE  17.1:PRINT  ’Numbar  of  DaU  Raeordad:  * 

780  LOCATE  19.1;PRINT  *Pnta  F2  to  9uk' 

790  LOCATE  20.1  :  PRINT  'ftaaa  F3  to  ehanga  toeoidino  rata* 

800  LOCATE  21,1:  PRINT  ’^aa  F5  to  chock  maHunetioning  of  LVDTa/Trancducara’ 

810  LOCATE  22,1:PRMT  *taaa  F6  to  bogin  1-0  conaoSdatkm  wdth  aaivo  control* 

820  LOCATE  23.1:  PRINT  USING  '1  in  ««#  raadbiga  am  being  facordad.’:L 
830  RETURN 
840  ' 

850  FOR  IH-OTO  7 
860  GOSUB  1360 

870  0ISP|l%l  >  0ISP|l%|  •  OISPINmiTi) 

880  NEXT  1% 

890  FOR  1%  «  0  TO  3 

900  GOSUB  1170 

910  Pll%l  •  PIIKI  -  PINITIAL(IK>| 

920  NEXT  1% 

930  LOCATE  5,1:  PRINT  USING  '  #d.«f«f':DISP(0):0ISP|1):DISP|2):DISP|3);0ISP|4):DISP|5):DISPI6l:0ISP(7l 
940  LOCATE  13.1;PRINT  USING  '  ###.#  •;P(01;P(1);P|2I;P(3» 

950  N  >  N  -f  1 

960  IF  L  >  1  GOTO  980 

970  IF  I  N  MOO  L)  <>  1  GOTO  850 

980  LSET  XFt  -  MKStlPlOH 

990  LSETYFt  >  MKStIPIIII 

1000  LSETZFt  >  MKS«IP(2II 

1010  LSET  BF«  -  MKS«(P(3II 

1020  LSETX1F4  -  MKSKOISPIOII 

1030  LSETX2F«  -  MKStlDISPim 

1040  LSETY1F*  «  MKS«I0ISPI2II 

1050  LSETY2F«  •  MKS«IDISP(3|| 

1060  LSETZ1F4  a  MKS«(0ISP|4)| 

1070  LSET  Z2Ft  >  MCStlOISPISII 
1080  LSET  Y21Ft  >  MKS«(0ISPI6)| 

1090  LSET  Y22F«  >  MKStlDISP(7|) 

1100  LSETTIMEFt  a  TIME* 

1110  M  a  M  -t-  1 

1120  PUT  f2.M 

1130  LOCATE  17,28:PRINT  M 

1140  GOTO  860 

1150  RETURN 

1160  ' 

1170  'SUBROUTINE  PRESSURE 

1180  LT%|0|a1:  LT«|1|a1:  M0%a1  'Locka  0ASH8  on  channel  1 
1190  CALL  0ASH8  |M0%,  LT»(OI.  FLA0%i 

1200  IF  FLA0«<>0  THEN  PRINT  ’Error  In  aatting  tha  charmal':  END 
1210  M>«al4 

1 220  CALL  0ASH8  (MD«,  IT*.  FLA0%| 

1230  MD%aS 

1240  TRAN%IO|aVARPTR(TEMP%|OII 
1250  TRAN«I1|  a  100 

1 260  CALL  DASH8  (M0«,  TRAN«IO).  FLAG**) 

1270  XaO  :  INDEX  a  0 
1280  FOR  J«1  TO  98 

1290  IF  ABS|TEMP%IJ)-TEMP«IJ-1||  <2  AND  ABSrrEMP«IJ)-TEMP%(J -f  1|)<2  THEN  X-X  +  TEMP%(J|:  INDEX  -  INDEX  +1 
1300  NEXT  J 
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1310  VOLT  -  Xy(M0EX*20.4n  'bi  m.Volt 
1320  P(l»)  <■  VOLT*CAL'mANS(l%)  ‘In  pal 
1330  RETURN 
1340  ’ 

1360  ' 

1300  'suenouTWE  displacement 

1370  LT%(0|>i0;  LT%(1l>0:  MO%-1  ‘Locka  OASHO  on  ehannal  0 
1380  CAU  DASH8  IMO«.  LTKIOI.  FLAO%l 

1380  F  FLAG«  <>0  THEN  FRMT  *Eiror  In  aottina  tho  etiannai’;  END 
1400  MO%-14 

1410  CALL  DASH8  |M0«.  1%.  FLAO«l 
1420  MOK«6 

1430  TRANIilOt  -  VARPTROEMPlilOti 

1440TRANH(1|  -  100 

1450  CAa  0ASH8  IMO«,  TRAN«|0|,  FLAQ«| 

1400  X  -  0:  INDEX  -  0 
1470  FOR  J-1  TO  08 

1480FABS(TEMP»|J)-TEMP»|J-1»  <2  ANOABSrrEMP«g) -TEMP«|J-f  1))  <2  THEN  X-X-f  TEMP«|J);  INDEX  -  INDEX  -fl 
1490  NEXT  J 

1500  VOLT  -  X/(204.8*M0EXI  'In  VoKa 
1610  F  ABS(VOLT)  >  16  GOTO  1370 
1520  OISP(IH)  «  VOLT/CALLVOT|l«|  ‘In  Inchea 
1530  RETURN 
1540  ‘ 

1560  ‘SUBROUTINE  TO  FRESHEN  THE  SCREEN  Fit 
1500  GOSUB  720 
1570  RETURN 
15B0  ‘ 

1590  ‘SUBROUTINE  TO  TERMINATE  THE  RUN  (KEY  F2| 

1600  CLOSE  #2 
1610  OUT  SH320,0 
1620  CLS 
1630  END 
1640  ' 

1650  -SUBROUTINE  TO  CHANGE  RECORDING  RATE  (KEY  F3) 

1660  CLS 

1670  INPUT  "What  ia  tha  recording  rata  (a.g.  1  In  LI  L':L 
1680  GOSUB  1550 
1690  RETURN 
1700  ‘ 

1710  'SUBROUTINE  TO  CHECK  MALFUNCTIONING  OF  LVDTa  and  TRANSDUCERS 
1720  CLS 

1730  INPUT  ’Do  you  want  to  ctwck  LVDTa  (yM)':AX« 

1740  F  LEFT«(AX«,1|  -  ’n’  GOTO  1900 
1750  F  LEFT$(AXt,1)  »  'v*  GOTO  1770 
1760  PRINT.PRINT  'Say  either  y  or  n  only':  GOTO  1730 
1770  LOCATE  10,1 
1790  PRINT  •• 

1790  PRINT  *• 

1B00  PRINT  •• 

1810  PRINT  ’• 

1820  PRINT  ’• 

1830  PRINT  •• 

1840  PRINT 
1850  PRINT  *• 

1855  PRINT  *• 

1860  PRINT  ’• 

1870  F  MKEY*-"  THEN  1870 
1880  CLS 

1890  INPUT  'Which  channel  (0  through  71  do  you  want  to  chock  ';l% 

1900  GOSUB  1360 

1910  PRINT  USING  '«#  ##.ff';l«:VOLT 

1920  F  mKEYto"  GOTO  1900 

1930  INPUT  'Do  you  want  to  chock  any  othar  channola  ly/nl  ';XXt 
1940  F  lEFTIIXXt.D  -  'n'  GOTO  1970 
1950  F  LEFTtlXXt.ll  >  'y'  GOTO  1880 
1960  PRINT:PRINT  'Say  oithor  y  or  n  Of.ly':  GOTO  1930 
1970  • 


CHECKING  THE  LVDTa 

Tho  uottago  biducod  bi  tho  LVDTa  and  coireaponding  channel 
numbora  will  bo  dieplayod.  A  ataady  fkiw  of  aoomingly  conatant 
nunnbara  (with  laaa  than  e/-0.02  volt)  bidicataa  that  tho  LVDT 
ia  functionbig  aotiafaetorly.  LVDTa  are  chackod  ona  at  a  tima. 
PRESS  ANY  KEY  WHEN  READY.  Good  Luekl 


1980  CLS 

1990  INPUT  'Oo  you  want  to  ohaek  prataura  trwiaducan  (y/nl':8X* 

2000  F  LEmiBXt.ll  -  'n*  GOTO  2230 
2010  F  L£FT*(BX«.1|  -  V  GOTO  2030 
2020  PRINT:PRINT  'Say  aWiar  y  or  n  only':  GOTO  1990 
2030  LOCATE  10.1 

2040  PRINT . 

2060  PRINT  '*  CHECKMG  THE  PRESSURE  TRANSDUCERS 

2060  PRINT  " 

2070  PRINT  '  *  Tha  pcoaaufo  and  eonaapondina  channal  numbart  wlH  ba  diaplayad. 

2080  PRINT  '*  A  ataady  flow  of  aaamlnoly  oonatant  vahiaa  of  prataura  {with 
2090  PRINT  '*  lata  ttian  ■¥!•  0.6  pal  varfatlonl  Mieataa  ttwt  tha  prataura 
2100  PRINT  "  trarMduear  It  functioning  aatlafaetorlly.  Thty  ara  balrxi 
2110  PRINT"  chackad  ona  at  a  tbiM.  PRESS  ANY  KEY  WHEN  READY.  Good  LuckI 
2116  PRINT" 

2120  PRINT . 

2130  F  INKEYt  -  "  THEN  2130 
2140  CLS 

2160  INPUT  'Which  channal  do  you  want  to  chack  ';■« 

2160  GOSUB  1170 

2170  PRINT  USING  '««  ««f.f':l%:P(l14| 

2180  F  INKEYt  -  "  GOTO  2160 

2190  INnjT  'Do  you  want  to  chack  other  chaniMla  (y/n|  ':YYt 
2200  F  YYt  -  'n'  GOTO  2230 
2210  F  YYt  -  'y'  GOTO  2140 

2220  PRINT:PRINT  'Say  althar  y  or  n  only':00T0  2190 
2230  RETURN 
2240  ' 

2260  GOSUB  1710 
2260  CLS 

2270  LOCATE  1,1;PRINT  'miTUL  READINGS  OF  THE  LVOTS:';  PRINT 

2280  LOCATE  13,1;PRINT  'INITIAL  READINGS  OF  THE  PRESSURE  TRANSDUCERS:';  PRINT 

2290  LOCATE  16,1;PRINT  'Ch.  No.':TAB(10);'lnliial  Prataura  (pail' 

2300  LOCATE  22,1;PRINT  'Plata  a  to  laconf  data' 

2310  LOCATE  23,1:PRINT  'Piatt  F4  to  chack  LVOTt/Trantducara' 

2320  LOCATE  3.1:PRINT  'Ch.  No.';TAB(10|:'VoKa' 

2330  RETURN 
2340  ' 

2350  GOSUB  1710 
2360  GOSUB  720 
2370  RETURN 
2380  ' 

2390  CLS 

2400  LOCATE  9,16:fWNT . 

2410  LOCATE  10.16:PRINT  " 

2420  LOCATE  11,1B:PRINT  " 

2430  LOCATE  12,16:PRINT  "  SERVO  CONTROLLED  ONE  DIMENSIONAL  CONSOLIDATION 
2440  LOCATE  13,1S:PRINT  " 

2450  LOCATE  14,15:PRINT  " 

2460  LOCATE  . . 

2461  KEYI3I  OFF 

2462  KEYI3)  ON:  ON  KEYI3I  GOSUB  4000 
2470  FOR  I  -  1  TO  60;  NEXT  l:BEEP;BEEP 

2480  LOCATE  9,16:PRWT . 

2490  LOCATE  10.16;PRINT  "  Watching  tha  aeraan  for  dlaplay,  adjuat  aigmaX. 

2600  LOCATE  11,15:PRINT  "  dgmaY,  and  aigmaZ  auchlhat  lhay  aia  equal  to  tha 

2510  LOCATE  12.16:PRINT  "  laqubad  TOTAL  vaitical  eenaoMation  piataun. 

2620  LOCATE  13.16:PRINT  "  aigmaX  and  tigmaY  ara  daciaatad  AUTOMATICALLY 

2630  LOCATE  14,16:PRMT  "  during  1-0  conaoBdatlon,  makrtalning  zoio  lateral  atrain. 

2660  LOCATE  16,16;PR»fr  " 

2660  LOCATE  17.16;PRINT  "  NOTE:  Al  diaplaeanionta  ahouW  ba  about  0.0000.  If 

2670  LOCATE  t8,15:PRINT  "  dlffaiant.  NEWITIALIZE  IPRESS  F2). 

2680  LOCATE  19,16:PRINT  " 

2690  LOCATE  20,16:PRINT  "  FOR  MANUAL  OVERRDE  DURING  CONSOLIDATION  PRESS  F7 

2600  LOCATE  21.16:PRINT  " 

2610  LOCATE  22.15:PRINT  "  Whan  laady  open  the  drainaga  valut  and  prata  any  key 

2620  LOCATE  . . 

2630  LOCATE  1,1;PRINT  TAB(6|;  'DXI';  TAB(13I;  '0X2*;  TABI21I;  *DY1';  TAB(29I;  'DYZ';  TAB(37|;  ‘OZI";  TAB(45I;  "OZZ"; 
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TABIS3I:  •0Y21’:  TABI61):  *DY22’ 

2640  LOCATE  4,1;PniNT  'SiamaX':TAB|12|:‘Sioin4Y’:TAB(21|:'SlgmaZ*:TAB(30):*Pore  Praiiure* 

2660  FOR  IK-OTO  7 
2660  OOSUB  1360 

2670  OISP(l%)  >  DISP(I%I  -  OISPINmiK) 

2680  NEXT  1% 

2680  FOR  IK-  0  TO  3 

2700  OOSUB  1170 

2710  PIIKI  -  PIIKI  -  PMITIALIIKI 

2720  NEXT  IK 

2730  LOCATE  2.1:  FRMT  USING  '  f«.«*«#';OISP(OI:OISP(1):OISPI2l;DISP(3):DISP|4|;OISP(5):DISPI6):DISP(7| 

2740  LOCATE  6,1:PRINT  USINO  '  «##.«  *:P(0I;P(1|:P(2|:P(3I 

2760  F  MKEYt  -  "  GOTO  2660 
2760  OOSUB  720 

2770  LOCATE  22,1:PRMT  'Pratt  F7  for  manual  owanida  during  1-D  contoldatlon’ 

2780  OUT  BH323.  AH80 
2780  OUT  BH320.6 
2800  DELOISP  -  .0006 
2810  FOR  IK  -  0  TO  7 
2820  GOSUB  1360 

2830  OISPIIKI  -  DISPdKI  ■  OISPINITIIKI 
2840  NEXT  IK 

2850  DISPX  -  OtSPIOI  -t-  0ISP|1| 

2860  DISPY  «  0ISP(2|  *  DISPI3I 

2870  IF  DISPX  -  DELDISP  >0  THEN  OUT  &H320,7:F0R  J-1  TO  20:NEXT  J:OUT  &H320.5 

2880  F  OISPX  -f  DELDISP  <  0  THEN  OUT  &H320.4:FOR  J-1  TO  120;NEXT  J:OUT  &H320.5 

2890  IF  OISPY  -  DELDISP  >0  THEN  OUT  &H320.13;FOR  J-1  TO  20:NEXT  J:OUT  &H320.5 

2900  F  OISPY  +  DELDISP  <0  THEN  OUT  &H320.1:FOR  J-1  TO  120:NEXT  J:OUT  &H320.5 

2910  FOR  IK  -  0  TO  3 

2920  GOSUB  1170 

2930  PIIKI  -  PIIKI  •  PINITIALIIKI 

2940  NEXT  IK 

2950  LOCATE  5.1:  PRINT  USING  *  ««.«««#' :0ISP(0):DISP(1I:DISP|2):DISP(3);DISP|4|;DISP|5):DISPI6):0ISP(7| 

2960  LOCATE  13.1:PRINT  USING  *  «#f.f  ■•.P|0>;P|1);P|2):PI3) 

2970  N  «  N  4  1 

2980  F  L  >  1  GOTO  3000 

2990  F  (  N  MOD  LI  <>1  GOTO  2810 

3000  LSET  XF*  -  MKStIPIOII 

3010  LSETYFt  -  MKStlPIHI 

3020  LSET  ZFI  -  MKS*(P(2II 

3030  LSET  BFf  -  MKS«(P|3)I 

3040  LSET  X1F*  -  MKS«IOISPIO)| 

3060  LSET  X2F«  -  MKStlOISPIIH 
3060  LSET  Y1F«  -  MKS«I0ISP|2II 
3070  LSET  Y2F«  -  MKSII0ISPI3II 
3080  LSETZIFt  -  MKS$(DISP|4|| 

3090  LSET  Z2F*  -  MKS«(OISP|SI| 

3100  LSETYZIFf  -  MKStlOISP(6ll 
3110  LSETY22Ft  -  MKS1I0ISP(7I| 

3120  LSETTIMEFt  -  TIME! 

3130  M  -  M  1 

3140  PUT  «2.M 

3150  LOCATE  17.28:PRINT  M 

3160  GOTO  2810 

3170  RETURN 

3160  ‘ 

3190  CLS 

3200  OUT  AH320.6 

3210  LOCATE  10.6:PRINT  -|1|.  Oota  ttw  plug  valva  IgraanI  lor  X.  Y  diractlont' 

3220  LOCATE  12,6:PRINT  *|2|.  Adjuat  tigmaX,  tIgmaY  ragulatora  to  tho  raadingt  ditplayad  on  tcraan* 

3230  LOCATE  14.5;PRINT  ‘Pratt  tny  kay  triwn  fMtliad* 

3240  LOCATE  1,1;PRINT  TABISI;  'OXI';  TAB(13};  '0X2':  TABI21I;  'DYI';  TABI28I;  *DY2':  TABI37I:  'DZI*;  TAB(45);  ■022-; 
TABI53I;  ■DY21’;  TABI61I;  ■DY22' 

3260  LOCATE  4.1:PRINT  ‘SigmaX*:TAB(12l;‘SigmaY*:TAB(21):'SigmaZ*:TAB(30l;'Pora  Prattura' 

3260  FOR  IK-0  TO  7 
3270  OOSUB  1360 

3280  OISPIIKI  -  OISPIIKI  -  OISPINITIIKI 
3290  NEXT  IK 
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3300  FOR  l«-  0TO3 

3310  QOSUB  1170 

3320  P(l%)  -  P(l%>  -  PMITULIIK) 

3330  NEXT  1% 

3340  LOCATE  2.1:  PRINT  USINO  '  «#.##««*:OISP(OI:OISI>(1);DISP|2):OISP(3);DISP|4);OISP|5I:DISP(6):DISP|7) 
3360  LOCATE  6.1:PRINT  USING  '  «««.«  ■:P(0I;P|1):P|2I:P(3I 

3360  F  MKEY«  -  **  GOTO  3260 
3370  CLS 

3360  OUT  &H320,0 

3390  LOCATE  10,e:PRINT  *Opwi  tha  plug  valvaa  for  x.  y  dirM^tiont' 

3400  LOCATE  14,6:PRINT  *tata  any  kaya  whan  tha  valwaa  ara  opanad' 

3410  F  MKEYI  -  '*  THEN  3410 

3420  GOSUB  720 

3430  LOCATE  22.1:PRINT  SPC(60I 

3440  GOSUB  860 

3460  RETURN 

4000  ' 

4010  GOSUB  1660 

4020  LOCATE  22,1:PRMT  'Praaa  F7  for  Manual  Ovarrido  during  1-0  ConaoKdation* 

4030  RETURN 
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Iff . 

20 ’• 

30'*  PM.BAS  -  To  Raad  and  Racord  OaU  (ftaMuni  and  DliplacaimnttI 

40  '*  during  Strain  ControRod  Loading 

50  •• 

6ff . 

70  CLS:ICEY  OFF 

72  DIM  DISPINrnt0I.OiSP(10),CAllVDT(10|,'T5MP«(n0t.l>MrnALI5I.P(S) 

73  DIM  CALTRANSiS).  TRAN%(1|.  DIO%(1|.  LT%(1) 

80  KEY(1)  ON:  ON  KEYI1I  QOSUB  1710 

90  ICEYI2)  ON:  ON  KEY(2)  QOSUB  1980 
100  KEY(3i  ON:  ON  KEYI3)  QOSUB  1B60 
110  OPEN  'b:daMl’  AS  *2  LEN  ■  SB 

120  FIELD  «2.  4  AS  XF4,  4  AS  YFI,  4  AS  ZF«.  4  AS  BF*.  4  AS  X1F«.  4  AS  X2F«,  4  AS  YIFt,  4  AS  Y2Ft,  4  AS  ZIFt.  4  AS  Z2Ft.4 
AS  Y21Ft.  4  AS  Y22F4.S  AS  TIMEF* 

130  M  -  LOFI2I/66 

150  REM  A  lautkia  for  loading  daahS.bln  outaida  batie  work  apaca. 

160  REM  May  ba  margod  at  tha  baginning  of  a  program. 

170  'DEF  SEQ  ■  BH1700  ‘Changa  tMa  load  addrait  to  auit  your  mamory 
180  REM  Loada  at  92k.  A  taro  addod  automatically  at  right  of  Bh1700  making  it 
190  REM  Bh17000. 

200  'BLOAD  *dath8.bln*.0 

210  OPEN  'dathB.adr*  FOR  INPUT  AS  #1 

220  INPUT  #1.  BASADR« 

230  CLOSE  <1 
240  0ASH8  -  0 
250  FLAQK  •  0 
260  M0%  «  0 

270  CALL  0ASH8  (MD%,  BASADRIL,  FLAQK) 

280  IF  FLAa%  <  >0  THEN  PRINT  'inatallation  arror* 

290  MD%>i10:  DIO«(OI«2;  0I0%(1|>2  'Rata  ganarator  with  countar  2 

300  CALL  DASHB  (M0%,  OlOltlOt.  FLAQ%I 

310  MD%>11:  0l0%(1|o2000 

320  CALL  0ASH8  |M0%.  0I0%(0|,  FLA0%l 

330  ' 

340  CALLVOTIOI >41.91:  CALLV0T(1|-41.66:  CALLV0T(2) >41.24:  CA'.,.V0T(3)>41.7e;  CALLV0TI4)  >40.91 : 

CALLVOT(S|>41.13:  CALL VDTI6I >41.05;  CALLVDT(7| >40.87  'Calibration  conatantt.  Vorts/inch 

350  CALTRANSIOI  >  2.027:  CALTRANSH I  >  2.02:  CALTRANSI2)  >  2.028;  CALTRANSI3I  >  2  31 4  'Calibration  Conatanta  in  PSI/m. VoK 
360  ' 

370  'INITIAL  READINGS  OF  THE  LVDTS 
380  FOR  1%  >  0  TO  7 
390  QOSUB  1520 

400  DISPINITIITLI  >  0ISP|l%|  'ki  inchaa 
410  NEXT  1% 

420  'INITIAL  READINGS  OF  THE  TRANSDUCERS 
430  PINITIAL(0l-3.6:  PINITIALdI >  1 .2;PINmAL(2l >-.3:PtNrnAL<3l>  1 1 
440  N  >  0:  L>6 
450  ■ 

460  'BEGIN  STRAIN  CONTROLLED  LOADING 
470  CLS 

480  LOCATE  9.1 . . 

490  LOCATE  10.15;PRINT  •* 

500  LOCATE  11,16;PRINT  •• 

510  LOCATE  12,15;PRINT  STRAM  CONTROLLED  LOADING 

520  LOCATE  13,16:PRINT  ■'* 

530  LOCATE  14.15:PRINT  •• 

540  LOCATE  1 5. 1 5:  PRINT . 

550  SOUND  1000,50:  SOUND  2000.50:  SOUND  3000.50 
560  OUT  AH323.  SH80 
570  OUT  SH320.21 
560  CLS 

590  LOCATE  7,10:PRINT  'Incraaaa  SigmaX  to  a  vahia  graatar  thar^  axpactad  final  Ifailuia)  value* 

600  LOCATE  9,10;PRWT  ’Praaa  any  kay  to  begin  attain  comrodad  loading* 

610  IF  MKEY*  >  **  THEN  610 
620  CLS 

630  LOCATE  9,10;INPUT  *What  ara  tha  X,  Y  wMtha  of  tha  apacimon  aftar  conaoMation  |in.|*:HX,HY 
640  LOCATE  1 2.10;INPUT  'What  la  tha  daabad  X-atrain  rata  (%  par  min.)*:  RATEX 
642  LOCATE  13.10:INPUT  *What  la  tha  doabad  Y-atraIn  rata  («  par  min.)  (NEG.I*;  RATEY 
650  LOCATE  14,10;PRMT  'Make  aura  that  tha  ifrain  la  eloaad  bofora  you  procaad* 
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660  LOCATE  16.ia.MI>UT  'Do  you  wMi  to  chango  tho  abovo  information  (y/n|';Xt 

670  V  IEFT*(X»,1J  -  ’y"  GOTO  680 

680  IF  LEFT$(X».1|  ■  "o’  GOTO  700 

680  PRINT:PRINT  *Say  oMiar  Y  or  n  only':  GOTO  660 

700  OELX  -  .01  *RATEX 

701  DELY  -  .OI'ABS(RATEY) 

702  DELOISP  -  0.0006 
706  NNN  >  400 

710  CLS 

720  GOSUB  1900  'Thnlno  Bagina 
726  TIME1  >  SS 
730  GOSUB  1200 
760  FOR  l%>  OTO  6 
760  GOSUB  1620 

770  DISPII%I  -  OISPdMI  •  OISPmiTdKI 
780  NEXT  l« 

782  DISPX  -  OISPIOI  OISPIII 
784  OISPY  -  0ISP(2|  -f  0ISP(3t 
790  DISPZ  ■  DISP(4|  -f  0ISP(5| 

782  STRAINX  >  |DISPX/HX)*100 

794  STRAINY  -  (OISPY/HYl'100 

802  GOSUB  1900 

804  TIME2  -  SS 

810  TIME  «  (TIME2  •  TIME1)/60 

820  EOOTX  >  STRAINX/TIME 

823  EDOTV  >  STRAINYH'IME 

825  IF  TIME  >  16  IREN  NNN  >  200 

830  IF  EOOTX  -  RATEX  >  OELX  THEN  OUT  &H320,23:FOR  J«1  TO  S:NEXT  J:  OUT  &H320.21 

832  IF  RATEX  •  EOOTX  >  OELX  THEN  OUT  &H320,20:FOR  J«  1  TO  NNNiNEXT  J:OUT  &H320.21 

834  IF  EDOTY  -  RATEY  >  DELY  THEN  OUT  AH320,29:FOR  J«1  TO  B.NEXT  J:  OUT  &H320.21 

836  IF  RATEY  -  EDOTY  >  DELY  THEN  OUT  AH320,17:FOR  J>1  TO  NNN:  NEXT  J;  OUT  &H320.21 

842  IF  DISPZ  -DELOISP  >  0  THEN  OUT  8H320,63:FOR  J«1  TO  20:NEXT  J:  OUT  &H320.21 
844  IF  DISPZ  +  DELOISP  <  0  THEN  OUT  8H320.6:FOR  J  .  1  TO  120:  NEXT  J:  OUT  6H320,21 
850  ■ 

860  -FOR  X  -  1  TO  30000:  NEXT  K 
870  FOR  l%»0  TO  7 
880  GOSUB  1520 

890  OISP|l%l  >  DISPIIK)  -  DISPINrT|l%l 
900  NEXT  1% 

910  FOR  0  TO  3 
920  GOSUB  1320 
940  NEXT  1% 

950  LOCATE  5.1:  PRINT  USING  '  ««.«#«f';DISP|0);DISP(1|:DISP|2);DISPI3);DISP|4);DISPI5):DISPI6l:DISP|7) 

960  LOCATE  13,1:PRINT  USING  '  9fi.t  'jPIOliPdllPIZhPISI 

970  LOCATE  21.28:PRINT  USING  'ff.Ffff'iEOOTX 

976  LOCATE  22.28:PRINT  USING  'F/./ftF'iEDOTY 

980  N  >  N  -f  1 

990  IF  L  -  1  GOTO  1010 

1000  IF  (  N  MOD  LI  <>  1  GOTO  750 

1010  LSETXF*  •  MKStlPIOII 

1020  LSETYFt  ■  MKSI(P|1)| 

1030  LSETZF*  ■  MKS«(P(2II 
1040  LSETBFI  >  IMCS«<P|3II 
1050LSETX1F*  -  MKS*(OISP(0)I 
1060  LSET  X2F*  -  MKS4IDISPI1II 
1070  LSETY1F*  >  MKS4IDISPI2)) 

1080  LSET  Y2F«  -  MKS*(0ISPI3II 
1090  LSETZIFt  -  MKS«I0ISP|4|| 

1100  LSETZ2F*  •  MKStlDISPISII 
1110  LSET  Y21FI  -  MKS«(0ISP|6)| 

1120  LSETY22FI  -  MKSt(0ISP|7)| 

1130  LSETTIMEFt  •  TIME* 

1140  M  «  M  4  1 

1150  PUT  #2,M 

1160  LOCATE  17,28:PRINT  M 

1170  GOTO  760 

1 1 80  END 

1190  ■ 

1200  CLS 
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1210  LOCATE  1,1:PRINT  'DISPLACEMENTS  IN  INCHES:* 

1220  LOCATE  3,1:PRINT  TABI6);  '0X1';  TAB(13):  '0X2';  TAB(21):  ’DYI':  TAB(29I:  *DY2';  TABI37I;  "DZI";  TAB(46I;  *DZ2'; 
TABIS3I;  *DY21’;  TABIBII;  'DYZZ' 

1230  LOCATE  9,1;PRINT  'PRESSURES  IN  PSI:* 

1240  LOCATE  11,1:PRINT  *SigmaX*:TAB(12l:'SigmaY*:TABI23);*SiginiZ*;TABI31):*PDra  Prsitura’ 

1250  LOCATE  17,1:PRINT  'Number  of  Data  Raeordad:  * 

1 260  LOCATE  1 9,1  :PRINT  'Praia  F2  to  quK* 

1270  LOCATE  20,1  :  PRINT  'Praia  F3  to  ehanoa  racordkig  rate* 

1275  LOCATE  21,1:PRINT  'Currant  X-Strain  Rata  la  K  par  min.* 

1 260  LOCATE  22,1  :PRINT  'Cunant  Y-Strain  Rata  la  %  par  min.’ 

1290  LOCATE  23,1:  PRINT  USING  '1  in  ###  raadinga  ara  being  racorded.*:L 
1300  RETURN 
1310  ‘ 

1320  'SUBROUTINE  PRESSURE 

1 330  LT%|0)  *  1 :  LT«t1 1  >  1 :  MO%  >  1  'Loeka  DASH8  on  channel  1 
1340  CALL  DASH8  (MOK,  LT%I0I,  FLAQ%I 

1350  IF  FLAG%  <>0  THEN  PRINT  'Error  in  cattbig  tha  channel*:  END 
1360  M0%>14 

1370  CAU  DASH8  (MOK,  l«,  FLAG%I 
1380  MO%«6 

1390  TRAN%(0)«VARFTR(TEMP%(0n 
1400  TRANKd)  -  40 

1410  CALL  DASH8  (MD%,  TRANKIO),  FLAGKI 
1420  X  =  0  :  INDEX  >  0 
1430  FOR  J  =  1  TO  38 

1440  IF  ABS|TEMP%IJI-TEMP%(J-1))  <2  AND  ABS(TEMP%(J)-TEMP%IJ-f  1))<2  THEN  X  =  X -t'TEMP%IJ):  INDEX  =  INDEX  -f  1 
1450  NEXT  J 

1460  VOLT  *  X/IINDEX '20.481  'In  m.Volt 
1470  P(H4I  =  VOLT'CALTRANSIHil  'In  pii 
1480  P(l%)  =P(I%I  -  PIN1TIAL(I%1 
1490  RETURN 
1500  ' 

1510  ' 

1520  'SUBROUTINE  DISPLACEMENT 

1530  LTK(0l«0:  LT%(1)»0:  MD%s1  'Loeka  0ASH8  on  channel  0 
1540  CALL  DASH8  (M0%,  LTHIOI,  FLA6%I 

1550  IF  FLAG%  <  >0  THEN  PRINT  'Error  in  catting  tha  channal*:  END 
1560  MD%  =  14 

1570  CALL  DASH8  (MD%,  1%.  FLAG%| 

1580  MD%i°5 

1590  TRANK>(0|  «  VARPTR(TEMP%(OII 
1600  TRAN%(1)  »  40 

1610  CALL  DASH8  IMD%,  TRAN%(0|.  FLAG%) 

1 620  X  >  0:  INDEX  -  0 
1630  FOR  J»1  TO  38 

1640  IF  ABS|TEMP%(J| -TEMP%IJ-1||  <2  AND  ABS(TEMPH(J)  -  TEMP%IJ  + 1)1  <  2  THEN  X  =  X +  TEMP%(J):  INDEX  =  INDEX  +1 
1650  NEXT  J 

1660  VOLT  -  X/(204.8*IN0EX)  'In  VofCa 
1670  IF  ABSIVOLT)  >  16  GOTO  1630 
1680  DISP(I%I  -  VOLT/CAI  LVOTdlil  'In  inchaa 
1690  RETURN 
1700  ' 

1710  'SUBROUTINE  TO  FRESHEN  THE  SCREEN  IF1I 
1720  GOSUB  1200 
1730  RETURN 
1830  ’ 

1840  ' 

1850  'SUBROUTINE  TO  CHANGE  RECORDING  RATE  (KEY  F3) 

1 860  CLS 

1870  INPUT  'What  ic  tha  recording  rate  (e.g.  1  In  L>  L';L 

1880  GOSUB  1710 

1890  RETURN 

1900  'SUBROUTINE  TIMER 

1910  TIMt  >  TIMEt 

1920  SS«  «  RtGHTt|TIM«,2| 

1930  MMt  «  MIOt(TIMt,4,2) 

1940  HH»  =  LEFT»(TIM».2) 

1950  SS  «  VALISStl'l  +  VAL|MM»I*60  ♦  VAL(HH»r3600 
1960  RETURN 


1970  ' 

1980  'SUBROUTINE  TO  TERMINATE  THE  RUN  (KEY  F2I 
1990  CIS 
2000  CLOSE  *2 

2010  LOCATE  7,10:PRINT  "Watch  the  gaga  and  turn  tha  prasaura  ragulator  to  bring  down  SIOMAX* 
2020  LOCATE  0.10:PRINT  "to  about  SIQMAY  and  SIOMAZ* 

2030  LOCATE  11.10:PRINT  "PRESS  ANY  KEY  WHEN  FWISHEO" 

2040  IF  INKEY*  -""  THEN  2040 
2050  OUT  aiH320.0 
2060  CLS 
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. . 

20  •• 

30  '*  RELAX.BAS  ■  To  Rood  ond  Rocoid  Doto  (ProMuro*  and  Oitplacamenta) 

40  '*  during  Ratoxation  Taat 

60  '• 

. . . . 

70  CLS:KEY  OFF 

72  DIM  DISPINmi0I.DISP(10|.CALLVDT(10).TEMP»ni0i.f>mrnAL|6).P(6| 

73  DIM  CALTRANS(6|,  TRANKI1).  OIOKI1),  LTKOt 
80  KEYdl  ON:  ON  KEY|1|  60SUB  1710 

90  ICEYI2I  ON:  ON  KEY|2)  QOSUB  1980 
100  KEYI3I  ON:  ON  ICEY(3)  QOSUB  1860 
110  OPEN  'b:data1‘  AS  «2  LEN  •  66 

1 20  FIELD  92,  4  AS  XFt,  4  AS  YF«,  4  AS  ZF«.  4  AS  BF*,  4  AS  XIF*.  4  AS  X2F«,  4  AS  YIFt,  4  AS  Y2F*,  4  AS  Z1 F*.  4  AS  Z2F$.4 
AS  Y21F«.  4  AS  Y22F«,B  AS  TIMEFF 
130  M  -  L0F|2)/56 

160  REM  A  loutlna  for  loading  dathS.bln  outaida  baaie  work  tpaca. 

1 60  REM  May  ba  margad  at  tha  baginning  of  a  program. 

1 70  *DEF  SEO  m  BH1 700  ‘Changa  thla  load  addrata  to  auit  your  mamory 
1 80  REM  Loada  at  92k.  A  xaro  addad  automatically  at  right  of  Utl  700  making  It 
190  REM  BhWOOO. 

200  *BLOAD  ’daah8.bln'.0 

210  OPEN  ’daah8.adr'  FOR  INPUT  AS  f1 

220  INPUT  91 .  8ASA0R% 

230  CLOSE  91 
240  0ASH8  <c  0 
250  FLAGKi  >  0 
260  MOIL  °  0 

270  CALL  DASH8  (MOW,  BASADR%,  FLAG«) 

280  IF  FLA0%  <  >0  THEN  PRINT  ‘inatallation  arror' 

290  M0%s10;  OI01L|0|s2:  OIO%(1)c2  ‘Rata  ganarator  with  counter  2 
300  CALL  0ASH8  (MOW,  DIOK(O),  FLA0%) 

310  M0%>11;  Dl0%|1|>2000 

320  CALL  0ASH8  (M0%,  OlOKIOI,  FLAQ%I 

330  • 

340  CALLVDWI  *41.91:  CAttVDTflf  »4f.66;  CALLVDTI21*4I.24:  CAIIVDTJ3)  =41,7&  CALLVDT|4)  =  40.91: 
CALLVDT(5)«41.13:  CALLVDT(6I -41.06:  CALLVDT(7|  >40.87  'Calibration  conatanta.  Volta/inch 

350  CALTRANSIO)  >  2.027:  CALTRANSd )  >  2.02:  CALTRANSI2)  -  2.028:  CALIRANSO)  >  2.01 4  ‘Calibration  Conatanta  in  PSI/m.  Volt 
360  ‘ 

370  ‘INITIAL  READINGS  OF  THE  LVDTS 
380  FOR  1%  «  0  TO  7 
390  GOSUB  1520 

400  DISPINlTdK.)  >0ISP(l%)  ‘In  Inchaa 
410  NEXT  1% 

420  ‘INITIAL  READINGS  OF  THE  TRANSDUCERS 

430  PINITIAL(0I  =  3.6:  PINITIALIH  =  1 .2:PINmAL(21  «-.3:PINrriAL|3|  =  II 

440  N  s  0:  L>S 

450  ‘ 

460  ‘BEGIN  RELAXATION  TESTING 
470  CLS 

480 LOCATE  . . 

490  LOCATE  10.15:PRINT 
600  LOCATE  11,16:PRINT 

610  LOCATE  12,1S:PRINT  RELAXATION  TESTING 

620  LOCATE  13,16:PRINT 

625  LOCATE  14,15:PRINT  (MAINTAVIMQ  NO  DEFORMATION  IN  ALL  DIRECTIONS) 

526  LOCATE  16.16:PRINT 
530  LOCATE  16.16:PRINT 

640LOCATE1 7, 1 6:PRINT . . . 

550  SOUND  1000.50:  SOUND  2000,60:  SOUND  3000,60 
560  OUT  &H323.  &H80 
570  OUT  &H320,21 
580  CLS 

600  LOCATE  7,10:PRINT  ’ftaaa  any  kay  to  bagin  Ralaxation  Taating* 

610  IF  MKEY*  -  "  THEN  610 
620  CLS 

650  LOCATE  14,10:PRINT  'Maka  aura  that  tha  drain  la  cbead  bafora  you  procaad' 

702  DELOISP  -  0.001 
706  NNN  -  400 


730  GOSUB  1200 
760  FOR  IK-  OTO  6 
760  GOSUB  1620 

770  DISPdKI  -  OISPIIKI  -  DISPINITIIK) 

7B0  NEXT  IK 

782  DISPX  -  OISPIOI  -f  DISP(1| 

7B4  DISPY  -  DISP(2I  DISP(3) 

790  DISPZ  -  DISP(4|  ■¥  OISPISI 

S30  W  DISPX  -  DEIOISP  >  0  THEN  OUT  BH320,23:F0A  J-1  TO  6;NEXT  J:  OUT  BH320.21 
S32  F  DISPX  DELDISP  <  0  THEN  OUT  4H320.20:F0R  J-  1  TO  NNNiNEXT  J:OUT  &H320.21 
B34  r  OISPY  •  DELDISP  >  0  THEN  OUT  liH320,20:FOR  J-1  TO  6:NEXT  J:  OUT  &H320,21 
B36  F  DISPY  +  DELDISP  <  0  THEN  OUT  4H320,17:F0R  J-1  TO  NNN:  NEXT  J:  OUT  &H320,21 
S42  F  DISPZ  -  DELDISP  >  0  THEN  OUT  BH320,63:F0R  J-1  TO  20:NEXT  J:  OUT  &H320.21 
S44  F  DISPZ  -f  DELDISP  <  0  THEN  OUT  4H320.6:F0R  J  -  1  TO  120;  NEXT  J:  OUT  BH320,21 
860  ' 

860  ‘FOR  K  -  1  TO  30000:  NEXT  K 
870  FOR  IK  -0  TO  7 
880  GOSUB  1620 

890  DISPIIK)  -  DISPdKI  -  DISPINlTdKI 

900  NEXT  IK 

910  FOR  IK  -  OTO  3 

920  GOSUB  1320 

940  NEXT  IK 

950  LOCATE  5,1:  PRINT  USING  '  #«.#«««':  OISPIOI:  OISPIH;  DISPI2I:  DISPI3):  DISP(4):  OISPISI:  DISP(6):  0ISP(7| 

960  LOCATE  13,1;PRINT  USING  "  itt.t  P(0|:  P(1|:  P(2I:  PI31 

980  N  -  N  >  1 

990  F  L  -  1  GOTO  1010 

1000  IF  I  N  MOO  LI  <>  1  GOTO  750 

1010  LSETXF*  =  MKSKPIOII 

1020  LSETYF*  -  MKS»(P(1I| 

1030LSETZFI  -  MKS«(P(2II 
1040  LSET  BF«  -  MKS«(P(3II 
1050  LSETX1F*  -  MKSHOISPIOII 
1060  LSETX2FI  -  MKSHOISPdll 
1070  LSETY1F1  -  MKS4|DISP(2II 
1080  LSETY2F*  -  MKS«(DISP(3II 
1090  LSETZ1F1  -  MKS«(0ISP|4|| 

1100  LSET  Z2F«  -  MKSIIDISPISII 
1110  LSET  Y21F*  -  MKS«IDISP(6|| 

1120  LSET  Y22F»  -  MKSt(0ISP(7ll 
1130  LSETTIMEF*  »  TIME» 

1140  M  -  M  1 

1150  PUT  »2,M 

1160  LOCATE  17,28:  PRINT  M 

1170  GOTO  750 

1180  END 

1190  ' 

1 200  CLS 

1210  LOCATE  1,1:  PRINT  ‘DISPLACEMENTS  IN  INCHES:’ 

1220  LOCATE  3,1:  PRINT  TABISI:  ’0X1*:  TAB(13|:  ’DXZ":  TABI21I:  ’DYI’:  TAB(29I;  ’DYZ’;  TAB(37I:  ’DZI’;  TAB(45I;  ■DZ2’, 
TABI53I:  ■DY21’;  TABISI  I:  •DY22’ 

1230  LOCATE  9,1:PRINT  'PRESSURES  IN  PSI:* 

1240  LOCATE  11,1:PRINT  ■SiomaX’:TAB|12l:’Slgm«Y’:TABI23l:’Sigm«Z’:TABI31l:’Pore  Preitura’ 

1250  LOCATE  17,1:PRINT  ’Number  of  Data  Recordad:  ' 

1260  LOCATE  19,1:PRINT  ’Praaa  F2  to  quit’ 

1270  LOCATE  20,1  :  PRINT  'Praaa  F3  to  ehatiga  racording  rata* 

1290  LOCATE  23,1:  PRINT  USING  '1  In  «««  raadinga  ara  baing  racordad.':L 
1300  RETURN 
1310  ' 

1320  ‘SUBROUTINE  PRESSURE 

1 330  LTKIOl  - 1 :  LTKI1 1  - 1 :  MOK  - 1  ‘Locka  0ASH8  on  channal  1 
1340  CALL  OASHS  IMOK,  LTKIOl,  FLAGKI 

1350  F  FLAGKOO  THEN  PRINT  'Error  in  aattbig  tho  cftannal':  END 
1360  MOK-14 

1370  CALL  DASH8  IMDK,  IK,  FLAGKI 
1380  MDK-5 

1390  TRANKIOI-VARPTRITEMPKIOII 
1400  TRANKI1I  -  40 

1410  CALL  DASH8  IMOK,  TRANKIOI,  FLAGKI 
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1420 

1430 

1440 

1460 

1460 

1470 

1460 

1400 

1600 

1610 

1620 

1630 

1640 

1660 

1660 

1670 

1680 

1690 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1830 

1840 

1860 

1860 

1870 

1880 

1890 

1980 

1990 

2000 

2030 

2040 

2060 

2060 

2070 


X-0  :  INDEX  -  0 
FOR  Ja>1  TO  38 

F  ABSrrEMP»(J)-TEMP%(J-1)l  <2  AND  ABS(TEMP%U)-TEMP%(J-f  1)|<2  THEN  X>X-fTEMP%|J);  INDEX  «  INDEX  -t-l 
NEXT  J 

VOLT  -  X/IINOEX* 20.48)  'In  m.Volt 
P<l%l  >  VOLT*CALTRAN8<l%l  'In  pti 
P(l%)  >PII«l  -  PINmAL|l%l 
RETURN 


'SUBROUTINE  DISPLACEMENT 

LT%(0l-0:  LT%(1|>0:  M0%-1  'Locks 0ASH8 on climM 0 
CALL  DASH8  (MO%.  LT%(0),  FLAGH) 

r  FLAQ%  <  >0  THEN  PRMT  "Error  In  sotting  ths  channol":  END 
MO%-14 

CALL  OASH8  |MD«.  I«.  FLAGTLI 
MOK-6 

TRANHIOl  «  VARPTR(TEMP«iOI| 

TRAN%(1I  -  40 

CALL  DASH8  (MO%.  TRAN%|0|,  FLAa«) 

X-0:  MDEX  -  0 
FOR  J-1  TO  38 

IF  ABS(TEMP«|J) -TEMP%(J-1||  <2  AND  ABS(TEMP%|J|  ■  TEMPTLIJ-r  11)  <2  THEN  X-X-fTEMP%|J|:  INDEX  -  INDEX  -fl 
NEXT  J 

VOLT  -  X/(204.8*IN0EX|  'In  Votts 
V  ABS(VOLT)  >  16  GOTO  1630 
DISPIIIi)  -  VOLT/CALLVOT|l%)  ‘In  inclios 
RETURN 

'SUBROUTINE  TO  FRESHEN  THE  SCREEN  |F1| 

GOSUB  1 200 
RETURN 


‘SUBROUTINE  TO  CHANGE  RECORDING  RATE  (KEY  F3) 
CLS 

INPUT  "Whst  is  ttw  rtcordino  rsts  le.g.  1  in  L)  L*:L 

GOSUB  1710 

RETURN 

‘SUBROUTINE  TO  TERMINATE  THE  RUN  (KEY  F2I 
CLS 

CLOSE  «2 

LOCATE  11,10;PRINT  "PRESS  ANY  KEY  WHEN  FMISHED" 

F  INKEYt  -""  THEN  2040 

OUT  SH330.0 

CLS 

END 


# 


Ml  .BAS  to  Porfomi  Numorieai  kitagrathMi  by  Mothod  I 


DECLARE  SUB  integ  (} 

COMMON  SHARED  beta#,  alpha#,  epsodot#,  D#,  epso#0,  con#,  n#,  P#0,  clp#0 
COMMON  SHARED  low#,  upp#,  stp#,  eps#0,  epsdot#(},  sumd,  feps#(},  q#(} 
CALL  inputdat 
CALL  integ 


SUB  inputdat 
beta#  =  .1 
alpha#  =  .01 
epsodot#  =  .1 
D#  =  (1  /  500) 
sm#  =  0 
n#  =  1000 
END  SUB 

SUB  integ 

DIM  epso#(201),  eps#(3),  epsdot#(3),  feps#(3),  P#{201),  q#(201),  dp#(201) 

DIM  nr#(3),  dr#(3) 

CLS 

INPUT  "Output  File  Name  a",  outS 
OPEN  outs  FOR  OUTPUT  AS  #2 
stpl#  s  (.1  •  IE-10)  /  201 
FOR  kX  =  1  TO  201 
epso#(kX)  a  (kX)  *  stpl# 

con#  »  (epsodot#  •  <1  ♦  epso#(kX)))  /  (epso#(icX)  *  (2  ♦  epso#(kX))) 
low#  *  IE-10 
IF  kX  >  1  THEM 
low#  s  epso#(kX  -  1) 

END  IF 

upp#  z  epso#(kX) 

stp#  =  (upp#  -  low#)  /  n# 

FOR  jX  »  1  TO  n#  /  2 
FOR  iX  »  1  TO  3 

eps#(iX)  *  low#  ♦  (iX  ■  1)  •  stp# 

epsdot#(iX)  ■  con#  •  (eps#(iX)  •  <2  ♦  eps#(iX))  /  O  *  eps#(iX))) 
nr#(iX)  «  (1  ♦  beta#  •  (LOG((epsdot«(iX))  /  alpha#)  /  LOG(IO))) 

IF  epsdot«(iX)  <  .001  THEN 
nr«(iX)  «  (1  ♦  beta#  *  (-1)} 

EM)  IF 

dr#(iX)  «  ((1  *  epe#(IX})  •  (2  >  epe#(iX)}  •  (.09  *  sps#({X))) 
feps«(iX)  >  nr«(iX)  /  dr«(iX) 

NEXT  iX 

SUR#  «  sunr  ♦  (((epe«(3)  -  epe#(1))  /  6)  •  (fepe#(1)  S  *  fepa«(2)  *  fepe#(3))) 
low#  ■  low#  2  *  stp# 

NEXT  jX 

PRINT  “low*",  low# 

P«(kX)  «  sun# 

NEXT  kX 
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'Nunerical  Oifferntiation 

FOR  IX  =  1  TO  201  STEP  3 

H#  =  (epso#(lX  *  2)  '  epso#(lX)} 

dp#(lX)  s  (-3  *  P#(IX)  ♦  4  *  P#(IX  ♦  1)  -  P#(lX  2))  /  H# 

dpOdX  +  1)  .  (P#(IX  2)  -  P#(IX))  /  u# 

dp#(lX  +  2)  *  (P#(IX)  •  4-*  P#(IX  ♦  1)  ♦  3  *  p#(ix  ♦  2))  /  H# 

NEXT  IX 

FOR  kX  =  1  TO  201 

q#(kX)  =  epso#(kX)  *  (1  +  epso#(kX))  *  <2  ♦  epso#(kX)>  *  (dp#(kX)) 
NEXT  kX 

FOR  mX  =  1  TO  201 

PRINT  «,  epsoOdnX);  P#(inX);  q#(niX) 

CLOSE  #2 
END  SUB 


M2.BAS  to  Rirfomi  Numorieal  Intogration  by  Mthod  H 

DECLARE  SUB  inputdat  () 

DECLARE  SUB  integ  () 

CONHON  SHARED  beta#,  alpha#,  cpsodot#,  D#,  epso#(},  con#,  n#,  P#(),  dp#() 
COMMON  SHARED  lou#,  upp#,  stp#,  aps#(),  apsdot#(),  tunB,  feps#(),  q#() 
CALL  inputdat 
CALL  integ 


SUB  inputdat 
beta#  >  .1 
alpha#  *  .01 
epsodot#  «  .1 
D#  =  (1  /  500) 
'sum#  =  0 
n#  »  1000 
END  SUB 


SUB  integ 

DIM  epso«(201),  eps#(3},  epsdot#(3),  feps#(3},  P#(201),  q#(201).  dp«(201) 

DIM  nr#(3),  dr#(3) 

CLS 

INPUT  “Output  File  Name  out$ 

OPEN  outs  FOR  OUTPUT  AS  #2 
stpl#  «  (.1  -  IE-10)  /  201 
FOR  k%  »  1  TO  201 

f 

suntf  =  0 

epso#(kX)  s  (k%)  •  stpl# 

con#  =  (epsodot#  *  (1  epso«(kX)))  /  (epso#(kX)  •  (2  epso#(kZ))) 
lOM#  =  IE-10 
upp#  3  epso#(kX) 

PRINT  USING  "low  >  ##.##r<;  low# 

PRINT  USING  “upp  s  ##.###";  upp# 
stp#  *  (upp#  -  lou#)  /  n# 

FOR  jX  =  1  TO  n#  /  2 
FOR  iX  =  1  TO  3 

eps#(iX)  >  low#  (iX  -  1)  *  stp# 

epsdot#(iX)  »  con#  *  (eps#(iX)  •  <2  +  epa#({X))  /  (1  ♦  eps#(iX))) 
nr#(iX)  »  (1  ♦  beta#  •  (LOG((epsdot#( iX))  /  alpha#)  /  LOC(IO))) 

'IF  epsdot#(iX)  <  .001  THEN 
'nr#(iX)  •  (1  ♦  beta#  *  (-1)) 

'END  IF 

dr#(iX)  »  ((1  ♦  eps#(IX))  •  (2  ♦  #pa«(fX))  •  (0#  ♦  eps#{iX))) 
feps«(iX)  X  nr«(iX)  /  dr«(iX) 

NEXT  iX 

sue#  -  sui#  ♦  (((ept#{3)  -  eps#(1))  /  6)  *  (feps#{1)  ♦A  *  feps#(2)  ♦  feps#(3))) 
low#  X  low#  *  2  *  stp# 


'PRINT  “lows",  low# 

P#(kX}  >  sun# 

NEXT  kX 

'Numerical  Differntiation 

FOR  IX  =  1  TO  20T  STEP  3 

u#  3  (epso#(lX  *  Z)  -  epso#(lX)) 

dp#(lX)  =  (-3  *  P#ClX)  +  4  *  P#(IX  ♦  1)  -  P#(IX  ♦  2))  /  w# 

dp#<tX  ♦  1)  «  (P#(lX  +  2)  -  P#(IX)}  /  u# 

dp#(lX  ♦  2)  »  {P#(IX)  -  4  *  P#<IX  +  1)  +  3  •  P#(IX  ♦  2>)  /  w# 

NEXT  IX 

FOR  kX  »  1  TO  201 

q#(kX)  *  epso#(kX)  *  (1  ^  epso#<kX))  *  (2  ♦  epso#(kX))  *  (dp#(kX)) 
NEXT  kX 

FOR  mX  =  1  TO  201 

PRINT  #2,  epso#<na);  P#<inX);  q#(nX) 

NEXT  nX 
CLOSE  #2 
ENO  SUB 


